i;as.\ 0(32294 


TECHNICAL REPORT 


rNASA-CB-13229) SOHLIHEAB AHD DIGITAL 
EAN“HACHIHE COMTBOL SYSTEMS HODELIHG 
Final Bepoit (City Oniv. of New *<>**•> 

112 p HC $7.75 CSCL 05H 


K73-29060| 

Unclas | 
G3/05 13457 ) 


NONLINEAR AND DIGITAL 
MAN-MACHINE CONTROL 
SYSTEMS MODELING 



Langley Research Center 
National Aeronautics ana S 
Hampton. Virginia 23365 


Administration 


FINAL REPORT No. 72-447-01 
Grant No. NGR-33-0So-0c3 
November 1972 



THE CITY COLLEGE of 

THE CITY UNIVERSITY of NEW YORK 



AHftN0WLqDGM3NT 


Tha author wishes to express his appreciation to 
Mr. James J. Adams of Langley Research Center for his useful 
suggestions throughout the duration of this study.. Thanks 
are due to COOT students Messrs. Theodore Rapban and Gideon 
Cialan for their help in obtaining digital computer results. 


1 



TABLE OF CONTENT S 


1, INTRODUCTION* 

2, STATEMENT OF PROBLEM . * * ' 

3, HUMAN OPERATOR MODELING - CONTROLLABLE FORM . 

3.1 Introduction 

3.2 Model-Reference Error Equation ...» 

3.3 Formulation of a Liapunov ^unction * « 

3.4 Adaptive Identification Controller . « 

3.5 Nonlinear Mathematical Model 

3.6 Linear Mathematical Model 

4 , MODIFIED MODELING TECHNIQUE - CONTROLLABLE FORM 

5, HUMAN' OPERATOR MODELING - OBSERVABLE FORM . . . 

5.1 Introduction 

5 ,o Model-Reference Error Equation 

5.3 Formulation of a Liapunov Function 

5.4 Formulation of Controller and Matber.ati.cal Modal 
(Nonlinear Case). •♦»•**** f 1 * ' * * * 

5.5 Linear Mathematical Model ....»* 


o • • • 


» t » • 


6 . EXAMPLES 

6.1 Example l.s Identification of Two Parameters . . . 

6.2 Example 2.: Identification of Four Parameters. . . 

7. CONCLUSIONS AND RECOMMENDATIONS 


7.1 Conclusions . • 

7.2 Recommendations 


**•••*• 


1 

4 

6 

6 

6 

8 

10 

11 

13 

16 

22 

22 

24 

26 

27 

31 

35 

35 

41 

4S 

43 

51 


ii 



06 

Lb 

9 L 

69 

£9 

9S 

£5 


, . • 1 Z*Z1 

. . * X ©xdcrex£. 1*21 

SXTflSEH VE1MH00 TVXIOIQ *2t 
,,«.*••• X^HTiOId 
.***•*• 9 XIGKSduV *0X 

V XlGKSdaV *6 

.*•••*• sSDK'dV.SiSH *3 


&Ss<i 




FijLs 


Title 

Page 

1. 

Vo , J{ o . Reference System Configuration with Hunan 
Operator' s Response Applied to Controller 

70 

2* 

Model - Reference System Configuration with _ . 

Mathematical Model’s Response «ppUau to -on.* 

71 

3. 

Logic Tlow Die gran - Nonlinear Model in 

Controllable 

72 

u * 

Logic 

Fora 

Flow Diagram - Linear Model in Controllable _ 



73 

5. 

Logic 

Fon 

-e- Die nr an _ nonlinear Modal in Observable 

7° . . : 

74 

6 * 

Ferns 

*j r «a r v Ci i^ in Observable 

T": n'-r Cl A ^T<i7d. — jjlr.»3r 

7 . : 

75 


iv 



abstract 


This study examines the utility of an adaptive modeling technique 
by which controllers can be synthesized to provide corrective dynar-ics 
to a human operator's mathematical nodal in closed loop control system. 
The technique utilizes a class of Liapunov functions formulated for this 
puroose, Liapunov's stability criterion and a model-reference system 
configuration. Toe Liapunov function is formulated to possess variable 
characteristics to tales into consideration the Identification dynamics. 
The time derivative of the Liapunov function generates the ider.tif icatio 

and control laws for the mathematical model system. These laws penait 

- - . + r>r>,iia v ~ updates the hirr^n cc>3 rater’s 

the realisation of a controller rf.-ic.* 

,, + < _, ^^-sters so that model and human operator produce 

ma t ner.e t r. ca j- _ y*-* 3 — ww 

th# same response when subjected to the same stimulus. 

A -.■ary ’useful feature of this study is the development of a 
digital computer program which is easily implemented and modified con- 
current with experimentation. The program permits the modeling process 
to interact with the experimentation process in a mutually beneficial 

way* 


v 
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1, INTRODUCTION 

This study present* an adaptive modeling technique which one 
nay utilise to identify and control a h^an operator's wtheaatieal 
node! . The technique is based u P on a nodel-reference system configu- 
ration, as shown in ?lgs. 1 5c 2, a Liapunov function formulated for 
this nodel-reference system and Liapunov's stability criterion. The 
reference system in the configuration represents the actual human 
operator. The model system in the configuration is constructed Ini- 

, t t * cel representation of th3 reference 

tially es an approximate n- - — - y 

system and is derived from consideration of the physiological processes 

evident in the actual human operator when performing a given tracking 

c k formulated to possess variable chars o 

task. The Liapunov function is io*....uiat-u * y 

^ the -*. 3 ui«i identification dynamics. The 

teristics in order t.o cerxve tna 

, r. on -ene rates the identification 

derivative of tr.a ^ispu.o / =>*'•* 

and control law. for the mathematical model. These identification and 

control laws permit the realization of a controller which updates the 

mathematical nodal parameters so that model and h«an operator give the 

„ cvhia-tpd to the same input sign 3 -* 

sans response wnen suoje^'-e^ ^ — 

. t _ y»q a theoretically consistent 
T’ne identification teennirue Oxxers a J 

. u i .j + ! pe fron s xoe r utertsl data. 

. human operator activities iru 

aopreve cn iox reme 11 -.5 *- J — a y~ 

That is, the approach is theoretically capable of operating on the 
available experimental data and identifying a model system that can be 

used in place of the human operator. 

A very useful feature of this study is the development of a dig- 

, , , 1 - p3< .n Y imuleaiented and modified concurr-n^ 
ital computer program whicn is easily irapi-menw- 

• ♦ r n this way, the modeling process interacts with 

vith. experimentation. in j» 
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4 f a * i nn oroces c in a mutually benefical way. 
tha ex:>sri^ 3 nta *ion proce^- 

. u , „ .... , tu 4 v that the Dueling technique con- 

It was the intent of trd= stu-y 

■ •• unities 'U order for the technique to be versatile 

tain "learning" capabilities .no.- 

ooe^tor activities 

enough to be applies to as broil a range • 

. T^is ’‘learning" feature 
as possible, including nonlinea* - 

. . iva qAohtivo controller structure* 

is incorporated into ohs a-a,uu 

a. . +v^= ^od^Tin^ technique is to synthe- 

Cns of the basic tasks of thi* n - 

j viatical model representing tne 

size an adaptive controller an. - 

4a , ~^-ence system. The realization of such 
human operator m *-.js no^i— 

^ .v, ^sir^d choice of the state variab.es. 
a model system de penis upon 

, ^ 0 -c33 which lead to different realizations. 

There exist several uoe.\~ 

p n< *- 'w.n vrn-.o n operator mati-^— 

. . n r - 3 ^ v. zaticns for 

In this study two 0 - 

. i r- r=t is the controllable realization 

w-a t '• ca 1 model is preser.o^, - - -- - 

^ observable canonical realizs- 

, . . j t-i-q cacOtll .o '■•••*' ^ w w ' 

of tie nc os j. sys -r.;j “" w 

. r- n ^r a -Va~es and disadvantages of both 

tion of the model system. — — — ° 

.. , . i V ;„ r*=> '~> 0 "+. . Particular attention 

^ ,^ 9 H 0 r S 3 -e describe m *m=. re i '°--- 

these realise lIg.io - 

_ *^ a of the identifying dynamics. It is 

is given to the conve* v 

, *. . * Gratification time for the 

. , <r , + v.A g study that tr.s pars-e^e* 

c ho A n i x n t jiio o L- •A-s-j 

v , caused nodal is nuch lass than the parameter 

observable ca^o-i- 

, C 0 - t v 9 controllable realization. The cis- 

identification tine requ^r^d - 

,, ^K, 9 t‘o'' is that one cannot directly 

advantage of the observable re-liza..o.. l 

. » _; s _t 3S except the first state. This be.-Cu-- 

neasure all of the s.a^s -* ,j ~~ 

aooerent in the report. 

* consist* o" two basic parts. The first part of the 
This report consists 

n i .-fvtai^sy a o crouches which lead to a 
report discusses two theoretical r.ouel- ^ ap„. 



controllable realization and to an observable canonical realization 
of the human operator mathematical model. The second part presents 
examples and experimental evidence to substantiate the validity of the 
mathematical results. In this second part of the report the reference 
system is represented by reels of tape containing recorded digital 
innut-output data of an artificial human operator when performing 
a given tracking task. A digital center is used to generate the 

identification dynamics for the mathematical Eodsl * 

The report is organized in the following manner. First, a nathe 

natical description of the problem is presented, illustrating the for- 
mulation of the model-reference system and the error differential 
equation. Se*t, the formulation of a Liapunov function and its time 
derivative is given. Imposing constraints on the Liapunov function 
and its time derivative permits one to obtain the identification and 

control laws which realize an adaptive controller. After establishing 

,„ho-' 3 the mathematical nodal reoresen 
a relations hip among several equations 

the human operator is derived. 
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2. STATIST OF P303L3K 

Before considering the derivation of the modeling approaches, 
a statement of the basic problem is given. Consider the model-re fere nee 
system configurations depicted in Figs. 1 & 2. The reference system 
may represent the actual human operator* The stimulus -response data 
for the human operator are available from experimentation* It is as- 
sured that the human operator when performing a control task may be 
described by a vector difisrential equation 

2 = A z + Hr (l) 

v Yat 9 z denotes the response state vector, z = dz/dt and r represents 
the stimulus vector. The square matrices A and H are unknown. The 

J s to find a set of elements for matrices A and H so that the 
f t g * ' p 11 -atones very closely the response of the numan 
operator when both Tq* (1) and human operator are subjected to the same 

stimulus . 

As a starting point for the modeling approaches, one may choose 
a tentative mathematical model system of the same form as Fq* (1). Let 
tbs model reoreser.t ing the controllable realization be described by 

u 

x = Bx f C_r ^ 


where r is the same stimulus state vector as for Fq* (1), x denotes 
the mathematical model’s response vector and x r dx/dt . The square 
matrices 3 and C are also unknown except that one may assign a set of 
initial values to the elements of these matrices. Let these initial 
matrices be denoted by 3 q & C Q . The initial values must be chosen to 



+ontativQ mathematical model system 

be ««. «1» ”«*“ “ f , „ ,,„ ur . 

. The problem is to idenxny 

as defined by a formulated Liapunov function. 

„ t . t ftt tha end of the Identification interval 
matrices B and C so that at the enn 

(B B) _Aand(C 0 -C) — H. The mathematical mode sys 

colidered identified and representing the t-n operator ^ 

There are several approaches one may take in order o 

.. Thi3 s tudv considers two approaches and points out 

modeling problem. This stuay JA ,., nn the 

, e either approach, in addi » 
the advantages and disadvantage ■ suitohi ng from one 

atudy discusses tte trade offs one must consider ^ switching 

approach to tha otter. . , to 

The next section presents the modeling technique v 

. »«— i — -« « — - r “- rint ; “tn 

.bb.l i. 4.^.1, ». »*■«“* *° • U "** r "° i, ‘ ”* ■' 



6 


3. HUMAN OPERATOR MODELING - CONTROLLABLE FORM 


3.1 Introduction 

There are several methods for realising relationships such as 
given by Eq. (2). The realisation depends upon the particular choice 
of the state variables. The state variables chosen in this section are 
referred to as the phase variables, a name that stems frcn the coordi- 
nates of the phase space. This choice of state variables is a natural 
one for the engineer, as these variables have a ready physical inter- 
pretation. For example, for the human operator the output (response) 
may be chosen as the first state variable, the rate of change with 
respect to time of the first variable then is chosen as the second 
state variable, the rate of change with respect to tine of the seccna 
state variable is further chosen as the third state variable etc. For 
an n-th order system n state variables necessary and sufficient to 
represent the dynamic behavior of the system. Note that the second and 
third state variables represent the velocity and acceleration of the 
human operator’s response. 

In order to derive the identification dynamics for the control- 
lable form of the human operator’s mathematical model system depicted 
in Fig. 1 or Fig. 2, one must first formulate the error vector differ- 
ential equation for the model-reference system in terms of these phase 
variables * 

3 • 2 Model-Refere n ce Error Equati on 


It is obvious from Fig. 1. that initially the response of the 
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, s „, n not be the same as the response - 
human operator (reference system) will 

„ r »«-«- —1. ““ *" ! ” bi *'*' a “ ‘ 

. This i3 dus to the fact that the tentative model’s 

same stimulus. This is cue «.« 

parameters are initially cashed. As long as these parameters s» 

J «. — t, «-» «• -“•* “* ‘ " “ 

* \ „ er ,™ S e is the model-reference system 

operator's (reference system) response 


(3) 


e = x - z 


, Differentiating Sq* (3) 

Note that e, _x and j are time varying vectors. 

with respect to time yields 


U) 


e s x - a 


,* do. denotes «. .i» 

(2) ,nto «) Tleld*. .« 

, ^ for +b© model-reference system 

error differential equation for 

T T ( 5 ) 

e = Be + bu z * dw r 


where 


bu 


l T = (B - a) 


dw T s (C - H) 


( 6 ) 


(7) 


,„t. superscript t ds.of. th. tr.»P— ™= <5 > “ 

„„,.,d clop • P*~ — “* 4 *“ 4,4 ‘" ’ 
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system. The choice for this type of model representation was governed 
by the desire to obtain a controllable realization of the final mathe- 
matical model system. Note' also that Eqs. (6) & (7) contain the iden- 
tification dynamics for the mathematical model's parameters. It should 
be pointed out that other realizations of the model system are possible. 

A different realization, based on a different choice of state variables, 

is presented later in this report. 

Note that vectors e, u and v denote the model-reference system 

error and parameter misalignment vectors. One may now formulate a 
Liapunov function and its time derivative that incorporate the error 
and these misalignments. Before one proceeds with the formulation of 
a Liapunov function and its time derivative it is useful at thi. point y 
to become familiar with the detailed forms of matrices and vectors used 
in this section. This information is presented in Appendix A. 

3,3 Fonnul r tion of G, L iajsu r i p v F ^c Uon 

Equation (5) Hiay be viewed as consisting of three perturbational 
vectors, namely e, u and _w. An appropriate Liapunov function should be 
positive definite in the error as well as in the parameter misalignments 
Therefore, one may choose a Liapunov function of the form 

t T T (8) 

Y — e Me u Nu y Qlf 


where matrices M, N and Q are sy^etric square matrices. The elements 
of these matrices may be constants, time varying and/or functions of the 
state variables (3) . '-men considering a linear mathematical model of the 



human operator one may restrict matrices M, N and Q to consist 
constant elements only* 

Differentiating "q. (8) with respect to tine and then substituting 
Sq* (5) and its transpose given by 


e T s q T 3 T f z T ub T * r T w d T 


(9) 


one may obtain 


the time derivative of the Liapunov function in the form 


V = - e D e 4- 2 


u T M t- iu X N z'(b‘Ke) 


T- T,,T, 


U V 


* 2 


w Q *• 


to Q 


T T 

* r (d Mg) 


v 


( 10 ) 


where 


D x = (B«l KB * A) 


UD 


and the dot denotes the time derivative. 

Liapunov's criterion for stability calls for V ? 0 and V f 0. 

One way to comply with Liapunov's criterion for stability is to constrain 
the elements of the D x matrix, denoted by and d ij , to satisfy the 


conditions 


and 


d u > 0 


d lj * d jl = 0 


( 12 ) 

( 13 ) 


where i and j denote row 


and column respectively and let 
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v T Q * 4w T Q 4= r T (d T hfe)J = 0 ( 15 )' 

The conditions given by Sqs* (12) & (13) enable one to evaluate the 
elements of D 1 and M matrices^. At this point one must perform a test 
to insure the positive definitnsss of these matrices* Matrices N and Q 
must also satisfy the conditions N y 0, Q 0 and are used as free 

l 

design parameters in order to influence the convergence time of the 
identification dynamics* Equations (14) & (15) csosssti tafc* the basic 
equations from which the identification dynamics and the adaptive con- 
troller are realised. 

, * 

It should b 3 noted that the resulting V is negative semidefinite 
because it depends only upon the model-reference system error and does 
net depend on the vectors u and This suggests that one should expect 
scene oscillations in the model parameters at the end of the identifica- 
tion interval, even if the model-reference sy3ten error e is zero. 

3*4 Adaptive Id e ntification Controller 

As mentioned before, Eqs* (14) & (15) give rise to the identifica- 
tion dynamics. In order to realize the adaptive identification controller 

T T 

one has to rearrange Eqs . (14) & (15). Since (b Me) and (d hie) are scalar 
quantities, one may obtain fresn Sqs* (14) & (15) the following relation- 
ships 

U T = - 4u T !«' 1 - /N _1 (b T He) (16) 

w T = - tw T QQ _1 - r T Q _1 ( d T Ke ) 


(17) 



- n 


The transpose 


of Sqs. (16) & U7) 7^ lds 
■ = . ^- 1 T H T u - H- 1T x(b T K-) 


(18) 


_ _ ig _1T Q T w - Q* 1T l(dV.e) 


(19) 




S,„,lon. (13) * <M> "f""* -Tt 

H u ecmputed « *. CU, - > - « ~ l 

, . w m an H Q must be positive ae- 

, . _ w n +e that matrices M, *1 ana * 

fication dynamics, hote 

j n +he Liapunov function, 
finite Since they appear in the hiap 

3.5 Menu near Ha tte me^icajJfoM 

, 4 iwl -od“l system one must find 
in order to formulate the rathe.atrcal ,..o--l 

il\ (7) (16) & (17). Asswing that the 

. relationship ar.ong tqs. (6), (7), U . .... 

V slower than the identification 
changes in the human operator are much slo 

one -ay consider matrices A and 
, . for the model's parameters, one ...aj 

time required f -- thp 1<tert .« cation interval. Tterefore, 

H as being time-invariant during the identic* 

„ ,ts & ( 7 ) .with respect to time yields 
differentiating Sqs. W * 1 " " 

( 20 ) 

* »T 
3 s bu 


and 


* .-T 

C : dv 


(a) 


nee b and d are constant vectors 
s. (20) & (21) ona ° btains 


Substituting 3qs. (lb) & (l 7 ^ int0 
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t. -1 T -1, T ^ 
3 - _ ^bu NN - bz N (b Ke • 


C = - Hv QQ” 1 “ dr T Q“* 1 (d T Me) 


„ ^ e f00 \ p. OV 0 r the identification time interval 

Integration of iqs . 122) St Ud) 0V9r 


yields 


t 

r B 0 {-buh" 1 +■ bzVc/Ke) dt 


C = C n - id 


4 - dr T Q“ 1 (d i Me)J dt 


. « u ,> ii„ o coiled ci ces for the mathematical 

where 3 0 and C 0 are the initially assure 

B oael. These oatrioas are giver, in Appendix A. The integrals in 

Sq s. (24) & (25) represent the nonlinear identification dynamics for 

J D , p v, ntp th a t these integrals consist of parameter mis- 

matrices B and b« .^oue xnao 5 

. as well as the error and 

alignments, model-reference input a..d re- r o..se as wax 

elements of the Liapunov function. 

To complete the formulation of the mathematical model one substi- 

tutes Sqs. (24) & (25) into Sq. (2). This yields 


“ r *j /n 

*-bu T NN~" 4- bz'ir^yMe) 
c - - -- J 


(26) 
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Equation (26) represents the h^an operator-s nonlinear mathematical 
model with its identification dynamics. Note that the identification 
dynamics depend upon the human operator-s response and stimulus. 

Since Sqs • (24) , (25) , (26) ( 10 > are ^ inear > time V3r7lng 

and interrelated, then the difficulties that arise in their computation 

may he alleviated by solving these equations iteratively using a digital 

computer. The digital computer program is easily implemented and ca 

modified concurrent with experimentation. This permits the identifica- 
tion process to interact with tta experimental data in a mutually bene- 

ficial way* 

The mathematical model is considered identified (representing 

N h 3=0. Due to the senidefinite V, the model's 
the hucnaa operator) //hen _ 

. _■> With this modeling 

parameters oscillate about sane nominal valu-s. 

-■Mo.- th= nominal values as being the identified 
technique one may consider the nominal vai 

to. «■ 111 “* “■* 

has identified ..trie > n«d »■ *■*•«■““““ *" ,h °” 

proper oh.lo. of » .* « »««„, ia ““' 

, . This modeling technique is illustrated 

fication of the model parameters. This moaeii „ 

via an example where all the details are shown. 

3,6 Tdnflar Ma thsniatical Modal 

The equations derived in the previous sections may be simplified 
ute n one considers that the human operator may be viewed as a linear 
system For certain tracking tasks with appropriate conditions this may 
he the case, prefers, it is useful and beneficial to investigate the 
simplifications and the resulting equations leading to a linear 
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mathematical model. 

Let us then presuppose that the human operator (reference system) 
as shown in Fig. 1. performs in a linear fashion. Since a linear 
mathematical model of the reference system is desired, one may let 
matrices M, N and Q in the tentative Liapunov function (Eq. 8) to con- 
sist of constant elements. Therefore, 


M = 0 
N = 0 
Q = 0 


(27) 

(28) 
(29) 


In view of Eqs • (27), (28) & (29) one obtains from Eq. (10) the 
time derivative of the Liapunov function as 


V = - e De 4. 2 


r.T„ t, ' 

u K <■ s (b 


■1 

.T T T 

w Q *■ r (b lie) 

) U 4r 2 

J 

- J 


w 


(30) 


where 


D = (BM * MB) 


(31) 


From Eqs . (16) & (17) the simplified controller equations are 

u T = - zV^b T Ke) (32) 

w T = - rV^t^Me) 


Following the same procedure as in the previous section, matrices 
B and C are identified by the equations 
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B = 3 0 - / bz T N _1 (b T fto)dt 


C = C o - 


dr T Q _1 (d T lie)dt 


(34) 


(35) 


where B 0 and C 0 are the initially assumed matrices for the mathematical 
model (see Appendix A) and the integrals represent the linear identifi- 
cation dynamics for matrices B and C. To formulate the linear mathemati- 
cal model one substitutes Eqs. (34) & (35) into 3q. (2). This gives 


x - 


B rt - / bz^V" 1 ( b T Me ) dt 


x + 


+ C n - 


TIT 
dr Q (d ife)dt 


(36) 


Equation (36) represents a linear mathematical modal of the human operator. 
Tte integrals in 3q. (36) are the identification dynamics. Although, 

Sq. (36) may be solved much easier than its nonlinear counterpart, many 
difficulties are still alleviated when one uses a digital computer. An 
example illustrating the identification of two parameters for a linear 
mathematical model is presented in the later part of this report. 
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4. MODIFIED MODELING TECHNIQUE - CONTROLLABLE FORM 

A review of the modeling technique and Eqs* (26) &. (36) show 
that the mathematical model of the human operator is a function of the 
M matrix in addition of being depended upon otter variables* But the M 
matrix which is computed either from Eq. (11) or Eq* (31) requires the 
knowledge of the B matrix. Note that the B matrix changes (getting 
updated) during the identification interval. This in turn, suggests that 
the M matrix will also change during the Identification interval* In 
order to take into consideration these changes, one must set up a digital 
computer subroutine which will compute at each iteration a new set of 

M and B matrices. Since the M matrix appears in the Liapunov function, 

« 

therefore one must also at each iteration interval monitor V and V so that 
Liapunov’s criterion is not violated. This requires another subroutine. 

In order to circumvent these additional computations , one may 
consider Fig. 2. and modify the model-reference error differential equation. 
From Eq* (3) one obtains 

z = x - e (37) 

Substitution of Eq. (37) into Eq. (5) yields a modified model-reference 
error differential equation 

e - il e + bu*x 4- dwyr (38) 

where bu^ and dw*^ are as defined by Eqs • (6) & (7). Note the modifica- 
tions in Eq. (38) when compared with Eq. (5). Equation (38) contains 
the human operator (referenco system) matrix A and the mathematical 
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models response x* 

Squation (38), as before, consists of three perturbational vectors 
namely e, u and Therefore, one may proceed with an identical develop- 
ment as previously. The Liapunov function is chosen as given by Bq. (8). 
The time derivative of the Liapunov function, which incorporates the 
error differential equation beccnea 



T 

e * 



4 . 2 ! 


*T 

W Q 4- 


T/ ,T. 


iw Q * X (d Me) 


v 


(40) 


where 

D, = (A T M <■ MA <- H) 


(a) 


Note that the dot denotes the time derivative and superscript T denotes 
the transpose. The nonlinear identification control laws obtained from 


Eq. (40) are 


*T 

u 



x T tr 1 (b T Me) ' 


(42) 


w T = - ^w T QQ - 1 - r T Q -1 (d 1 Ke) 


(43) 


Tte transpose of Sqs. (42) & (43) yields the controller equations 


u = - *r 1T H T _u - N -1T x(b T Ma) 


(44) 


jj s - 1T Q T v - Q- 1T r(d T Ms) 


(45) 
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Substituting Sqs. (42) & (43) into Sqs. (20) & (21) and integrating the 
resulting equations over the identification interval, one obtains 



c 


T -1 T 
bx N X (b Me) 


dt 


(46) 


t 


C = 



dr T Q"* 1 (d T Me) 


dt 


(47) 


whore B Q and C 0 are the initial value mtrices* Substitution of Sqs. 
(46) & (47) into Eq. (2) produces the nonlinear mathematical model of 

the human operator 



, T * -1 , T -l /v T N 

4bu NN bx i' (b ho ) 


dt 




*1 C 0 - 


-Ld^QQ’ 1 * dr T Q* 1 (d T Ka) 


dt V r 


(48) 


Note that the identification dynamics for this modified mathematical model 
depend upon the stimulus and the model's response as compared with the 
previous model which depends upon the stimulus and the human operator's 
response . 

As before, Sqs. (40 - 48) may be simplified if one desires to study 
the human operator as a linear system. This is possible for certain 
tracking tasks with certain conditions. It should be noted, however, that 



if there are inherent nonlinearities in the human operator that cannot . 
be separated, than one should use the nonlinear mthesatical ,odel to 
represent the hu^n operator. If the nonlinearities can be separated 

then one may simplify 3qs * (40 - 4B). . 

As shown in section 3.6, ( linear case) one say let K = 0, S = 0 

and Q = 0. In view of these simplifications Eqs. (40 - 48) beccoe 


V = - e *■ 2 


u T ti 4- /(b‘Mg) 


T,, T, 


U 4- 2 


w T Q + r^d'Mg) 


T, ,T, 


w 


(49) 


whara 

= (A T H * MA) 


(50) 


= . /rV*) 

f = - rWMe) 

IT T 

u r - N~ x(b Kg) 


w = - Q _1T r(d T Ka) 


z 

B a B 0 - ( b3^N _ ^(b T ‘i§) | * t 


(51) 

(52) 

(53) 

(54) 

(55) 


c = c 


i 

dr T Q* 1 (d T Ke)dt 


(56) 
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£ =|b 0 - j bx T M -1 (b T te)dt j x * 


t 

|c 0 - f dr T Q- 1 (d T lb)dt j r 


( 57 ) 


The problem arising in this modified modeling procedure is that 
matrix A. is unknown initially. Therefore, one cannot cccapute matrix M 
from Sq. (U) or Sq. (50) by satisfying the conditions given by Sqs . (12) 
and (13). Since matrix A represents the reference system which is 
actual Inman operator, then one may assume a priori that matrix A is 
stable. This will always be the case since the human operator is as- 
sumed to be stable when performing a tracking task. According to LaSalle 
and LefschetJ 181 if A represents a stable system, then there eari* 
a positive definite symmetric M matrix for every positive definite D 2 
matrix in Sq. (50). Therefore, in order to generate the identification 
dynamics , one may choose a set of elements for D 2 and M matrices as long 
as D 2 and K are positive definite. An alternate way would be to use the 
Chosen initial values of matrix B 0 for matrix A and then choose a matrix 
D 2 and solve for the elements of the M matrix. In this way, one obtains 
a preliminary set of constants which may be looked upon as weighting 
factors for the different controller identification loops. In either case 
this is not a systematic way of computing the elements of the M matrix, 
but it is more likely a trial and error approach. This is the disadvan- 
tage of the modified modeling technique. Work is being continued in this 
direction to establish a systematic way of computing the H matrix elements 
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for this modified modeling approach. Experimentation has shown meanwhile 
that an intelligent choice of D 2 and M matrices leads to very good iden- 
tification of the human operator's model parameters. Digital compute 
simulation studies have been performed using this modified modeling ap- 
proach and the results are included in this report. 

A review of the presented moieling technique (controllable form) 

indicates that the adaptive identification dynamics require the measure- 

. „ __ iro -i i a3 the entire stimulus vector* 
ment of the entire response vector as well 

This is due to the fact that it is desired to realize the mathematical 

model based upon a phase variable representation of Sq. (2). If pha30 

variables of the response and stimulus are not available from measurements 
than the identification dynamics are inadequate to update the model's 
parameters. One may generate the non-neasurable phase variables using 
differentiators, then the generated phase variables are either noisy or 
have discontinuities. This produces oscillations in the identified model 

parameters and it becomes difficult to determine when the model is iden- 

, . + p that the identified model parameters 

tified. Experimental results indica 

- i „ Thflrpfore in- this modeling tech- 
oscillate around some nominal valua * 

. , values as feeing the identified parame- 

nique one may consider these nominal values as being 

tens for the human operator in the reference system. 

To overcome this problem a different modeling technique is considered. 

This modeling technique generates identification dynamics which do not 
require the measurement of the entire stimulus and response vectors. All 

„ the actual stimulus and response. The improved model- 

ono has to measure is the actudj. 

ing technique is shown in the next section. 
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5. HUMAN OPERATOR MODELING - OBSERVABLE FORM 
5.X Introduction 

A 3 mentioned before, the desired realization of the mathematical 
model system for the hman operator depends upon the choice of the state 
variables. Let us define, in this section, state variables which differ 
from the phase variables presented in the previous sections. In order 
to generate tbs new set of state variables let the mathematical model be 
described by the following n th order differential equation 

dVdt n * d n - 1 (h Q . 1 >:)/dt n ' 1 * .... * d(hjx)/dt h 0 x = 


d n r/dt n * d rl ' 1 (g n _ 1 r)/dt n - 1 4. ... ♦ d( Bl r)/dt * g Q r 


= 8 , 


(58) 


where x denotes the response and r the stimulus. One may rearrange 
Eq. (58) so that the stimulus and response derivatives of the same order 

arc combined* This yields 


dVdt n * " g n-i r)/dt “"" * *“* 


n-1 


4 , dthjx - g x r)/dt 4- h o x = g o r 

The high order differential equation given by Sq. (59) is equivalent to 
E s-etem of n first order differential equations. Let the states be 


defined as 
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= x 

x 2 = ii* ( h n-l X " S n-l r) 

• r i (60) 

x 3 s x 2 * (h n-2 X ' g n-2 r) 

# 

= ^n-l * (b l X ' g l r) 

where the dot denotes the time derivative. Rearranging the set of 
equations given by Sq. (60) one obtains the state equations for the 

mathematical model 


x x = - ^n_i x * x 2 * g n-l r 
x 2 = - ^2 X * X 3 * g n-2 r 
x 3 = - V 3 x + x 4 g^ 3 r 


( 61 ) 


' h l x * *n 4 g l r 
Xn = - h Q x * g 0 r 


where x and r are scalar quantities and can be treasured. In matrix 
notations Eq» (6l) has the form 


x = Hx 4. gr 


( 62 ) 


where x, x and g are vectors, r is a s 


calar and H is a square matrix 
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The expanded forss of these vectors and matrix H are given in Appendix 3. 

As a starting point for this modeling approach, let the human 

operator be described b y 

j = Az «■ br (63) 

where z, z and b are vectors, r is a scalar and A is a square matrix. 

Note that Eq. (63) has the same form as Eq. (62). The problem is now to 
identify a matrix A and a vector b such that the response z of Eq. (63) 
matches the experimental output data of an actual human operator when 
performing a tracking task and is subjected to a known stimulus. In order 
to proceed with the identification of matrix A and vector b, one must 
first formulate the error differential equation for the composite model- 
reference system as shewn in Fig. 2. Then, using Liapunov's criterion, 
one must obtain the identification dynamics and a controller which will 
update matrix K and vector s , so that at the end of the identification 
interval the response x matches the experimental response data of the 
actual human operator. The mathematical model is considered identified 
when H — f £ ar ,d s _ b. The realization of Eq. (62) then represents the 
reference ays tern and is the mathematical representation of the actual 
human operator in ths observable form* 

5.2 Model-Refe^nca Error Equation 

Let the error describing the difference between the actual htman 
operator’s response and the tentative model system's response when both 
are subjected to the same stimulus be defined as 


e r z - X 


(64) 



25 


Differentiating Eq. ( 64 ) with respect to tine yields 

• • » 
e s z - x 


Substituting Eqs. (62) & (63) into Eq. (65) gives 


(65) 


e = Az «• br - Hx - gr 


( 66 ) 


The above error differential equation nay be rearranged in terms of 
parameter misalignments by letting z = e * x. Tte resulting equation 
assumes the form 

e : Ao (A - H)x 4- (b - ^)r (^7) 


Let the parameter misalignments be defined as 
(A - H) = uc T 

and 

(b - jj) - v 


( 68 ) 

(69) 


where ^ and w are column vectors and c T is a constant row vector. These 
vectors are shown in detailed form in Appendix 3* 

Substitution of Eqs. (68) and (69) into Eq. (67) yields the error 

diffential equation 


since 


s Aq 4 * u(/x) * wr 



(70) 
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then Eq* (70) may he written as 

0 = A e * ux ^ _vr ^ 70a ^ 

vter9 x and r are scalars and denote the response and stilus respectively. 
Note that the error differential equation contains the model- response 
without its derivatives and the model-reference system stimulus without its 
■derivatives. This is a very desirable result, since both of these quan- 
titles, x and r, can be measured. 

Equation (70) contains three perturbational quantities, namely e, u 
and _w. These quantities, as stated before, denote the model-reference 
system error and the parameter misalignments respectively. 

5 .J Pfimnlation of a Liapun ov Function 

An appropriate Liapunov function for Eq. (70) should be positive 
definite in the model-reference system error as well as in parameter mis- 
alignments. Therefore, one may choose again a candidate Liapunov function 

of the form 

V = e T Me + u T Nu 4. w J Qw 

where matrices K, N and Q are symmetric positive definite matrices whose 

elements my be constants, time varying and/or functions of the state 

_i_ » j. , n / rr i \ ^ + v, to tine & n cl then substi— 

variables. Differentiating Eq. (71) with respecu vo 

tuting Sq. (70) and its transpose given by 


e T = e T A T 4 . u T x 4- Jr 


(72) 


one obtains the time derivative of the Liapunov function 
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V - e T D^e 2 u T (Ku * Max 4 - Nu) 4- 2 w T (Qw 4 - Mer *- Qw) 

vhsre 

D 3 = (A T M 4 . MA 4 . M) 


(73> 

(74) 


To satisfy Liapunov’s criterion V > 0 and V ~ 0 one may solve 

Eq. (74) in tha same manner as described in the previous modeling technique 

and let 


(Nu 4- Hex * Nu) = 0 

• « 4 . 

(Qw * Kgr * Qu) - 0 


(75) 

(76) 


Equations (75) and (76) constitute the basic relationships fras which one 
realises tha adaptive controller and the identification dynamics. It 
should be noted again that the resulting V may be negative senidefinite 
or negative definite depending on the condition given bySqs. (75) ana 

(76). 


5.4 Formulation of Controller ar.d Mathematical Model 

(Nonlinear Case) 

Rearranging Eq 3 * (75) and (76) one obtains 

u = - N-^lex - r 1 Nu (77) 

i - q _1 Ker _ Q -1 Qw (78) 


Equations (77) & (78) represent the identification dynamics and the 
realization of ttese equations produces the adaptive controller, Matrices 
M, N and Q are determined in the same manner as explained in the previous 
sections# Note that r and x are scalar quantities here and represent the 
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measurable stimulus -response data* 

In order to continue with the formulation of the mathematical model 

system, one must establish a relationship among Eqs. (68), (69), (77) and 
(78). As before, one nay assume that the changes in the human operator 
during his performance are much slower than the identification time re- 
quired for the model's parameters. This restriction permits one to consider 
matrix A and vector b as being time-invariant during the identification 
interval. Therefore, differentiating Eqs. (68) & (69) with respect to time 

yields 


H = -uc J 


(79) 


and 


(80) 

£ = - W 

Substituting Eqs. (77) and (78) into Eqs. (79) and (80) respectively, 
one obtains 


H = fTHfexc 7 *■ N" 1 Nuc T 


£ = Q _1 Ker 4 . Q -x Qs 


-11 


(81) 


( 82 ) 


Integration of Eqs. (81) and (82) over the identification interval yields 


t 

H = H 0 + J (N -1 >!exc T «■ N _1 Nue T )dt 


(83) 


£ - £, 


t 

* f ( Q - Hter a Q '''Qu)dt 


(84) 


where H and jj are the initially assumed values for matrix H and vector £ 
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respectively. Matrix H Q and vector z o are given In Appendix B. The 
integrals inEqs. { 83 } and ( 84 ) constitute the identification dynamics 
and may he considered of the memory or "learning type. 

To complete the formulation of the mathematical model one substi- 
tutes Eqs. (83) and (84) into Eqs. ( 62 ). The resulting equation is then 
the mathematical model which at the end of the identification interval 
represents the actual human operator. The mathematical modal is given by 

t . 

^ = [ H 0 * J (ir 1 Max/ <■ !r 1 !iuc T )dtj X 4. 

-t 

A Zo* { (Q_1 ^ r * r (85) 

Since Sqs. (83), (84), (85) and (74) are nonlinear, time varying and 
interrelated, then one must utilise a digital computer and solve these 
equations iteratively. 

The mathematical mcdel may he considered to represent the actual 
human operator when e = 0. Digital computer simulation has shown that 
with prefer choice of N,Q and M matrices one obtains a rapid identifica- 
tion of parameters for the mathematical model. This is illustrated via 
an example later in the report* 

It is desirable for the identification dynamics to possess memory 
("learning") capabilities so that the mathematical model’s initial parame- 
ter values are permanently updated to eventually eliminate the need of the 
error information being supplied by the model-reference system configuration. 
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The interval of time required for the mcrfel to update or "learn" a set of 
parameter values which produce a model response that is equivalent to the 
response of the h«» operator being modeled could be referred to as a 
"learning" period. This "learning" period will vary for different systeias 
and will depend upon the characteristics of the stimulus signal. 

One possible means for determining when a set of nodal parameters 
have been identified is by examination of the model parameter waveform 
throughout the "learning" interval. During the initial se^- 
"learning" interval, the "learned" parameter waveforms usually fluctuate 
over a wide range of values as the model tries to conform to the reference 
system. In the later portion of the "learning" interval, the fluctuations 

, i tend toward their desired levels. There- 

die down and the model parameters tena wwa. 

t , nuc^uanons have ceased, it can be assumed that a set of 
fort, once ths flue -aa viwri 7 

E odei parameter values have been established. If a change occurs in the 
hucan operator's tracing tash, the model parameter waveforms -will again 
begin to fluctuate and a new "learning" period is initiated, un.il the 
fluctuations cease and the model establishes another set of parameter values. 
This adaptive property is inherent in the integrals- given by 3q. (83) -d 
(84) and therefore, one may consider these integrals as t*ing of the memory 
or "learning" type. Note that the mathematical model of the human operator 
contains these integrals. Hence, the model has also this adaptive "leaning" 

property • 

Th. Idor.tlfiMtlon of ,c..r , lb (85) =V b. BUM V 
„,m = ln. tb, >«™ W^tbr'. S '- <HI 


z = A z 4- br 



31 


, . (0 * . 0 and r = 1 . From Sq. ( 63 ) on,, obtains 
At t = 0, let z(0) - 0 ana r 

AX ’ (86) 
z( 0 ) - b 

(87) 


or 


z(A t)/A t = b 

des j as at tta end or the identif: l-ti« — * t.t vector , 
Since one desires au 

in tba computations 

approach vector b, then one nay set 

( 88 ) 

z(£> t)/& t = £ 

. A . ,.tnizi"e the initial conditions 

Note ttat vector £ may be detem^ _ 

* VTjatio.n (88) suggests that -s- - ■» 

a i rwPf'ren" 8 system. ->pi3 

of the modex-refer r ^ ^ ^ ^ riffl en te X 

- ^ntined durang the ^ ^ ^ tJdJ ** has been 

response data of the bunan ope curator. 

0 „ A of tape serves as — - 

recorded on magnetic tape a.« — »' 

m view of Sq. (88) one nay modify Sq. (So. as 


i = [» 0 . pw * ^ T >«]» .[* * *<** ,/ ‘ ‘l r 


(89) 


, „„ th „ simplified mathematical model of tae 

Squat ion (89) represent tne si_p 

„ , , „ t _, cki n* task. Note that tne ider.vi.-ea 

human operator when per.orm.n 0 - “ the first iteration 

ti on time has been reduced appreciably, «xnce a - 
one has to identify only the K matrix. 

5 5 T.inaar Mathe maSj^alJ^l 

’ ft ions and for specific tasks the human operator 

Under certain condxUons and f pe . In this 

3 .*» ration niay operate as a lines j 

in a closed loop control option may pe 
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case tte human operator may be described by a linear mathematical model; 
In order to formulate a linear model, one has to simplify the equations 

presented in section 5»4* 

Since a linear mathematical model is desired, then in accordance 
vith Liapunov's stability theory, M, H and Q matrices appearing in the 
Liapunov function may consist of constant elements. Hence, M = 0, 

N = 0 and Q = 0. In view of these simplifications, the time derivative 

of ths Liapunov function ba cones 


y _ _ «. 5 j T (Ku <r Hex) t- 2 w T (Qv *■ 


where 


(90) 


(91) 


= (A H * KA) 

Liapunov 1 s criterion for stability calls for V ? 0 and V 6 0. One way 

to satisfy this criterion is to impose the following conditions: M > 0, 

K > 0, Q > o, 0 and let 

(92) 


(Hu +■ Hex) = 0 
(Qw 4 . iter) = o 


(93) 


M (92) - (93) «« 

liM.r natheiiatical =«i.l ot U» "t"* 10 '- P*™"* 11 * 99*' 

(92) and (93) yield 5 

i = _ rte, 1,41 

w = - Q -1 Her 


(95) 
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(«, - («> ■» «* m ,ta ~ u “““ ° f 
llnoar S,b,tl««U»S «>• (94) — <» 1 »>” 

Eqs. (79) and (80) respectively, one obtains 


* 1 T 

H 2 N~ r£X£ 


(96) 


ter 


i = Q 

roA> and (97) ov*r identification interval gives 

Integrating 3qs. (96) ana KVU ^ 

t 

H = H 0 «■/ (H - V S xo T )dt 


(97) 


(98) 


I s -So "• 


(Q'H'ar)dt 


(99) 


»„ - fc «. ««* — - " - a *‘Z,Z7 

tMr . H„U„ (90) *»- I”) !■» < 62 > 9*" 1 


X = 


/ 1 T 

H * j (N Maxc )dt 


x + 


-1* 


£ q ± (Q~ y ^ r ) dt 


(100) 


mathematical medal of the human. 

be simplified by utilizing 


Equation (ICO) represents the linear 
o aerator. Tte identification of vector * ^y 
Sq . (88). Hence, the linear nautical nodel >» «• 


i = |^H 0 4. r (N-W/)dt| x *[i5 0 * z(A t)/^t] r 


( 101 ) 
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Tte realization of this mathematical modal is illustrated in example 2. 
which illustrates the identification of four parameters. 

In order to illustrate both identification techniques two "artifi- 
cial" examples are considered. The reason for presenting these "artificial" 
examples is to show the usefulness of the identification methods and the 
information ore must know a priori about the reference system so that the 
resulting model is a true representation of the human operator when per- 
forming a particular task. By an "artificial” example is meant that the 
reference system doss not represent a true human operator performing a 
prescribed task, but instead the reference system repressnts a dynamic 
differential equation with known coefficients. This differential equation 
is then subjected to a known stimulus and. the generated response with the 
stimulus are recorded on tape. Only the stimulus and response data are 

utilized in the identification process. 

The purpose of illustrating the identification techniques via these 

examples is twofold. First, it gives the designer a means of verifying 
tbs validity of the controller or control laws. Second, one has a check 
on ths identified mathematical model parameters. In addition, one obtains 
an insight and experience in choosing the M, H and Q matrices in order to 
shorten or lenghten the identification time. Once this knowledge is gained, 
one may undertake the modeling of the actual human operator data provided 

it is recorded on the tape in the required fora. 

The first example illustrates ths identification of two parameters 
for a second order mathematical model expressed in the controllable form. 
The second example shows the identification of four parameters for a second 
order mathematical model expressed in the observable form. 
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6. EXAMPLES 

6.1 

. ivi,,qtrat0S the identification of 
Consider Fig. 1. This example ill^tra 

. , p . ond war mathematical model. The desired re- 
two parameters for a second c.der ma 

, moiol i S a state determined controllable 
alisation of the mathematical mod.l 

• nf the state variables is in accordance 
form. Therefore, the choree of the state 

with the desired realisation, ^ ri3r _ 

- n f tv reference system, denoted by *. » 8“ 

Let the response 01 t..e re 

ated by the differential equation 


S * 22s V 121* = 102f * 1£7r 


(n) 


. „ i The initial stimulus is 

uWe r „ _ y and * * y *s »ho« l * TbS 

- tK . plant. The information recordea on 
applied to the last stage o. t.„ pL9 

. tv initial stimulus , the reierence syste- 

magnetic tape consists o. • - e^baok 

*-.» that the plant’s response is fe-«w 
response and plant’s response. - ' 3 the re f e renoe 

to the input of the reference system. If one -.re P ^ 

system with an actual hunan operator whose task would be o con 

th „ ,* «.* “““ “ “• 
tor hov to move the stic* 
plant to conform to a desired perforoa 

P ^ a , chosen for the reference sys-ecu 

A tentative mathematical model is chosen 

The model has the form 


Z 4. V * V = 102f * lS7r 


(32) 


j k it this ooint one assumes 
at tand is to identify b 2 and bj_. 

The problem at ha assign initial 

» 2 and d, are e»,l.W -*»». "»»■* lh “ 
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values to b 2 and b^ The restriction on the initial values is that they 
must lie within the stability region of the tentative mathematical model 
as determined by a Liapunov function generated for this model system. 

The next step in the modeling process is to realise a controller 
or control lavs which update the initially chosen values of b 2 and bj_ so 
that the final mathematical model represents the reference system. In 
order to accomplish this task, one makes use of the derived equations. 

Using a phase variable representation of Sq. (32), the B-matrix has the 

fora 

0 1 
- b l ** b 2 

Therefore, in view cf (20) 

i 

j 0 0 

3 =| 

_ bi - b 2 

where 

0 

1 

and 

"1 _ “h (E6) 

"2 -h 

Note that the dot denotes the time derivative and superscript T denotes 
the transpose* 

Since a linear model system is considered, then using hqs . (32) and 









V 


(S6) yields, after some 
tion dynamics 


. the following ider.tifica- 

alrebwic Efini?ul a -- rns ' 


\ = U/n u ' 


■ z ‘ ( a 12 e * °22^ 


Integra 


b 2 = (1 /r* 22 ) ' ' (c 12 e V 0 Z2 9) 

tion of 5qs. (37) and (38) g ives 

£ 

■ Ol =\0>J U/^ • * ■ 12 S * E22< 


e)dt 


(S7) 


(S8) 


(39) 


b 2 51 b 20 


0 * / U/»22> • 1 • ( "l2 e * E2: 


,e)dt 


(510) 


. b .. the initially ^ ^ 

“ ° 10 20 • . ard 4 are entities. Note that 

tively. The variables z, a, • _ are 

or . the taco. "lerer.vS _ 12> “22’ il 22 

z and z are record • are orap uted 

t . Lian-onov function. The elements - 12 • 22 

obtainea fron - - 2bd n are free design carters 

{*, o). Sler-ents n-n ana ri 22 

frca 5q. (31) °- — 7 0 . Using 5q3 • (31) and (53) 

and =wi t satisfy the conditions n u ? 0, n 22 

with a diagonal D-rcStrix of t..e 

0 I (3U) 


D = 


- d 


11 


0 


- d 


22 


and a syn^ne 


metric Matrix of the font 
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K = 


*11 *12 
*12 ra 22 


(E12) 


0 re obtains 


ci 2 = d u /2b i 




22 


= (d 22 * d 11 /b 1 )/(l/2b 2 ) 


anc 


*11 


(b 2 d ll * b l d 22 * t 1 d il )/(2b l b 2 ) 


(S13) 

(EH) 

(S15) 


Without less of insight into the nodsling procedure one nay let d n = 1 


end d 22 = 1. Equations (S13 to 315) reduce 


P-12 = l/2bi 


to 


ana 


= (1 * l/b^/U/^) 


= (b 2 <- b* u b 1 )/(2b 1 b 2 ) 


(E16) 


(E17) 


(S18) 


. . , 1 11 ^ *v«, =■) =-=-ts e, , , a. - and E-- are dependant 

Note that if one uses 00 . UH t..e eie— ..i» -j_i> X 2 22 

upon the b *s. Therefore, during the identification process, one must 

ccccute r.ev sets of a^'s at each iteration step. Since the m^'s appear 

in the Liapunov function, then one nust also compute V and V at each 

iteration step using 3qs. (8) and (30) and check if Liapunov's critericn 

is satisfied at each iteration. 

Substitution of Sqs. (39) and (E10) into 3q. (32) yields the uatfce- 
Eatical eodel for the tape recorded data. The model is described by 
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* *l> • l 




4* * 


’£ b 10 *■ j O/n^)- z-(°i2 e * a 22 s ^ dt ] x " 102 ^ + 187r 


(319) 


vh9 „ the Bij 's «. given by .,s. (316 - 318) and «. -*«- *" ^ ** 
diagonal N-natrix used in the fetation of the Liapunov function. 

As mentioned before, or* «.y also use 3q. (50) to find a set of 
B . 8 needed for the controller identification dynamics. Squation (.0) 
contains the reference system matrix A which is ^nown. Ho information 
is recorded on the tape. In this case, th, generation of the n's 

, n . \ c »q tvs chOS C I"i S C C OT d_L H e, wO "th© 

a trial and error approa ch* ibs j 

_ , T r _ l£) , T'nis theorem states that if * 

the ore ni given by LaSalle and Lefscns-2, * -* 1 

, *v*„ the re exists a positive definite symmetric 

represents a stable system, then tnere ex^sr ✓ 

+ ■ ^ sv^-stric Do matrix* (See -»Q.« (50))« 

JUiratrix for ovary negative defini- 2 

U7 \ /-/l\ and fA-17) one obtains after some algebraic 

Using Bqs* (t3), (obj ana (a * j 

manipulations the identification dynamics given by 


where 


by = (SB) ♦ ( n 22 x “ n 12*^ 
b 2 = (SB)-(n u x - r^ z x) 

(SB) = (e^jS * m 22 e)//(n ll' n 22 ” n 12 ) 


(320) 

(321) 


Integration of Sqs. (320) and (321) gives 


\ ~ ho v / ( S3) ' (n 22* " n 12 x)dt 


(322) 
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x 

= b 20 */ ( S 5 '-< n U* " D 12 * ,dt 


(E23). 


vhsre a 


that a norwi 


and / s are 


nc b„ are the initial values of ^ and b 2 respectively. Note 
diagonal S-catrix was used in tbs above equations. The B y ’s 


obtained Trcs. the Liapunov function. 


Substitution of Ecs. (S22) and (S23) into Sq. (E2) yields the 
nathenatical node-1 for the tape recorded data. The nodal is described by 


•a, — 

x *|b 23 ♦ J( S2)-(r. u i - n 12*' d j * * 


= 102r * l£7r 


(524) 


b 10 4- J !S3)-(n 22 x - n 12 xUt 
0 

. q -»_._* ■* r r- o ^ ^ ta X c 0*^ nut 

Folic vine is a list of nvn^ricsu values -s-u ... t — -© 


program: 


V 5 * r* ^ V ^ 1 Q T 3 nt S 

B n = 100 

K-r>a t.riy sler.er.ts 


k ll 


12 


- 27 


= .9 


Doa = 27 


“22 


n 22 = 1 


Q-rnatri'* elements 


*11 


= 1 


q 12 = * 95 


mm .1 Parameters of_MoM 


b 20 = 50 


b 10 = 50 


n „ = 1 
*22 


I te ra tlon Interval 


At - 0.0X1 sec. 

The digital ccrr.puter results are show, in section 12. 
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6.2 2. THpnt.i fl cation of Four gaj3Bgi*M 

Consider Fig. 2. This example illustrates the identification of 
four parameters for a second order mathematical model. The desired re 
alization of the model is a state detained observable form. Therefore, 
one chooses the states for the model in accordance with the desired 

realization* 

~ + u p rpfprence system, denoted by a, be generated 
Let the response of the re-ter-nce j > 

by a differential equation 


d 2 z/dt 2 * 22dz/dt * 121 z = 102dr/dt * l 87 r 


(3255 


where r = - y and z = d 2 y/dt 2 as shown in Fig. 2. The initial stimulus 

„ ,P P u.d «. .« ,u»< «• pi*"''. "*•"“"** « «*** “ *“ •'* 

of the reference system. Fig, 2. may be considered as a representation of 
a tracking task, in which a human operator observes on a visual display 
tte slant's response y and adjusts a manipulator, joystick, handwheel or 
similar output device in such manner that the plant's response conforms to 

_ Thp information recorded on magnetic tape consist 

a desired performance, *he 

of the initial stimulus, the reference system's output and the plant's 
response . Note that the reference system's output, denoted by z, represents 
the' adjustment which the human operator imparts on the plant by manipulating 

the stick* 

In order to describe the recorded data via a mathematical relation- 
shlp , one may choose first a tentative mathematical model. Let the tentative 

mathematical model for the recorded data be 


d 2 x/dt 2 4- h^dx/dt 4- h Q X = g^dr/dt * g Q r 


(226) 
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where 1^. h Q , Bl g 0 are to be identified frco the recorded tape data. 

Koto that the problem at tend is the identification of four paranetere 
based on experimental recorded data. Rearranging 3q. (326) yields 


d 2 x/dt 2 v dlh^x - g x r)/dt h Q x = gj 


(E27) 


Defining the state variables as 


X 1 = x 


X2 - x-i * (^l x l - g]_ ] 


(E28) 


the state equa 


te eouations for the tentative mathematical aodel are 


*! * " Vl " X 2 * S i r 


Vl * g o r 


,, x j x f rc -« this choice of state variaoleo is 

The H-astrix resulu^rj*, irc^ b.u.- 


- h 1 


_ Vi 0 

- **o 


md the _g- vector is given by 


(329) 


(330) 


(331) 


note that h x , h of gl and g 0 are unknown quantities. One may assume initial 
values for h,, h D , g, and g Q with the restriction that the assumed initial 
values lie within the stability region of the tentative mathematical model. 
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This stability region may be determined by a 


Liapunov function formulated 


for the model system* 

The next step in the nodeling process is to find the control law 
which update the initially chosen values of hj_, h Q , g 1 and g Q , and perait 
the realization of the adaptive identification controller. 

Using the equations derived in section 5. one obtains 


where 






(532) 


(333) 


and 


0 

1 


(234) 


Note that the dot denotes 


the time derivative 


and superscript T denotes 


the transpose* 

Since a linear mathematical model is being considered, then using 
Sqs. (96) and (97) one obtains, after seme algebraic manipulations, the 
identification dynamics for the model parameters. 


hj_ = (x/n 11 ) ■ (n-jje «- a 12 ®) 
h Q = (x/n 22 ) • (n 12 e * nrf) 

§1 = ' (n ll s V *12® } 


(335) 


(236) 


(337) 
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, •, (E38) 

g 0 = (r/q 22 ) • ( m 12 e * m 22 e) 

Note that matrices N and Q are assumed to be diagonal for simplicity. 
Integration of Eqs. (E35 - E38) yields the nodel paraoeters 

(S39) 


h l = h 10 *■ ( ( ^ n ll' ' V n 12 s)Qt 


t 

= h oo * J (x/ n 22> ' (m 12® * B 22 e)dt 


(E40) 


V 

gl = h0 *f (r/o u ) ■ (m n e , nrf)* 


(su) 


" & 00 


z 

- j (r/q 22 ) * ( K i2 e * “22 


g'ldt 


(E42) 


x* + pec,— .x>d values for the so^el s 

v v cr q r,d " are tne initial assu-^i 

vnere h 1Q , ^ 00 > g 10 a ‘ °oo 

4 vo e r* Y G and e are measured quantities and re- 
parations. The variables r, x, e ana 

j_ + Mn j a i r^> f'er^nce error and its time 
present the input, nodal's output, nodal 

derivative respectively. Elements » 12 , “ 22 ’ "ll' °22’ q H q 22 

are obtained from the Liapunov function. Elements n u and q u are 

design parameters and must satisfy the condition n u 7 0 and q u > 0. 

- +V . Q vt , nd o ratrices which are chosen for this 
These elements arise from the - and 

, + 4 . Thp 's are chosen in accordance with 

example to be diagonal matrices. The m ij are 

(18) 

the previous stated theorem 


«> (^39 - E42) into To. (S2&) yields the 

Substitution o* ^qs • V-iV 
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mathematical nodel 


xJ h 10 * |(Ai)'(vv ,dt ] i ' 

l n 




(x/n 22 )' (n 12 e 4. n 22 ' 


sHtjx = 


• f !l5 ./h/V".rv4* 


* s 


00 


* j (r/q 22 )>(n 12 s * ffi 22® )dt ] r 


(343) 


, . . -1 nj-i -if ons assuces that tr.e 

Equation (343) =a 7 *>e si 2 ?xifiea 11 on 

systsa and the aatheaatioal nodal are at rest at t = 0. 

* (0) = 0 and Xj(0) = 0 at t = C. Then 3q. (329) S^' 3 


\ = g l 

*2 “ S o 


rafersnos 


= 1 1 


(344) 


(345) 


x^(A t)/A t = g]_ 
x 2 (&t)/At = g 0 


(346) 

(347) 


Since one desires 
interval, then, using the 


b = g Q at the end of the identification 
orded tape data, the b's r_ay be ccaputed at 
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the first iteration by 


!,(ot)/At = bj = g. 


z 2 (At)/A t = b 0 = g 0 


(S43) 

(346) 


Note that vector s can be determined from the initial conditions. Thus, 
at the first iteration interval one determines vecvor ,g« The renain_n Q 
identification time is then utilized to identify the elements of the 
H-raatrix. This simplification reduces the time for identifying tne mcdex s 

parameters . 

The simplified mat heretical model has now the form giren by 


t 

t’ J h 10 4. f (*/ »!,)*(»!!• 


[lice * / (x /n 22 ] ' {r ~12 e * a ~J )i: 




>10 


+ a n (A t ) / At r 4- 


4 . z,(At)/At r (250) 


where z. and I, are defined in the sane manner as x x and x 2 . 


Following is a 
computer pro 
M-matrix elements 


list of nune ri ca 1 values used in toe digital 


“u = 730 


n, 0 = 100 

JL 


m 22 = 33 
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H-ratr i x elements 


n n = °- 5 



22 


0.05 


Q-matrix elements 

q X l = 1 °*l 2 = ° q 22 = 1 

Initial parang tor s of Model 

h 10 ~ 2 h oo = 200 g 10 “ 0 g oo ~ 0 


Iteration interval 


t = 0,0001 sec* 

The oror-Traii; and computer results are shown in section 12* 



7. CONCLISIO'vS A'O RZGCSOS NDATI0K3 


7*1 Conclusions 

Tvo adaptive modeling techniques to identify and control a human 
operator's riatbematical model has been presented* The modeling tech- 
niques offer a theoretically consistent approach for identifying a mo^sl s 
parameters free experimental data, A very useful feature of this study 
is the development of a digital computer program which is easily imple- 
mented and modified concurrent with experimentation. In this way, the 
modeling process interacts -with the experimentation process in a mutually 
beneficial vay* This is shown via the logic flow diagrams in rigs. 3, 

4, 5 and 6, !<ote also that it is tie systematic and logical use of the 
digital computer that permits one to effectively apply the class of 
Lisouncv functions no the modeling of a human operator when performing 
a tracking task. 

One of the objectives of both modeling techniques has been the 
synthesis of a controller or control lavs that provide the identification 
dynamics for the model's parameters. An important feature of these 
ident if icat ion dynamics is their "learning" capability* This memory 
or "learning* capability provides a means for determining when a set of 
model parameters has been identified. All one has to do is to examine 
the model parameters throughout the "learning" interval, As it can be 
observed from the digital computer printouts, during the initial segment 
of a "learning" interval, the model parameters usually fluctuate over 
a vide range of values. In the later portion of the "learning" interval, 
as the parameters approach the correct values , the fluctuations in the 
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parameters diminish and the model parameters tend toward the correct 
values. Therefore, once the fluctuations have ceased, one may assume 
that a set of model parameters has been identified for the human operator 


performing a given tracking task* 

The "memory" or "learning" capability of the identification dynam- 
ics has another important feature. It permits one to establish a cor- 


respondence between reference system’s disorders and parameter value 
changes in the mathematical model. This may be observed from the follow- 
ing experiment. If a change occurs in the human operator's tracking 
task, the model parameters begin to fluctuate again and a new "learning" 
period is initiated, until the fluctuations cease and a new set of pa- 
rameter values has been identified for the mathematical model which rep- 


resents the human operator for a changed condition. 


This adaptive property 


is one of the main advantages of both modeling techniques. 

Perusal of the modeling techniques indicate that the identification 
dynamics developed for the mathematical model in the controllable lorm 
require the measurement of the entire response vector as well as the 
stimulus vector. This is due to the fact that the mathematical model has 
been synthesized based upon a phase variable representation of nq. (2). 

If some of the phase variables of the stimulus and response are not mea- 
surable, then the identification laws are inadequate to update the model's 
parameters. To alleviate this situation, one may generate the additional 
required phase variables by successively differentiating the last measured 
phase variable. The resulting phase variables are then noisy and may con- 
tain discontinuities. This produces oscillations in the model parameters 

to determine the identified set of model parametes. 


and it beccmes difficult 
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\ 

Since these oscillations are steady and oscillate about sense nominal . 
values f then one may consider these nominal values as the set of iden- 
tified parameters for the mathematical model. 

To overcome this disadvantage, the second modeling technique 

which realizes a mathematical model in the observable form has been 

» ro,l.v T **“““* "" “ d ~ “* 

tb. •»««. *" a "'“ r “ i,r 

«, identification «—-• » “ 

the scalar re.ncn.. and ati.nl™ " « « ““**» b? E< >" ‘ 8 ” 

and («). in addition t. thla «- ““““ ay ”* 1 “ 

of tll , seed nodding »«pl» *“ 

Tharafore, no. nay aon.ldar tb. id.ntlfio.tidn t.obni,™ lading to tb. 

nathacatioal -odd of tb. ob.ary.bl, torn „ find — «» “> 11 “- 

tifi cation technique which prcduces a mathematical model of the control- 
lable form. 

The identification demies have been derived from a class of 
Lia duhov functions which possess variable characteristics. An important 
feature about Liapunov’s Direct Kethod is that Liapunov functions are 

arbitrarily ebeen for the identification dynamics. The condition these 

, t , is that H > 0, H > 0 and Q 7 0. Experimental 
matrices must satisfy is inax n * 

. m thflt matrices K, H and Q affect the rate of model 
observations have shown that matrices , 

parameter convergence to their correct values. 

A review of the identification dynamics equations show that the 

factors on the model- 

elements of K, N and Q matrices serve as weighting 
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reference system error as well as on the response and stimulus variables- 
Several sets of values for the elements of M, N and Q matrices have been 
considered. Two of the better sets are shown in the computer printouts 

in section 12. 

The digital computer results given in section 12 are self- 

explanatory. I 

| 

i 

7,2 Recommendations 

It should be emphasized that the adaptive controller is intended 
mainly to aid in the identification of the mathematical model's parame- 
ters. In general, analysis of the results obtained by using the adaptive 
controller does not necessarily dictate that the actual human operator 
under study must be modeled in the proposed manner. Instead, analysis of 
the results suggest certain modifications which should be considered in 
future studies. One of this considerations is a better and clearer 
delineation of the experimental data. Note that the recorded experimen- 
tal data needed to formulate a mathematical model of the controllable 
form is different from the recorded experimental data needed to generate 
a mathematical model of the observable form. Therefore, one of the recom- 
mendations of this study is to first decide on the desired model structure 
and then perfora the experimentation with the human operator and record 
the data needed for identifying the model's parameters. Should the experi 
mental data be difficult to obtain, then one can modify easily the model 

structure to accommodate the experiment. 

The following areas of research are logical extensions of this 


study: 
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1. Develop a design criterion that would permit the mathematical 
model parameters to become identified in the shortest time* 

2m Investigate new forms of Liapunov functions from which dif- 
ferent identification dynamics may be obtained, 

3, Study the usefulness of the presented modeling techniques to 

the development of diagnostic techniques. This recommendation may b3 

| i 

useful for investigating abnormalities in a human ope rat or T s performance 
with a particular plant. It may also be useful for investigating human 
operator performance with different plants. 
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9. APPENDIX A 

This appendix gives the detailed forms of matrices and vectors 
used to derive the mathematical model in the controllable for*. 
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( A— 18 ) 


(A-20) 


(A- 21) 


(A-22) 
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Additional „ InforpaVioj} 

. + vr> ^rpr^n-'e system's scalar response 

z = z(t) represents the reference 

(measurable). It is the stick output. 

x = x(t) represents the mathematical D odel system's scalar response 
(measurable ) . 

. = «• QO.aUraCor.nce ooopo.i.. •!»««'■ ■«>" 

error (measurable) * 

y = y(t) represents the plant's scalar response Measurable) 

_u = u(t) and _w = v(t) -present the parameter misalir-nt vectors, 
r » r(t) represents the scalar stimulus. 


M = dM/dt 
H =■ dN/dt 
Q = dVdt 
u = du/dt 
v = dv/dt 
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XO. APPENDIX B 

This appendix gives the detailed forr* of ^trices and vectors 
^ed to derive the statical .odel in the observable for.. 
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where x and x are given by Sqs. (60) and (61) respectively. Vectors 
z , i , e a 


nd e are obtained in the same mnner 
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h (n-l)o 1 ° * ’ 

^-2)0 0 1 * ’ 


0 

0 




h 10 

^oo 


0 0 

0 o 


« • ♦ 


• • • 


0 

0 
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g (n-l)o 

g (n-2)o 


(3-20) 
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*00 


<*11 d 12 *13 * • * *!» 

d a d 22 d 23 • * * d 2 “ 


d nl d n2 ^3 * ' ‘ 


(B-a) 
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D/. = 


*11 


12 


d 12 

d 22 


13 


23 


In 


2n 


(B-22) 


In 


2n 


3n 


nn 


or 


D , = 


l ll 


22 


. 0 

• 0 
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0 


0 


. d 


nn 


(B-23) 


scalar response (measurable-. 

Z = z(0 represents the reference system sc- 

It is the stick output* 

x B x( t) represents the mathematical node! sysWs scalar response 

(BeaSUrable !* ode , r , f ^, oe coneosite system's scalar 
e = e(t) represents the model-relerenc 

grror l measure ole ) » 

_ . J(t) represent. tk, stent- ™J“* (»»«»“•)• 

„ . sW .pd j *«> represent the .->« 
r = rlt) represents the scalar stimulus. 

M s dM/dt 
N = dN/dt 
Q r dO/dt 
u - du/dt 
y = dv^dt 
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FIGURES 
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^ _ rt a, ti * on V* t Xsth6u3ti csl 

FIG, 2* Model - Hsferencs oysters uoim-gu-^.on 

Kndcl's Response Applied to Controller and e = x - z 

Uote: To realise the as t hectical ro.isl in the observable 


2 - X 


























FIG. 4. Logic Flow Diagram - Linear Model in Controllable Form 

















Update 
for next 
iteration 

















Head: r, N, Q, jc 


Ccomute x using Eq. ( 62 ) 


Compute e using Eq* ( 64 ) 


Choose > 0 and occiputs M 

using Ea. (91)* (3eo Ref, 1 


Compute K~ x , O' 



Compute u, v, u, w 
using Eqs* (94) & (95) 


Update H and compute £ and updated 
x using Eas, (98), (88) &(101) 


If e = 0, u = 0, w = 0 

see Eqs, (68), (69) & (70) 


Compute V and V 
using Eqs, (71) & (90) 


. If 

V £ 0 

V > 0 



Update 
for next 
iteration 


Select new B Q , £ 0 


Logic Flow Diagram - Linear Model in Observable Form 












12 * DIGITAL COMPUTER RESULTS 



F0R TRAN- IV- C-LEVEL-^ZO MAI-N OAT.E~*~723 13 

0001 REAL IN1.IN2 _ 

CC02 t DATA Z1.Z2.Y1.Y2/0. *0* .0* *0. / 

. 0003 ’ ..DATA R1.R2/O..0./ — 

npr>4 . • REWIND 9 

0305 READ ( 5. 1 CCO) Al » A2.* HL..» H2. 

or 0 A WR ITE< t, 1000)A1 ,A2 »H 1 *H2 

00 0 7 . 1000 FORMAT! AG 10. 3) 

00 OR D T= l • £- A 

DOOR c ' C ” INPUT SIGNAL — * 

... 03 ID 1 INU1.G 

001 1 ' 00 IOC J«1 *2 COO 00 

00 13 . P 1= IN 1- V 1 

00 1 A R 2»(R l-R 10 J /DT . „ 

c : REFERENCE SYSTEM OUTPUT 

... 001 5 _Z 2D»~Aini-A2*Z2+Hl*EU«-H2*^ 

0016 . Z2«Z2*0T*Z20 

00 17 J Z 1 C Z 1 *0 T *Z 2 

• • . c OUTPUT LIMITER AND DOUBLE INTEGRATOR 

nni« IFIAOSIZ D.LT.1..0R. ABSIZIJ .GT..ODGO TO 2500 

.3319 1 FI ARS< ZLI .3 T. 1. I 2 L = S I GN( l . ,Z U 

_ 0320 ; : _1F(AUS(Z1).LT...C1) ZL-0.0 

CO? 1 . • GO TO 2501 

00? 2 „ 25CO 7L“Z1 .. ... ... .. 

03 Pi 250 Y?»Y2*IC*ZL)*DT 

0!3?A , . Y l-Yl*CC*V2)<‘DT .. 

m?5 ••. . 2 F GRKATI I 6. 6E 1 5. 7) 

. 00?6— — WR ITE( 9 . 2) J.R1.R2.Z1 ..Z2.Yl.Ot2 

00 ? 7 . • IF (MC)Ot J . 20C0) .EO. C) wRI TE 16.2) J.R1 .R2 ,Z1 .Z2.Y1, Y2 

00 2 8 ICO CONTINUE „ ._ . ..... 

OO ? 9 FNO FILE 9 

0330 REWIND 9 . _ 

0031 STOP 




FORTRAN 

000 1 
000? 
0003 
0 0 0 '« 
0<;05 
0008 
00^7 

000R 
r.or.Q 
00 1 0 
coil 

CO 1 ? 

00 \ 3 
00 1 ** 


C 015 
CO 1 6 
001 7 
0.0 1 3 
001 9 
0 0 ? 0 
017 1 
no?? 


00? 3 
00 ?* 
00? 5 
00 ? 6 
00 7 7 

C.0 7 8 

017*5 


0030 
00 3 1 


003? 

0033 

003* 


00 3 3 
00 3 0 
013 7 
00 3 3 

0030 

nv.f) 

< 1- 1 


iv G LEVEL . 20 


MAI N «. 


OAT E 


77.313 


16/52/32 


I MPLICI T REAL (M,N) 

INTEGER mud 

DIMENSION JO T l 3 201 
OA TA XI . X? . 2 1 * l 2 /C » \ 

_ DATA ft 1.R2/0. • C. / 

data Y1.Y2/0* *0. / 

REWIND 15 

INPUT CP 

READ! 5 • 1 0 CO ) B 1 *U 2 *C 1 *C 2 

R FAOl 5* ICC U Mil . M2 2 ♦ Ml 2 

.r [A(H 5 , 1000 > N U ♦ N2 2 * CL1 »Q2Z ._ ... 

’ RE Al> l 5. 10C2)0T 

• READ! S? tl JOT 

SlJCTf-tN CARDS (OR COMMENT MUST BE 
U l lr. SORT ( . 0*Nll*N?2 > 

Q 1 2 = SORT ( . S«uU*0?.2> 


i 0* * o* / 


REPERENCEt MODEL AND 


MAT RU PARAMET ERS — - 


PLACED AT DATA'S END. PLUS 2 BLANKS 


CHECK 

W R 1 T F I 6 * 7-0 )i) UB2.C 1 »C2 

kR1Utt.23»OT 

wR I T E t 6 » 2 5 ) Mil * M2 2 t Ml 2 

WR l TC- ( 6 . ? 1 ) Nil • N?2 *N12 *011 *0^2 * 12 

v,R ITE1 t.2C» JOT 

' WR 1 r E ( 6 ♦ 7 . 2 ) 

T = C . 

DO ICO J»l. 100000 


OP COMPUTER READING 


20 

22 


PR IN T OUT 

1f< .NOT.I J.eo.l.CR.«CO<J.|«).E0.3IICO TO so ... 

1. ' X2 ' , TBS. 'I 2' « T ^ 7 . ' VI ' . TL 09. Y2 / I ^ 


.SO CCJNT1 NUE 

t«t+dt 

READING the data tape 

READ! 15.2>J«R1*R ?* 11 » Z2 ’ Yl ,Y2 
2 F URMA T(I 6. 6E 1 5* 7) • 

model system state EQUATIONS 

X2D=-BUX1'B2*X2*C l*Rl+C2*R2 
X?*/2 + 0T*;<2D 

X1=X1 + DT*X2 .... ... .. -• 


adjustment calculations 


E 1 

r ? 


x 1-2 l 
X 2- 2 2 


Mi»IKl?*U+^22*E7) /( NIl*N2Z-NU**2) 
\ jDasn-X N??o XI-N12 *X21 

U ?!)■* SO * t NU*X2-N12^Xl> 
h i .iitiu^.n 
n .» w ?i n i ■''' .’ i 


ORT fi AM— IV-C—L 6VEL 


PAGE -0002- 


20 FORMAT! ' 1 LYAPUNOV STABILITY MODELING Oh SEC CM) CKOER SYS 
l • STARTING PR A METERS CP THE MODEL SY $T £,M. . . ■ // I 10 . 1 U 1“ '.1 
?30.*B2“ '.El C. 3. T60. 'Cl » • ,E10 . 3 .Tti3 . ‘C2* 'ffcl0.3///l 

Jl FORMAT!' MATRIX ELEMENTS . '// * NIL 8 1 ,lPE10.J t 10X r *N22 

_ 13, 10X. 'N12** * * E 1 C. 3 / ' Oil* 3 * ,E10.3 .LOX . '022* • . E 10. 3. 10X, ' 

210 - 3 Jill \ 

23 FORMAT!* DELTA T= '.1PEL0.3//J . 

25 FORMAT!' M MATRIX ' // ' ML 1 B ♦ , l PE 13 . 3 » 10X » M2 2» ' , E10 

1M 12* ' .E 1C. 3//) 

26 FORMAT! 16! 20AA/) //I 

.1000 FORMAT! AGIO. 3) — 

1CC1 FORMAT! 3G 10. 3) 

.ICC2 FORMA T! G10.3) - - - —• - • - 

STOP 


TEM' //// 
PE 10. 3, T 

> * , E10 . 
Q12» ',£ 


.3. 10X, » 


lYAPUNCV STABILITY MODELING OF SECOND ORDER SYSTEM 

ST ART 1 AG FRAME TER S OF THE MODEL SYSTEM... 

31* 5.CCCE 01 B 5* 5.Q0CE Cl _ ... C1-. 1.0706 02 


DELTA T* l.OOOE-OA 
M MATRIX...... 


Ml 1 = 

1 .030 F n? 

M2 2- 

2. 7 CCE 

01 

M12 » 

2. 700 E 01 

MATRIX 

ELEMENTS. .. 

• m • 

.... - 

- 

- 

... _ 

M \ r. 

1.C0CE CO . 

N 22* 

1. OCCG 

oc 

N12- 

3.9V.E-0 1 

Cl l* 

1 .00 0 E 00 

022= 

1. CCCE 

oc 

0 12 = 

9. A87E-01 










C 2- 1.020E 02 



TS5~ 

O 



1 ' 

0.0 

5.C00E 

Cl 

200 

1.990 6-02 

5.8536 

01 

800 

3. 9901-02 

5.7365 

01 

603 

- 5.990F-Q2 

— 5.65CE 

Cl 

aco 

7.9H9F-02 

6.7726 

C 1 

1000 

9. 9886-C2 

7.3336 

0 1 

1200 

1. 199F-CI 

1. 269 E 

02 • 

- 1800 

1. 399E-01 

6.7766 

01 

1600 

1 . 599 F-C 1 

5. 9866 

0 1 

1800 

1.798E-C1 

— 6.8836 

CL ... 

2000 

1.553 E-0 1 

6.8816 

Cl 

— 2200 

2 . 193E-C 1 

6. 9 2 8 E 

01 

2800 

2. 35EE-C1 

7.0676 

01 

2600 

2. 598E-C 1 

7. 3 C 8 E 

01 

2R00 

? . 75 B 6-0 1 

7.88 6E 

Cl 

— . 3000 . 

2.558E-0 1 

7.878E 

Cl __ 

3200 

3. 158 F-0 1 

7.685E 

01 

— 38 00 

3.398 E-OL 

8 . 18 2E 

01 

3600 

3. 59 86-Cl 

8.0 6 8E 

01 

3800 

3.798 E-0 1 

7.8886 

01 

8000 

3.598F-0 1 

7.630E 

01 

— 8 200 _ 

8. 1 57F-01 

. 7. 1176 

01 

<.800 

8. 3576-C1 

6.6856 

01 

- 8630 . . 

8 .557F-C1 

, 6.2185 

Cl .. 

8800 

8 . 757f-0 1 

6. PC EE 

Cl 

— 5000 ... 

8.597F-01 

. 6. 7 5 EE 

01 

5700 

5. 197 F-Ol 

9 .8 20b 

0 1 

5800 — 

5. 35 7E-C1 

... 1. 3 9 9 E 

02 _- 

5600 

5. 597F-0 1 

8 . 1 H 5E 

01 

5 HOG 

5. 757E-C1 

i. me 

02 .. 

6000 

5. 5976-C 1 

8. 785E 

Cl 

.. 6200 . 

6. 15 7 f-01 

. 1 . 2 3 3E 

C 2 .... 

6800 

6. 35 7E-C1 

5 . 7 8 8 E 

01 

6600 

6.596F-01 S. 8536 

01,_ 

6 800 

6. 7566-01 

1 . 182E 

02 

7000 

6 .996E-0 1 

. 1.085E 

C2 . 

7 700 

7. 156E-C1 

5.8 13E 

01 

. 7800 

7. 356 f-C 1 

1.072E 

02 . 

7 60 0 

7. 5 9 6 E~f 1 

1 .05 IE 

02 

— 7800... 

7. 796C -01 

_ 1.CC3E 

C2 __ 

H 000 

7.996F-C 1 

1 . 0556 

02 

8 700 

8. 1 96 F-C 1 

1.0056 

C 2 

6 8 00 

8 . 386F-C 1 

1.0365 

C 2 

8 6 00 

8. 596 F-C 1 

1 .00 25 

0 2 

8 800 

8. 795F-0 1 

1 . ( 0 36 

C 2 

9000 

8 • 9 9 *> t - C 1 

9.9 956 

0 1 

9 203 

9 . 1 9 « f - n 

5. 7l-9f- 

Cl 

9 8 00 

9. 395F-01 

5.5 UOf: 

01 

9600 

5. 595 F-f 1 

9.8865 

0 l 

9 800 

9. 795E-C 1 

9. 5 3 9 fc 

01 

10000 

5.9S5E-0 1 

1. C10C 

c? 

10 200 

1.0196 CC 

1. 1 135 

02 

10 8C0 

l ♦ C 3 9 f CO 

1. it 26 

02 

10600 

1.059E 00 

1. 1 3 5fc 

02 

10800 

1.079F 00 

) .030E 

02 

UOOQ 

1 . C ‘J f 00 

1*1036 

02 


B2 XI 21 


5. OOOE 

01 

0.0 


0^0 


8. 5086 

Cl 

1.688 6 

GO 

1.668fc 

30 

8.2 CCE 

01 

2. 731 E 

00 

2 .7 33 E 

00 

. 8.327E 

01 

.... 3. 3536 

00.. 

_ 3 . 350 t 

00 

2. 9b 86 

Cl 

3. 6356 

00 

3 .680 6 

00 

2. 728E 

Cl 

3. 71 0 E 

00 

. 3 .63 2 E 

00 

-2.629E 

01 

3. 628 E 

00 

3. 5806 

JO 

3.725: 

Cl 

3. 31 95 

00 

3.2536 

33 

8. 6 52E 

Cl 

2. 938 E 

00 

2 .8 78 E 

JO 

3.7896 

01 

.. 2. 887E 

00 

2.822 6 

00 

3. 6156 

Cl 

1. 981 6 

00 

1.9 16 6 

00 

3. 6856 

Cl 

1.8826 

00 

1 .378 fc 

00 

3. 5356 

Cl 

8. 3896 

-01 

8 .250 6 

-01 

3.276E 

01 

3. 3836- 

-01 

2.707 fc 

-01 

' 3. 128E 

01 

-2.0986- 

-01 

-2.727E 

-0 1 

3.098E 

01 

-7. 3286-01 

-7 .9856 

-0 1 

2. 90 EE 

Cl 

-1.225 E 

00 

-1 .2056 

00 

2. 385 = 

Cl 

-1.675 E 

03 

- 1 .7 33 c 

33 

2. 8256 

Cl 

-2. 0886 

00 

-2 . 185 E 

00 

3. 187E 

01 

-2.8856 

00 

-2.5016 

JO 

2. 5286 

Cl 

-2. 7966 

00 

-2.8026 

UO 

. 3. 582 E 

Cl 

... -2. 9 96 E 

00 

-3 .096 6 

00 

8.1206 

01 

-3. 1876 

00 

-3 .2306 

00 

8.732E 

01 

- 3 . 32 56 

00 

-3.353 E 

00 

5.280E 

Cl 

-3. 8136 

oo 

-3 .815 6 

00 

8.930E 

Cl 

-3.8676 

00 

-3.8 19 C 

00 

2. 7276 

Cl 

—3.985 6 

00 

-3 .365 E 

JO 

■ 1. 702E 

Cl . 

._ - 3 . 8 80 6 

00 

-3.258 E 

00 

8. 8 72 E 

Cl 

-3.2836 

00 

-3 .10 IE 

DO 

5.983E 

00 

-3. 06BE 

00 

-2 . 09d L 

00 

3.9796 

Cl 

-2. 8 00 E 

00 

-2.6556 

00 

1 . 28 76 

oc 

-2. 5866 

00 

-2.3/76 

00 

2.9336 

01 

-2.225E 

00 

-2 .0 73 fc 

00 

,2. 8266 

01 

..-1. 0 99 6 

00 

-1 .7906 

00 

1 . 109E 

01 

- 1 . 5 5 1 E 

00 

-1 .3 93 t 

uO 

2. 1736 

Cl 

-1.1 B8E 

CO 

-1 .0376 

00 

2. 8836 

ct 

-8. 221 E-0 1 

-o .766fc-0l 

1.8826 

Cl 

-9. 651 fc- 

-01 

-3*1066- 

-0 1 

2. 1186 

Cl 

-1.1156- 

-01 

3.1516- 

-0 2 

2.6886 

01 2. 3C2G- 

■01 

3.6816- 

-0 1 

2.0 5 76- 

01 

5. 9826- 

•01 

6.861 £- 

•01 

2. 6256 

01 

8.8 096- 

■01 

9.8126- 

•0 1 

2 . 2 7 7 E 

Cl 

1.1176 

oo 

1.29/1 

Oj 

2.6685 

Cl 

1. 3606 

00 

1 . 9 8 / fc 

JO 

2.0866 

01 

1. 56/5 

00 

1.6916 

00 

2. 6 92 5 

01 

1 . 7396 

00 

1 . 6 8 J L 

33 

2.9826 

01 

1 . 0 7 7 r 

00 

1 .9 92 C 

JO 

3. 2096 

Cl 

l . 9005 

CO 

2 .Oil ft. 

00 

3.91 9C 

Cl 

2. 09 9 C 

00 

2 * 1 9 ; t fc 

u J 

3. 82 66 

01 

2 . COOfc 

00 

2. 165 L 

a J 

2.9576 

01 

2. C 5 9 E 

00 

2 . 1 99 L 

00 

1 . f 7 95 

Cl 

2.Q01C 

00 

2 . O'/a c 

00 

8. 5836 

00 

2.0206 

00 

2.016c 

0 J 

1.019E 

01 

1 . *J 96 6 

00 

1 .9096 

JO 

2. 98 5E 

Cl 

1, 79 HE 

00 

1 .7 7.5 6 

00 

1.2 976 

01 

1*6005 

00 

1 .60 JE 

00 


— 12 - __ 


Y2 . 


0.0 

0. 6y7fc 01 
8.1&7E 01 
2.208 E 0 1 
8 . 1 1 L 0 0 

-3.722E-02 
-l.o J6£ 01 
- 1.8836 01 
“2 . lu9E 0 1 
. -2 .5786 01 
-2.788E 01 
-2.7986 01 
-2.7/7E 01 
-2 .725E 01 
-2 « 6 85 E 01 
.-2.539E 01 
-2.9C8E 01 
-2 . 2 ui E 01 
-1.901E 01 
-1 ,6tj 3 6 0 1 
-1 -3 /6E 0 1 
„-l . 1 J7E 0 1 
-tf.iudC 00 
.-8 .30 l E 00 
-J.558E 00 
-I.U77E 00 
1.7236-01 
5. 3 50 6 00 

1 . u 9 6 E 01 
9 . 29 3 E 00 
1.8256 01 
I . 9 7 7 6 01 
1.895C 01 

__ U » 8 u 3 E 01 
l.UuOE 01~ 
. .1. I 38 E 0 1 
1 .8296 01 

i.hjoe oi 

1 .7oOb 0 1 
1 .0076 0 1 
1.3396 0 1 * 
1 .'.20 6 01 
1.2 89c 0 1 

1 . 1 29 t 0 1 
9 . »» J 7 fc 00 
/./(ME 00 

0 .0/. 2 6 00 

9 . 2 a 81 00 ■ 

2 .u 76 l 00 

1 .2 -.9 1; 00 
“i .110 E-G 1 
-1 fc GO 
-9 .8 98 I; 00 
-7.9296 00 
“7 .08 9 L Oil 
—8 . 3 7 1 00 


0 .0 

6.702 E 01 
'♦.HOE 01 
.. 2.19 7E 01 
7 . 7 8 2 E 00 
-2.900E 00 
-1-0U7E 01 
-1 . 6 80 E 01 
-2. USE Cl 

— 2 .8226 Cl 
-2.6256 01 
-2.7806 01 
-2.779E 01 
-2.7585 Cl 
-2.672E Cl 
-2.5806 01 
-2 . 3 6 5 E 01 
-2.1586 01 
-1.913E Cl 
-1.6806 01 
-1.366E 01 
-1.070E 01 
-7 *60 2E CO ' 

- 9 . 6 8 0 E CO 

-1.652G 00 

1.237E CO 
8 . 0 3 9 E CO 
..6.6 386 00 
9.032 E 00 
1.U9E 01 
1.3076. 01 
1 .<.67 E 0 1 
1.3976 Cl 
_ 1 .69 5 E 01 
’l .76 26 01 ’ 
1.7976 01 
1 .80 1 E 01 " 
1.775E 01 
1.721E 01 
1.6616 01 
~ 1 . 5 3 7 E Cl ' 
1.8116 01 
1 .26 76 0 1 
1 . 1 O 7 E 01 
9 . 3 3 6 E CO 
7.38 2 E 09 
9.69 3E 09 

3 . 6 nr- oo 

1.939 C 00 
1 .0636-01 
-1.0396 00 
-3.3296 CO 
-8.6826 00 
-6.2 V6E 00 
-7.393E 00 
-0.630E 00 


0.0 

. 1.109E-03 
0 .0 5 IE- 0 3 
_ 2.8 70E- 02 
5 . 3 3 2E- 0 2 
9.897E-02 
1.89QE-01 
2U79E-01 
2 .966E- 0 1 
. 3. 862 E- 01 
8.8896-01 
. 5 . 896E- 0 1 
6.996E-01 
8. 1276-01 
9.2676-01 
..1.080E 00 
1.1 50E CO 
1.2555 00 
1.358E CO 
. 1 . 8 8 6 E 00 
1.S2HE 00 
...1.60 IE 00 

1.662E 00 
. 1. 712E 00 
1.7896 00 

1.7 75E 00 
1.708E 00 " 
_1 . 7895 00 

1.7 795 00 ’ 
1.7876 00 
1.7256 00 
1.6836 00 
1.6336 00 
1.575E 00 
1.5115 00 ’ 
1.882E CO 
1.369E 00 
1.2986 00 
1 . 2 1 8 E 00 
1.1826 00 
*1.0676 00 
9. 9 8 8£- 0 1 
9.2616-01 
8 .6 196-01 
8. 030E- 01 
7.80 36-0 1 
7 . C 8 2 6- 0 1 
6.65 36- 01 
6.3385-01 
6. 1016-01 
5.9826-01 
5 .U60E-0 1 
5.0 5 26- Cl 
5.9196-01 
6.0585-01 
6 .23 16-01 


‘ C. 0 

5. 33 9: - 02 

1 . a 8 OE - 01 

_ 3 . 7 28E - o L 
5. 635E-C1 
.6.0936 -Cl 
1.0225 00 
„ 1.226= CO 
1.8116 00 
. 1.5 70= 00 
1.7C0E 00 
1.7996 00 
1.8656 00 
1.8975 00 
1.897E 00 
1.8655 00 
1. 802c CO 
1.711: CO 
1. 5986 00 

1.8 585 00 
1.29SE CO 
1.119c 00 
9.3036-01 
7.3256-01 
5.2915-01 
2. 2385-01 
1. 2006 -Cl 
“7. 8 925 - 02 _ 

- 2. 6996 - CL 
-8.5005-01 
-6.1685-01 ‘ 
-7.6785-01 
-9.01 35-01 
-1.016E 00_ 
*1.1106 00 
-1.182E CO 
*1.2386 CO 
*1.2635 00 

- 1 . 2725 CO* 
*1.2596 00 

* 1 • 22 8 6 00.' 
-1.1776 00 
*1.1 ICE 00 
*1.0286 CO 
-9.3236-01 
•«. 2 5 5: -Cl 

• 7.O9/E-01 “ 
•5. 8725 -01 

- 8.6CCE-01 

- 2. J06E-01 
•2.0100 01 

• 7.38 05-02 _ 
5. Cl 86-02 
1.6795-01 

2. 7 ill C - Cl ~ 

3. 7926 - Cl 



f.YCl E 

T IME 


. . ft 1 


11 70 0 

1 . 1 19F 

00 

1.C92E 

02 

1 l 400 

1 . 1 38 f 

00 

1 .0996 

C 2 

t 1600 

1. 158? 

OC 

1 .1406 

07 

1 \ MOO 

1.178? 

CC 

1.0/66 

02 

1 2000 

1.148? 

00 

1. 13 56 

C? 

17 700 

1 .ziar 

00 

1 . ( 8 7 6 

C 2 

17 400 

1.738? 

no 

1.1146 

C 7. 

1 ?r.r>0 

) .?57f 

CC 

1 . 1006 

0? 

l ? a no 

1.771? 

00 

1.0896 

02 

1 i*09 

1 . 7 5 7 E 

OC 

1 . 1066 

0 2 

i izr.o 

i.nn 

00 

1 . C 9 1 6 

C2 

1 3400 

1.327? 

CO 

1,084? 

C 2 

1 3600 

1.757? 

Ci c 

1 .087? 

0 2 

1 3P V) 

1.376? 

CO 

1 .0 866 

02 

\ 4 000 

1 . 796 ? 

00 

1 .0806 

02 

14700 

1 . 4 1 6 f 

00 

1 . 0 7 26 

02 

1 400 ■ 

1.436? 

00 

1.0686 

C2 


1.4609 1.4561: CO 1.068S C2 

14 p.no 1 .'4 7 6 F OC 1 -078E C2 

16 0 CO 1.496? 00 1.10 26 02 

162 00 1.515 f 00 1 . 1 2Bc 0 2 

16^00 . 1 - 5 35 E CO 1.167E 02 

1^600 1.555F-CO 1.1 5tE C2 

15600 1 - 575 F CO ) . HOE 02 

16000 1*595? CO 1.099E 02 

1 6 200. 1.614? CC 1.127E C2 

16 400 1 .6 2 4 F 00 1.133E 02 

16*00 1 . 6 5 4 F CC I.110E C2 

1 4 600 1.674? 0 Q 1.1076 C2 

\ 70 00 1-694? GO . 1 * 1 20 E CZ 

172C0 1.714? 00 1.117? 02 

17 400 1 . 7 2 0 E CC 1.106E 02 

17600 1.753? CO 1.1C6E C2 

17800 1. 77 3 F CO 1*1096 02 

IflOPO 1.793F CC 1.1G7E C2 

137*0 1.613? OC 1.102? 02 

18400 - 1 . 83 3? OC 1.C 986 02 

18600 1.F57 F CC 1.094E 02 

18800 1 .P7?F 00, 1. C51E 02 
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1.1236 02 
1. 17 36 C 2 
1.1236 02 
1 . 123E C 2 
1.1236 07 
1.12 36 0 2 


1 . 

12 36 

C 2 

1 . 

12 36 

02 

1. 

1 2 26 

02 

1 . 

1 2 36 

02 

1 . 

1 2 26 

02 

1. 

12 36 

02 

) . 

l 2 36 

0 7 

1 . 

1 2 36 

C 2 

1 . 

12 36 

C2 

1. 

l 23E 

02 

1 . 

12 36 

0 2 

1. 

1236 

02 

1 . 

12 36 

02 

1 . 

l 2 36 

02 

1 . 

12 35 

02 

1. 

1236 

02 

1. 

1 7 3E 

02 

1 . 

.12 36 

02 

1 . 

.1236 

02 

1. 

, 1236 

02 

1 . 

, 1236 

02 

J . 

. 12 36 

0 2 

. 1. 

.12 36 

C 2 

1 . 

.1236 

C 2 

1 . 

. 1 ?3E 

02 

1. 

. 1236 

02 

. 1 

.1236 

02 

1. 

.1736 

02 

1 

.12 36 

C2 

1 

. 1236 

C 2 

1 

.1236 

C 2 

1 

. 1236 

C2 

1 

.1236 

02 

. 1 

-1236 

C2 

1 

.1236 

02 

1 

.1236 

02 

1 

. 1 2 26 

02 

1 

. 1 2 2c 

02 

1 

.1226 

C 2 


B 2 


2.2206 01 

2.2206 Cl 
2.27C6 Cl 
2.7206 Cl . 
2. 2 206 01 

7.2206 CL 
2.220C 01 
7.22CE Cl 

2.2206 01 

2.2206 Cl 
2.220C 01 
2. 7206 01 
2.770C Cl 
2.22 Ct Cl 
2-2196 01 
2.2196 01 

2.719E Cl 
2. 2.1 BE Cl 
2. 21 BE Cl 
2. 21 BE 01 
2.21 'IF Cl 
2.2176 CL 
2. -2 176 01 
2.2176 01 
2.2176 Cl 
2.2176 01 
2.2176 01 
2.2.176 01 

2.2176 Cl 
2.2176 Cl 
2.7176 Cl 
2.2176 Cl 
2.2176 Cl 
2.2176 01 
2.2176 Cl 
2.2176 Cl 
2.2176 Cl 
2.2176 01 
2.2176 01 
2.2176 CL 
2 . 2 1 76 Cl 
2.217E Cl 
2.2176 01 
2.21 7E CL 
2.2176 01 
2.2176 Cl 
2.217E Cl 
2. 2176 Cl 
2.2176 01 
2.2176 01 
2-2176 Cl 
2.217E Cl 
2.2176 Cl 
2.217E Cl 
2.2 176 Cl 
2.2176 Cl 


XI ■ ■ 

-6. *156-03 
-*.2266-03 
-1. 9*9 G -03 
*. 5C9F-0* 
2. 8696-03 
5. 0796-03 
6. 0*96-03 

0. *196-03 
9. 6 9ft ( -03' 

1. 0 72 t -02 
1. 1*86.-02. 
1.2016-02 
1.2306-02 
1.2366-02 
1 . 222 6-02 
1 . 1 ft 7 6-92 
1. 1.33L-02 
1.0*36-02 


9. 

7 75 E- 

03 

8. 

7 9* E-l 

33 

7. 

7 05 6-' 

03 

5 . 

5 32 E~ 

03 

5. 

2 56\- 

03 

*. 

019 6- 

03 

2 . 

7 2 7 F - 

03 

1. 

**K> 

G3 

1. 

7806- 

0* 

-1. 

0**6- 

0 3 

-2. 

1 806 - 

03 

-3. 

160E- 

03 

-*. 

0056- 

■03 

-*. 

7*5003 

- 5. 

* 02 6 * 

•03 

-5. 

,809L“03 

-5. 

,9116- 

-03 

-5. 

, 81*6- 

-03 

-5. 

, 5 92 OH 3 

-5. 

.3006- 

-03 

-*, 

.9716- 

-03 

-*, 

. 6 32 6- 

-03 


.2986-03 

-3. 

. 9 80E- 

-03 

-3 

.6836- 

-03 

-3 ' 

.*116 

-0 3 

-3 

. 1 ** E 

-03 

-2 

. 65*6 

-03 

-l 

• 966 E 

-03 

-1 

. 16* l 

-03 

-3 

. Q65 E 

-c* 

* 

. 7196 

-0* 

l 

. 0326 

-03 

1 

.*196 

-03 

1 

. 6806 

-03 

1 

. 8 3 ft £ 

-03 

1 

.9136 

-0 3 

1 

. 9326 

t 

o 
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•5. 

2*2 6- 

-03 

3 . 

.0756- 

-0 3 

■8 . 

.389 U- 

-0* 

1 . 

>50*6- 

-0 3 

3 . 

,8516- 

-0 3 

5. 

. 925 L-0 3 

7 . 

.667 L 

-0 3 

9 

.1 10 6 

-J 3 

1 ■ 

.0 2/6 

-0 2 

1 

.117 6 

-J 2 

l . 

.182 6- 

-0 2 

1 

.2226 

-0 2 

1 

. 2 39 6 

-0 2 

1 

.2 3*6 

-0 2 

1 

. 2 09 L 

-0 2 

1 

.18*6 

-0 2 

L 

.1016 

-02 

l 

.0236 

-0 2 

9 

.3126 

-0 3 


8 .2766 -0 3 
7 . lWt-J 3 
5.9*16-03 

* .68 7 6-0 3 
3 .*05 6-0 3 
2 . 1 1ft 6-0 3 
ft . * ftl L - 0 * 
-3 . ft 82 6 - 0 * 
-l .5 7 5 L-0 3 
-2 • 66ft 6-0 3 
-3 . 5 9<j t -0 3 
-*.387 6-03 
-5.0*7 6-03 
-5 . 6 6 1 6 - J 3 
-*.0 0*6 —0 3 
-£> .0*16 “0 3 
-5 .ft 0 1 6-0 3 
-3 .600 6-0 3 
-5.252 6-0 3 
-*.07*6-03 
-*.*926-03 
-*.1226-03 
-3.7756-03 
-3 .*56 6-0 3 
-3 .16 7 E“3 3 
-2 .003 E~03 
-2.3916-03 
-1 .699 6-0 3 
-8 .98 1 1-0 * 
— * .*3*6-05 
7.21* £-0* 
1.2666-03 
1.631 6-03 
1 .867 E-0 3 

1 .9 9 3 u-j'3 

2 .0*7 £-0 3 
Z.03- 6-03 


X2 


1.1 186-01 
1 .0 85 6-01 
1.1796-0 1 
1.2136-01 
1.1*9 I: — 0 1 
9.9736-02 

8 . **36-0 2 
/ .06*6-02 
5 . 7*06-02 
*.**56-02 
3 .2 J* 6 -0 2 
2 .0 35 6-0 2 
8.8256-03 
-2.1*96-03 
-1.48/6-02 
-2.24 5F.-02 
-3 . 10*6-0 2 
-3.89/6-02 
-*.59*6-02 
-5. 1 91 6-02 
-5.**96-02 
-*. 0306-02 

— 6 . 30 1 6-0 2 

— 6 .**0 6-0 2 
-* .* 756-02 
-o .3*86 — 02 
-6 .2376-0 2 
-5 .9 ** 6-0 2 
-3.29*6-02 
-*.3376-02 
-3 .9*56-02 
-3 .*806-02 
-3 .0 20C-0 2 
-1.17*6-02 

3.7 1 U-0* 
ft .** 66-03 
1 .3 196-02 
1 .8 76 E-0 2 
1 .*8 3 6-0 2 ’ 
1.692 6-0 2 
1.8356-02 
1.5*06 -02 
1 .*236-02 
1.2976-02 
1 . 7 70 E-0 2 
3.0476-02 
3.7826-02 
*.1*66-02 
*.352 E-0 2 
3.3126-02 
2.32*6-02 
1.589t“02 
1.0 316-02 
5.u*36 _ 03 
2. 1 596-03 
-3 .7*96-0* 


LZ 


1.112E-01 
1.0706-01 
1.15*6-01 
1.181E-C1 . 
1 . 1 306-01 
9.51*6-02 
7 .9*7 6-0 2 
6 . 80 3 E -02 
5.1* flC-02 
3.85* 6-C2 
2- 6221: - 02 
1 .* 3 16-0 2 
2.9826-03 
-7- 71*E-C3 
-1 .7 7*6-02 
-2 • 65*6” 02 
-3.5266-02 
-* .26 U-02 
-* .8996-02 
-5 . * 3 5 E- 0 2 
-5.88 1 G- 02 
-6.1596-02 
-6 .3606-02 
-6. **06-02 
-*.*206-02 
-6.2616-02 
-6.082 6-C2 
-5. 7 * 7 6- C 2 
-5.060E-C2 
-*.2716-02 
-3.b5*E-02 
-3.l72F.-02 
-2 .7096-02 
-ft .5136-03 
3* 7 7 6 E-0 3 
1.15*6-02 
1.8086-02 
1.8386-02 
1.91*6-02 
1.3096-02 
1.7986-02 
1 .670E-C2 
1.5226-02 

I. 3 67E-02 
1 . 8 1 6E-G 2 
3 .0536-02 
3. 7686-02 
*. 1 71E-C2 
*.3126-02 

J . 2* 5E-C2 
2.2 30E-02 
1* *726-02 
8.932 E-03 
*.*266-03 
7 .07 6 E- 0* 

-2 .U05E-C3 


Y 1 


1. 

C026 

00 

1. 

0026 

00 

1 . 

00 1G 

00 

1. 

0016 

00 

1. 

0006 

00 

9. 

9976- 

01 

9 . 

99 36- 

01 

9 . 

>9 9QE- 

01 

9. 

9 8 76- 

OL 

9 . 

,90*6“ 

0 l 

9. 

,9826- 

01 

9. 

, 980E-0V 

9 . 

.9796- 

■01 

9-. 

.9786- 

-01 

9. 

.9776- 

'01 

9. 

.9776- 

■01 

9, 

. 9 7 76 - 

-01 

9 

.9 786- 

-01 

9 

. 9 7 9 E - 0 1 

9 

.9H96-01 

9 

.98 26-0 l 

9 

.9 8*6 -01 

9 

. 9066- 01 

9 

.9896' 

-01 

9 

.99 IE- 

-01 

9 

.99*6 

-01 

9 

.9966 

-0 1 

9 

.9996 

- 0 1 


l.OOOE 00 
1.0006 00 
l.OOOE 00 
1.0016 00 
l.ooie oo 
uooie ca 
1.0016 00 
1 .00 IE 00 
1.0016 00 
1.0C16 00 
1.0016 00 
l.ooie no 
1.0016 00 
1.0016 00 
1.001E 00 
1 .00 16 00 
1.001E 00 
l.ooie oo 
1.0006 00 
1.0006 00 
1.0006 00 
l.OOOE 00 
9. 999E-01 
9.9996-01 
9.999E-01 
9 .9 9 BE ” 0 1 
9 .9936-0 1 
9 . 9 9 H k - 0 1 


. Y2 ... 


-6.3Z6E-03 
-6.5756-03., 
-C. 6 926-03 
-6.6726-03 
-6.5106-03 
-6.21*6-03 
-5.0 0*6-03 
-5.2596-03 ,, 

- * . 7 L 6E - 03 

- *.a 71H-03 . 

- 3. 380E-03 

- 2. 65 7E-03 . 

- U 9 L HE - 0 3 
-1.1 786-03 

- * . * 0 76 ” 0 * 

2 . 7 1 8: - 0* . 
9. 8206-0* *’ 
1 . 590G-03 
2.177F-03 

2. 7C5E-03 . 

3. 1 *76 -03 

3 . 56 CE - 03 
3. 8796-03 
*.1216-03 
*.2 86!: - 03 
*.3 756 - 03 . 
* . 3 ft BE- 03 

4. 32 HE - 03 
*. 1 996-03 
*.0106 -03 

3. 7 70E- C3 
2. *0t6-C3 
3- 1636-03 
2.81 IE-03 
2. **96-03 

■ 2.0906-03 

1 . 7 * 56 - 03 
1. *206-03 
1. 116E-03 
6.3*96-0* 

" 5.7666-0* 
2. 3996 - C* 
1.232E-0* 

- 7. 532E-C5 
-2.5726-0* 
-* . 16 7F - 0* 
- 5.*C0£-0* 
-6. 182E-0* 
-6. *6 56-0* 

- 6 . 2 * y 6 - 0 * 
-5.6*26-0* 
-*.7656-0* 

- 2 . 7096 -0* 
-2. 5**6 - C* 
-1. 32 66 -C* 

- 5 . ft* IE-06 



CYC.lE • 


67700 < 

67 403 t 
67600 < 

67R00 — < 
68000 
68700 

400 

fcftfcGO • 

hH HO 9 
60 000 . 
607 00 
69403 
60 fc 00 
6 ) 7.00 
70000 
70 700 
704 00 
706 00 

70 P 00 
7 1000 
71700 

71 4 00 .. 

71 f>00 

_ 71 8 CO . 

77000 
_ 77700 ... 
77400 

„ 77600 _ 
7 7 800 
.. 7 30 00 
73700 
.. 734(;0 
73600 

7 3 800 . . 

74000 
.. 74 7 00 . 
74 400 

.. 74 ft no 

74 AOO 
7 60 00 .. 

75700 
75400 , 
75600 

75 BOO 

76 0 00 
_ 76? 00 

76 400 
76600 

76 800 
7 7000 
7 / 700 

77 400 
7 7 600 
77 B 00 
78000 
78200 


T ImE 


ft 1 


>.fr73F CO 
b.fcSU CO 
b • 7 1 2 E CC 
b.732E 00 — 
6 » 7 5? f CO 
6 • 7 72 F CO 
6 . 7 9? E CC 
6 • 01 2 E 00 

6.831F CC 
6 . fl 5 l f 00 — 
6.B71E 00 
6.8916 00 
6.911 F 00 
6.9216 00 
6.950F 00 

6 • 9 70 f 00 
6 . 990 f CO 
7.0106 00 

7.0 30 f GO 
7.050E CO 
7.C69E 00 
7.089 E 00 

7 . 109 f 00 

7 . 1 29 6 DO 
7 . 1 49 E CO 

. 7 . 1 6 9 E OC .. 

7 . J 8 ft F 00 
_ 7. 2086 CC _ 
7.7786 CC 
7.7486 OC 
7.?6ftf 00 
7 . 2 ft ft E 00 
7.307F 00 
__ 7. 3 2 7 E CO.. 

7.341E CC 
. 7.367E 00 
7.3ft / f C Q 
7.4076 CO 
7.477E CC 

7 . 4 4 6 E 00 

7 ♦ 4 6 6 F 00 
4866 CC 
. 50 6 E CO 
, 5 2t» 6 00 
546 F 00 
, 5 6 5 E 00 
, 5 ?. 5 f CC 
.605 F OC 
.6756 00 
645 F CO 
,<656 CC 
,6846 C C 
.7046 00 
.724E 00 
7.744E CC 
7. 764 E 00 


7. 

7. 

7. 

7. 

7. 

7< 

7 

7. 

7. 

7. 

7. 

7 

7 


1.1226 
l .1236 
). 123E 
1. 123E 
1 . 1 2 IE 
1. 1 23E 
1. 1736 
12 36 
1 2 it 
123E 
1 2 31? 

1 . 12 311 
1.12 2 ? 

1 2'3t 
1236 
1 2 3 B 
1 2 26 
1 2 26 
1 . 1 2 3E 
1.12 36 
1 . 1 236 
) . 1 2 26 
1 . 1 2 36 
1 . 1236 
1736 
) 2 36 
12 36 
1 2 3 E 
1236 
1236 
1236 
12 36 
1 2 36 
. 1236 
. 1 2 31: 

. 1 236 
.12 36 
. 1 2 36' 
.12 36 
. 1? 36 
. 12 36 
12 36 
1 2 3F. 
1236 
12 36 
1 2 26 
l 7 36 
12 36 
12 36 
,12 36 
.1236 
,1236 
. 1236 
.1236 
. 12 JE 
.1236 


1 

1. 

1. 

. 1. 
I . 
1. 
1 . 

. 1. 
t . 
. 

1 . 
. 1. 
1. 
1 . 
1. 

- 1 ■ 

1 

1. 

1 

1 

1 


C2 
02 
C2 
C2 
02 
02 
02 
02 
C2 
02 
02 
02 
02 
C? 

C? ’ 
C 2 
0 2 
02 
02 
C 2 
C 2 

0 2 ... 
02 

0 2 . 
02 
C 2 . 
02 

02 _ 
0? 
02 , 
02 
02 . 
C 2 
02 .. 
02 
02 
02 
■ 02 
02 
02 
02 
02 
02 
C 2 
02 
02 
02 
C2 
C 2 
02 
02 
02 
02 
02 
C 2 
02 


B2 


2-2176 01 
2.21 76 01 
2.2176 01 
2. 2176 01 
2.2176 Cl 
2.2176 01 
2.2 176 Cl 
2. 21 76 
2.2176 


XL 


Z1 


X 2 


22 


Yl 


Y2 


Cl 

01 


Cl 

Cl 

01 

Cl 


2.2176 Cl 
2. 21 76 01 
2.2176 01 
2.2 1 76 01 
2.2176 01 
2. 2 176 Cl 
2.2176 
2. 21 76 
2. 21 76 
2 . 2 1 76 
2.2176 Cl 
2.2176 Cl 
2.2176 01 
2.217E 01 
2. 21 7c 01 
2.2176 Cl 
2.2176 Cl 
2.2176 01 
„2 . 2 1 76 01 
2.2176 01 
2.2176 01 
2.2176 01 
2.2176 01 
2 . 2 1 76 
.2-2176 
2 . 2 1 76 
2.2176 01 
2.2176 Cl 
2.2176 CL 
2.2176 01 
_ 2.2176 
2176 
2176 01 
2176 Cl 
2176 
21 76 
21 76 
2176 Cl 
2176 01 
>2176 01- 
,2176 CL 
.2176 Cl 
.2176 01 
.2176 Cl 
21 76 01 

2 1 76 Cl 
2176 Cl 


2 . 
2. 
2 . 
2 , 
2 - 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 , 
2 . 
2 
l. 


01 

ci 

ci 


01 

Cl 


Cl 

Cl 

01 


1. 09? 6-03 
1. 8096-03 
1. 6876-03 
. U 5386-03 . 

1 . 362 6-03 
l. 1696-03 
9. 6146-04 
7. 4 00 6-04 
5. 1 516-04 
3. 0466-04 . 
1.4/5 1-04 

5. !>15 C-3 5 
-1 . 08 t 6-05 
- 5 . 3 4 8 L - 05 
-7. 8886-05 
-9.2 32 6-05 
-9. 740 6-04 
-9.746 6-05 
- 9 . 3 84 6-05 
-8. 8126-05 
-8.1266-05 
-7. 3926-05 
-6.6546-05 
-5. 94 0 6-03 
-5. 2666-05 
— 4. 644 6 — 05 
-4. 0766-05 
_-3. 563 6-05 
-3. 1 C/6-05 
-2. 702 L- 05 
' -2.3436-05 
-2. 0296-05 
-1. 7546-05 
-1. 3136-05 _ 
-1 .303 6-05 
-1. 1241-05 
-9. 6 67 6-06 
-8. 3 1 l 6-04 
-7. 1406-06 

-ft.. 131 6-06.. 

- 5 . 2616-06 

-4.5136-04 

-3. 0 /0fc-06 
-3. 317 6-06 
-2. 84? F-Ob 
4 33 6-0 6 
2 .0 84 6 -06 
1 .7066-06 
•U 5296-06 
•1. 309 6-06 
-1. 1206-06 
- 9. 5 836-07 
-8.2026-07 
-7. 01 8 6-07 
-6. 0056-07 
-5.137E-07 


- 2 . 


119 6 4 6-33 

L .853 6-3 3 
1 . 7116-03 - 

1 .540 6-0 3 , - 
1.345 L-0 3 - 

1.1346-03 - 

9.158 6-0 4 
6.3 53 6-34 
4.3396-04 - 

2.4126 -0 4 
9 .3 8 26-0 3 ■ 

-V .5 336 -J 6 
-7 .3796-05 • 

-l . L 22 6-04 
-1 .3 266 -0 4 
-1 .405 6—0 4 
-1 .4 0 2 6-04 
- 1 .3 4 46 -0 4 
-1 .2556 -0 4 
-1 . 149 6-04 
-1 .0 35 6 -O'* 

-9 .2236 -0 5 
-il . 1 30 6-0 5 
-/ .11 16 -3 5 
-6.178 6 - J 5 
-5 .3 3^6-0 5 
—4 .5856 —0 5 
-3 .923 6-0 5 ... 
-3 .3446-05 
-2.8416-05 
-2 .406 6-0 5 
-2.0 326-05 
- L . 7 1 2 6 -0 5 

-1.4'39fc“05 . 

-1 .7.0 / L-0 5 

-1.0116-05 

-8.4 4 9 6 -0 6 
-7 .0 5 2 6-0 6 
- 3 . 8 / / 6 -06 
_ -4 .89 2 6-0 6 
-4 .04 6 6.-0 6 

-3 . 3 /i> 6-00 
-2.800 6-06 
-2 .3206-06 
-1 .3 2 16-0 6 
- 1 . 338 6 -0 n 
-1.313 1-06 
-1.0846-06 
-8 .940 6 -0 / 
-7 .371 6-07 
-o .0/36-0/ 
-5 . 0 0 1 l -9 / 
-4.1 136-0 7 
-3.3846-0 7 
-2 .7826-0 7 
-2 .286 6 -0/ 


3.244 6-03 
5 .1126-0 3 
■6 .8 ?9 E -0 3 
•8.160 6-03 
-9 .3 60 L-0 3 
*9.8 366-0 3 
*1 . 0826-02 ” 
-1.1 226-02 - 
-1 . 1 18 6-02 - 
-8 .6 33 6-0 3 
-6.1366-03 - 

-4.1186-03 “ 

-2. 0746-03 - 

-1. 692 6-03 - 

-9.3366-04 
-4.4246-04 * 

- 1. 1 016-04 
1 . 0 72 6 - 04 
2. 426L-C4 
3 .207 6-0 4 
8.5926-04 1 

8.7096-04 , 

3 . usO 6 -04 
3. 4 806-04 
3 .245 6-04 
2.5776-04 , 

2 .6976-04 
2. 4206-04 .. 
2.1556-04 
1 .90 7 6-04 
1 .6 /ft E-G4 
1.47LE-04 
1 . 2836-04 
1.119 6-0 4 .. 
9 .7^06-05 
8.4266-05 
7 .290 6-0 5 
0.2986-05 
5 .4 34 6-0 5 
4 .6826-0 5 
4 .0 3 16-0 5 
4 o 7 6-0 5 
9806-05 
5 oQ 6-0 5 
1 9 7 6 - C 5 
l . 8 u5f? - C3 

1.0 17 6-0 5 
1 . 3ft;> 6-0 5 
1 . idilL-OS 

1 .0 ItiE -0 5 
7206-06 
4oft t-C 6 
3946-06 
4 /4 6-0 6 
, 68 6 E -06 

4.0 106-06 


4.6106-03 9.9986-01 

6.3806-03 . 9. 99 BE -0 1 

■8.006E-03 9.9986-01 

•9.1686-03 9.9986-01 

•1 .022 E- C2 9.9996-0 1 
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-1 .9 9 16- 13 

2.21 76 

Cl 

-1.0496-11 

-1 .615 E-13 

2.2176 

01 

-8.9 74 £-12 

-l .3 106-13 

‘2 . 2 1 7E 

Cl 

-7. 6736-12 

-1 .062 £-13 

. 2.2176 

Cl 

-6. 561 E-12 

-8 .6 11 E- 14 

2 . 2 1 76 

Cl 

-5.611 6-12 

-6 .9 8 It- 1 A 

2.2176 

01 

-A. 798 6-12 

— 5.6 dO h — 1 A 

2. 21 76 

01 

-A. L 02 E-12 

-A .5 M 9 L - i 4 

2,2176 

01 

-3. 5086-12 

-3 . 720 £-14 

2.2176 

Cl 

-3.0006-12 

-3 .0 166-14 

__ 2 . 2 1 7E 

01 . 

-2. 365 6-12 . 

-2.444E-14 

2,21 ?G 

Cl 

-2.1916-12 

’* -1.98 1 t- 14 

... 2.2176 

01 

-1.0756-12. 

_ -1 .606 E-LA 

2.217C 

01 

- 1 .604 6-12 

~ -1 . 302 6 — 1 A 

2.2176 

Cl 

-1.3716-12 

-1 .055 £-14 

2.21 71: 

01 

-i. i m-i2 

-8 .6 AML— 15 

_ 2.2176 

Cl . 

„-l. C03r-12 . 

.... -6 .92/ 6- 15 

2. 2176 

Cl 

-tl. 3 726-13 

— 5 ,6 1 i L~l5 

2 . 2 1 76 

01 

- 7. 3 306 -1 3 

-A .5AM L - L 5 

?. 2176 

01 

-6. 260 P-1 3 

-3 .685 6-15 

2.2176 

01 

-5. 3606 - 13 

-2 . 9 M 6 h - L 5 

2. 21 ?£ 

CL 

— A. 58 4 6-13 

-2. A 19 L — 1 5 

_ 2 . 2 176 

CL 

-3. 9196-13 

_ -1 .960 c - 1 5 

2.21 7E 

01 

~ 3 . 3 5 L 6-13 

-1.5MMC-13 

2. 21 76 

01 

-2.8656-13 

-1 . 2 M 6 L — 1 6 . 

2. 2176 

01 

-2. 450 E-13 

-1 . 0426-15 

2.21 76 

01 

-2.Q95E-13 

-tl . 439E-L 6 

2.2176 

Cl 

-1. 7926-13 

-6 ,H 35 E-lo 

_ 2.2176 

01 

-L. 5326-13 

_ —5 • 5 3i» t — lo 

2.21 76 

01 

-1.3106-13 

' -A .AbAt-16 

2. 2176 

01 

-1.1206-13 

-3.631 E— 16 

2.2176 

01 

— 9.3 /ME - IA 

-2 . 9 A L fc - 1 6 

2.21 76 

Cl 

-0. 1 9 0 6 - 1 A 

-2 .3626-lft 

2 . 2 1 7E 

Cl 

-7. 003 £-14 

-1 .9 29 fc - 1 6 

2.2176 

Cl 

-5. 9UHE-L4 

-l .562 L-lO 

2. 21 76 

01 

-5.1206-14 

-1 .Zo'j t— 16 

2.2176 

01 

-A . 3 7 8 6 - 1 A 

-1 .0246-16 

2.2176 

Cl 

• -3. 7 A A E - 1 A 

-8 . 2 92 L - 1 7 

2.2176 

Cl 

-3. 201 6-1 A 

-6-7146-17 

2. 2L7E 

Cl 

-2. 7376-1A 

-5 .A 36 6-17 

2. 2176 

01 

-2. 3A1 6-1A 

-A .401 £-17 

2.21 7E 

01 

-2. 002 6-1 A 

-3 .563 £-17 


L2 


Y1 


Y2 


A l 


5 7d E - 10 2.0 72 6-11 

4 .598 E- 10 1.6826-11 

3 .9 32 t- 10 1 . 3 6 6 6- 1 1 

3 » Jo2.E“ 10 1.109E-L1 

2. 8/56-10 9.000 E- 12 

2.430E-10 _ 7 .306 E- 12 

2.1026-10 5.9296-L2 

1.7976-10 4 .0 1 2E- 12 

1 .6 J7fc- 10 3 .905 E- 12 

1. 31 At- 10 .. 3 . lc.9 £- 12 
l. 12A6-10 2.371 E- 12 

9.6106-1 1 ... 2.0M6 E- 12 

3 .2 17 E - 1 1 l.o92E-12 

7.026E-11 1-3736-12 

6.0GBE-11 1.1136-12 

5 .1376-11 9.0306-13 

A . 393 E - 1 1 7.3246-13 

3 . 7 !> 6 £ - 1 1 5. 9406- 13 

3.2126-11 A 17 E- 13 

2 .740 t-1 1 „ 3.9066-13 

2.3AHE-1 1 3.Lo7E-l3 

2.0uHC-ll .. 2.56RE-13 
1 . / 1 76- 11 " 2.0 8? E-13 
; .'.{.BE- 1 1 1 ,6886-13 

1 ..'556- 11 1.3686-13 

l.G/AE-11 1.1096-13 

v.luOL-12 U.V906- 14 
7. toO E- 12 .... 7.2117 6- 1A 
6 .7 1 2 E - l 2 S.VOoC-14 
5 ./ 4 3 l - 1 2 A . 7 b 0 E - 1 A 

A . 9 U 11 1 - 12 3 - M 7 9 6 - 1 A 

4.1966-12 3 . I A A E - 1 A 

3 .5M8 fc- 12 2.5476- 1A 

1 .Ottifc- 12 _ 2 .06 A E- 1 A 

" 2 .02 A E-L2 1 .6726- 1A 

2 - 2 A 3C - 1 2 1.3656- 1A 

1.9 let- 12 1 .090E-1A 

1 . 0 AO E- 12 _U .09 56- 15 
I .A0 3E- 12 “ 7 .2066-15 
1 . 159E- 12 5 .8 376- 15 

i .O 26 E - 1 2 A. 72 BE -15 

M./096-12 3. 0 306-15 

7 . A 9 9 fc - 1 3 ' 3.1026-15 

6.4 12 E-13 2.3136-15 

5. 4 8 36-13 2 .0 3 5 E- 1 3 

4.6896-13 1.6 486-15 

A .U09E-13 *' 1.3356-15 
3.4286-13 1.0UIE-15 

2.9 31 E-13 0.756E- 16 

2.5U6L-13 7 .090 G- 16 

2.1A3E-13 5.7A16-16 

i .8 33E- 13 A .OA96- 16 
1 .567 E- 1 3 3 . 764E-16 


l.OOOE 00 -i.612E-'9A 

l.OOOE 00 6.612E-0A 

l.OOOE 00 £. 6126- 04 

1 .0006 00 6.6126-04 

1 .00 OE 00 t. 6126 - OA 

_ 1 .0006 00 .. 6. 612E-CA 

l.OOOE 00 6. 61 2 1 - OA* 

l.OOOE 00 . 6.612E-0A 

l.OOOE 00 6.6126-CA 

_ l.OOOE 00 6 .612E-0A 

l.OOOE 00 6.612E-0A 

1 .0006 00 . 4. 612; -OA 

l.OOOE 00 6. 612E-CA 

l.OOOE CO 6.612E-0A 

l.OOOE 00 6.612E-0A 

l.OOOE 00 £.6126-04 

l.OOOE 00 £.6125-04 

l.OOOE 00 .. 6.612E-04 

l.OOOE 00 6.6L25-04 

l.OOOE 00 6.612E-0A 

l.OOOE 00 6.612E-0A 

l.OOOE 00 ... 6.6126- OA 

1.0006 00 6. 61 26 - OA 

l.OOOE 00 6.612E-0A 

l.OOOE 00 6.612E-04 

l.OOOE 00 6.612E-0A 

l.OOOE 00 6.6126-04 

_ 1 .OOOE 00 £.6126- OA 

1.0006 00 £.6126-04 

1.0006 00 6. 6126-04 

1.0006 00 6.612E-0A 

l.OOOE 00 6.6126-04 

l.OOOE 00 6.612E-0A 

l.OOOE 00 _ £.6126-04 
l.OOOE 00 6.6L2E-0A 

1 .OOOE 00 £ . 6 1 2 6 - OA 

l.OOOE 00 6.612E-0A 

1.0 COE CO 6. 612E-0A 

l.OOOE CO 6.612E-0A 

_ l.OCOE 00 6.612S-0A 

l.CCOE 00 *" 6.612E-0A 
l.OOOE 00 6.612E-CA 

l.OOOE 00 6 . 6 1 2c- OA 

1.0006 00 6 . 61 26 - OA 

l.OOOF 00 6.6126-CA 

1.0006 00 6.612E-0A 

l.OOOE 00 ' 6.612E-0A 

l.OOOE 00 6.612E-0A 

L.OOOE 00 6.612E-CA 

1.0006 00 6.612E-0A 

l.OOOE 00 £ . 6 1 2£ - CA 

l.OOOE 00 6.612E-0A 

l.OOOE 00 6.612E-0A 




DATA l 1,Z2, VI, Y?/ 0.0,0. 0.0.0, 0.0/ 
- — RE W I ND " 


RE AD ( 5 » i 1 A 1 » AO i 8 l » BO 


—1— FORM AT < A-G 1 0 * 3 I 

WR i T E ( 6 1 1 )Al»AO»Bl f BO 
*=1*0* 


4 


-c 


DT = 1 • E-4 

■ 

DO 100 J = 1 1 100000 
OBSERVABLE 


FORM -WITH -FEEDBACK 


Y 1 D= C * Y 2 

— Y2D=C*Zl 

Z1D = -A1*Z1 + Z2+P1MR-Y1) 

-Z-20=-A0*Z 1 + B0*( R-Y 1 I 


Y1=Y1+DT*Y1D - ' 

Y2=Y2+DT*Y2D 

__ -£1 = ZH-DT*Z ID " 

Z2=Z2+DT*Z2D 

— W RlTE(9-t2)JtRfZltZ2tYl»Y2 

2 F0RMAT(16» FA* L, AE15.7} 

100- CONTINUE 

END FILE 9 

REWIND 9 

STOP _ 


FORTRAN IVC LEVEL 20 


MAIN 


CATE « 72 3^6 


01/55/23 


PAGE 0001 


cccr 

CCC2 
CCC 3 

ccc*. 

CCC5 
CCC6 
C C C 7 ‘ 
CCC 8 
CCC9 
CC 1C 
CC1 l 
C C 1 2 
CC 1 3 
CCU 


START PRCGR am ““7 TT “”T — T 

IMPLICIT REAUK;NI ■ ’ 

INTEGER MOD 
REWIND S 

CATA 21 »Z2 . X1,X2,Y1 ,Y2/0-0*0.0 f 0 .0 ,0 -0 ,0 .0 .C .0/ 

P E AC < 5 » 2 0C >M»F0.C1 ,GC 

F E AC ( 5 » 2 C 1 )Mll»M22tM2 

*REAC{5,2C2)Ml»N22fM2"' 

PEACl 5 » 2C3 K11,C22»C12 . „ ... 

P E A C ( 5 » 2CA ) CT 

WRITE! C»2C)M, HO ,G1*G0*CT 

VRITEtt, 25 >M11 ,M22,Ml2 % 

WRITEIC.21)ML»N22»N12,C1L*02 2,C1J 

■WRITE(6,22) 

1*0.0 .. 


CC 1 5 

CC16 

* C C I 7 

C C l 8 

ecu 

CC 2 C 
CC 2 1 


CCZ 2 
CC 2 3 


CC IOC K * 1 » ECCOG 

IFCK.EC.l.CR.MCD(K,lOOJ.EC.O>WRITE(6*23)T,Hl*HC»''l»GO*Zlf22»Xl,X2* 

1Y l ,Y2 __ . 

T * T -♦ D T 

C REAOING OF REFERENCE SYSTEM 

REAC(S,1CCC)J,R,Z1.Z2.YI,Y2 
iCCO FORMAT ( 1 6 » F A . IUEL5.7) 


IF [ J.EC. 1>G Wl/CT 
IF (J.EC.l)GC*Z2/CT 

MODEL 

>IC»-M«X1 + X2 + GIMR-YH 
X2G*-FC’*X1+CC* ( R- Y 1 ) 


SYSTEM STATE EOJATICNS- 


CC2* 
CC 2 5" 


X 1 * X 1 C T * X 1 D 
“X2*X2*CT«X2C 


C C 2 6 
CC27 

CC 2 8 
C C 2 **' 


E1-Z1-X1 
E 2 ° Z 2“ X 2 

MC»-XlMMll*EUMl2*E2)/Nlt 

hOC a -Xl+{M12*El+M22*E2)/N22 


C C 2 C 
CC 3 I 
CC 2 2 
C C 2 3 


CC 2 4 


CC 3 5 


CC36 

CC37 

CC28 
CC 39 
CC4C 
CCAl 
CCA 2 


H-FHCUFIC 
fo«=hq*ct*i-oc 

FORMAT CM LIAPUNCV STABILITY MODELLING OF SECOND OROER SYSTEM*//// 
-- 1* STARTING PARAMETERS CF MC0CI... , /TI0,»H1» ’ f IPM0.3,T30 t »H0« ME 
210.3, TCC.'Gl* ME10.3, TCC, *C0» ME10.3///' DELTA T « . * r E 1 0 . 3 / / ) 

21 FORMAT I* MATRIX ELEMENTS.. *»•»*//* Nil* *,lPEIC.' , ,1.0X, , N22 o »ELO. 

3/’ CLh MEIC»3,1CX,*C2 2* • , E ' 0 . 3 , 1 0 X , ’ C 1 2» ME 


ICO 

20 


12 , 1CX, *M2* ME1C- 

22^FORMAT( ■ ^ T IME • , 9 X , *H 1 * ,1X, * HO * « 9 X » 1 G 1 * ,9X, 'GO* , 9 y » Ml* » 5X , M2* ,9Xt 
1*XIM5X,*X2MSX,*Y1M5X,*Y2'//) 


23 FORMAT ( T l , m IX , 1PE 10 .3 )) 

25 FORMAT t * M MATRIX *//* 

1M.1 2= * ,E 10 .2//) 1 
2CC FORMAT t AF10.2 ) 

2 C 1 FORMAT t 3F 10,3) 

“2C2 FORMAT ( 2F 10-3 ) 

2C3 FORMAT (2F10. 3) 

2CA FORM AT ( F 2C .7 ) 


Mil- « t lPei0.3il0X,«*22* * , E 10 .3# LOX, • 


\D 



FCfURAN IV G LEVEL 2C 


MAIN LATE - 723-6 

EKC CF PRCGMAM 


01/55/23 


PAGE 000? 


CCA3 
CCA A 
C C A 5 


p Gw INC 5 
S TCP 
E NO 



tliPINCV STABILITY HCUELLING CF SECOND CREEP SYSTEM 


STARTING PARAMETERS CF MCDEL . * . 

V-l«* 2.CQCE CC. FC = 2.C00E C2 


' G 1 “ 0*0 


GO a 0.0 


DELTA T« 1.0CCE-C4 


H MATRIX 

Hll* 7.C0C6 C 2 


M 2 2 * 3.C0CE 01 


M 1 2 * -l.CCOE C 2 


MATRIX ELEMENTS. 


Nil’ 

Oil 


5.CCCE-C1 
l.CCOE OC 


N 2 2 
Q2 2 


5.CCC6-C2 
1 . CCCE 00 


N 1 2 1 
Q 1 2' 


o.c 

0.0 


TIME 


C.C 

9.9CCS-C3 
1 .S9CE-C2 
2. 99C E-C2 
3-95C6-C2 
4.S9CE-C2 
5 .SSC t-C2 
6 . 5 9 C E-C 2 
7 . S il 5 e - C 2 

<; . <; f e E - c 2 
i.csse-c i 

1.155E-C1 
1.295E-CI 
1 . 2S5E-C 1 
1.4S5E-C 1 
1.5596-01 
1.65SE-C l 
1 .7S8E-C 1 
1 .£S£ E-C 1 

1 . 5 5 £ E - C 1 
C 5SE-C 1 
1S8E-C 1 
2SU-C l 
3S8E-C 1 
4 9 2 E -C L 
55E6-C 1' 

45 e s-c i 

2.752c -Cl 

2 . esee~ci 
2 . 5 5 2 E-C 1 

3. C52E-C 1 
3. 1586-Cl 
3.2SEE-C1 


HI 


HO 


;ccor 

,5’EE 
. 5 A 4 E 
. C 7 4 E 
.457E 
. 9 2 C E 
. 5 C 9 E 

. i ( -j t 

!<it7E 
. J79E 
. 2 ?2C 
.5C9E- 
. 3C 5 E 
.6 1 3E 
.CSCE 
.5 3CE 
.C75E 
. 5E1E 
.2 1 3E 
. e 3 3 E 
. 7 C 0 E 
. 5 1 1 E 
. 1C IE 
.C2 IE 
. 7 6 V E 
. 5 3 fi E 

. 2 e 7 G 

.6336 
. 7 A l E 
.6 let 
. 3 ICE 

.42 3E 

,5cce 

. 2 7 9 E 


cc~"2: 


CO 
CC 
Cl 
C l 
Cl 
Cl " 
CO 
Cl 
Cl 
C 1 
•Cl 

or 

C l 
C 1 
CO 
CO 
Cl 
Cl 
C l 
C l 
C l 
CC 
CC 
CC 
CO 
C 1 
Cl 
Cl 
C l 
01 
CO 
CO 
Cl 


1. 
1. 
1. 
0. 
2. 
1 . 
1 . 
1 . 
7. 
1 . 
1. 
1 . 
9 . 
1 . 
1 . 
1. 
1. 
e 
8 
1 
1 
1 
l 
1 
1 
1 
1 
1 
1 
I 
l 
1 
l 


occr'c'2 

5856 C 2 
8 3 2 £ 

3 5 5 E 
93 5 E 
445 6 
. i 5 p r 
.479 E 
. 1C2E 
. 0 P C E 
. 1 4 C E 

e 

. 3 5 4 E 
.22 16 
. I 2 5 E 
.9276 
. 6 2 (: E 
. 2 C 3 E 
.4416 
.549E 
.037 6 
. 25 26 
. 43 76 
. 5 5 3 6 
.9786 
.92 7£ 
.4466 
. 3 R C 6 
.35R6 
.3876 
. 4 6 4 E 
.969 6 
.6486 
. 0 C 8 6 


C 2 
C 2 
C l 
C 1 
C 2 
C 2 
C 2 
C 1 
C 2 
C 2 
C 2 
C 1 
C 2 
C 2 
C 2 
C 2 
C l 
C 1 
C 2 
C2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C? 
C 2 
C 2 
C 2 
C 2 
C 2 


G 1 ^ 

— 

GO 

— 

71 


0.0 


~0.0 


0.0 

* 

1.0206 

02 

1 .E70E 

C2 

9.1506“ 

01 

I.C2CE 

C 2 

1 .870E 

02 

1 . /; 6 4 6 

CO 

1 . C2C E 

02 

1 .8706 

C 2 

2-2 t 26 

00 

1.C206 

02 

1.8706 

02 

2.7316 

0 0 

1 . C 2C G 

C 2 

1.8706 

C 2 

3 . C 6 9 6 

CC 

"• 1 . C 2 C 6 

C 2 

•* 1.87C6 

C 2 ' 

3.3516 

OC 

1.0206 

02 

1.8706 

C 2 

3.5326 

OG 

1 .0206 

02 

1.8706 

C 2 

3.6416 

on 

l . C 2 C I: 

02 

1.8/06 

C 2 

3 . it 9 0 C 

CO 

1 .C2C6 

C2 

1 . 8 706 

C 2 

3.0656 

OC 

1 .0206 

02 

1.8706 

C 2 

3 . 6 2 4 6 

OC 

’ 1.C20C 

02 

1.870 E 

0 2 

3.54 31: 

no 

1.C206 

C 2 

1 . 8 7 0 e 

,C 2 

3.4186 

C 0 

! 1 .C 206 

02 

1.8706 

C 2 

3.2046 

00 

’ 1.C2CE 

02 

1 . 8 706 

02 

3 - C 6 4 6 

00 

’ 1.0206 

02 

1 .8706 

02 

2.GE2E 

CO 

! 1.C20E 

02 

1.8 70 6 

C 2 

2 .0 t 26 

CC 

. ” 1.0206 

C 2 

' 1.870E 

C 2 

2.4266 

OC 


?2 


'0 - n 

l.-’OOE 00 

1 . r B 7 6 0 0 
l.^NCE 00 
-1 . r 2lE-0l 
* 1 .*’826 00 

- 3 . * 9 1 E 0 0 

- 6 . ^ 0 9 6 0 0 

00 

- I . 1 * 6 G 0 1 
-1/716 01 

- 1 .'9 7E 0 1 

- 1 . n 7 o e oi 


XI 


X2 


Y1 


Y2 


0.0 ‘ 0.0" 
1.0056 00 0.6576-01 

1.9476 00 -1.2616-01 
2.5036 00 - l - 9 2 1 E 00 
2. 4786 00 -2.0496 00 
2.3776 CO -3.0046 00 
' 2.65 76 00 -3. DOGE 00 
3. 6086 00 -5.8826- 00 
4.2676 00 -9.4856 00 
3.6766 00 -1.1126 01 
3.28C6 00 -1.2306 01 
3.7436 00 -1.5666 01 
4.1296 00 -2-037E 01 


■o:o 

l. 3726-04 
1.077C-03 
' 3. 4576-03 
7.938G-03 
1.4026-02 
2.4406-02 
3.7146-02 
5.2986-02 
7. 2C 9 6-02 
9.4516-02 
1,2026-01 
1.4926-01 


" 0.0 

1.391E-02 
5.2896-02 
1.1206-0 1 
1.H71E-01 
2 . 7466-0 1 
3.7146-01 
4 .7486-0 1 
5.8266-01 
6.9276-01 
0.0346-01 
9.1336-01 
1.021C 00 


1.C206 02 
1.C2C6 02 
I.C2CE C 2 
1.02C6 02 
1.0206 02 
1.02C6 02 
1.C2CE C 2 
1.0206 02 
1.C206 02 
1.C2C6 02 
1.C2C6 C 2 
1 . C 2 C 6 C 2 
1 .0206 02 
1.0206 02 
1.0206 02 


1.8 70 6 C 2 
1 .870 6 C 2 
1.0706 C 2 
1 . 8 706 C 2 
1 .8 70E C 2 
1.8706 02 
1.8 70 6 C 2 
1 . 8 TOE C 2 
1.8701: C 2 
1 
1 
l 
l 


2-1786 
1 .9206 
1.6546 
I . 22 26 
1 . 1 C 6 6 


CC 

OC 

CC 

OC 

OG 


8.2666-01 
'5.5C66-C1 
2.7386-01 

-2 . 75 66-C4 

8 7 06 C 2 - 2 . 7 C 2 6 - C 1 
8 706 C 2 -5. 2496-01 
8 ?C 6 C 2 -7.9296-01 
8706 C 2 -1.C436 OC 


1.0706 02 -1.2656 00 
1 .8706 C2 -1.5176 OC 


- 2 r-’ 3 5 6 

Cl 

3 . 3 3 6 E 

00 -2.2966 

01 

1.8146-01 

1.1266 

00 

- ? . 9 o C: 

01 

2. 704 E 

CO -2.4396 

01 

2. 1676-01 

t. 2266 

00 ' 

- 2 . ‘M 1 G 

01 

2-7226 

00 - 2*. 6 5 2 6 

01 

2.5496-01 

1.3216 

00 

- 2 . n 5 7 6 

01 

3-09 3 6 

00 -2.987E 

01 

2.9596-01 

1.4116 

00 *■ 

- 3 - ' 6 5 E 

01 

3.C65E 

00 -3. 3 146 

01 

3.3546-01 

1.4946 

00 

- 3 . *> 5 3 6 

01 

2.3976 

00 -3.4706 

Cl 

3.0546-01 

1.5706 

00 “ 

- 3 . r 2 2 E 

01 

1.7086 

00 -3.5316 

01 

4 . 339C-01 

1.6396 

00 

-3. '-686 

01 

1.4546 

00 -3.5806 

01 

4.0366-01 

1 .701E 

00 

- 3 . ■» 9 3 E 

01 

1.3026 

00 -3.6446 

01 

5.3556-01 

1. 7546 

00 

-3- n 95E 

01 

1.2396 

00 -3.7336 

0.1 

5.0086-01 

l.CCOE 

00 

- 3 . n 7 4 6 

01 

1 . 1 9 4 E 

00 -3,8446 

01 

6.4346-01 

1-0376 

00 

-4 . rt 29E 

01 

, 1.1046 

00 -3.9646 

01 

6. 9896-01 

1.066E 

00 ~ 

-4.^426 

01 

9.2346' 

-01 -4 .0736 

01 

7. 552 £-01 

1.BC76 

00 

- 4 . C 7 l e 

01 

6.50 6 6' 

-01 -4. 1516 

01 

B. 1216-01 

1.099E 

00 

-4.05 8 6 

0 l 

3.240 6: 

-01 -4.1906 

01 

0.6 916-01 

1.9036 

00 

» 

-4.0226 

01 

- 1 . 74 9 E 

-02 -4 . 192E 

Cl 

9.2626-01 

1.8996 

00 

- 3 . 9 7, 5 U 

01 

-3.5816 

-01 -4.1586 

01 

9.0306-01 

i.ee7E 

00 

-3.8876 

01 

-7.0526 

-01 -4.0856 

Ql 

1.0396 00 

1.6676 

00 ~ 

-3.7836 

01 

-1 .0786 

CO -3.9636 

01 

1.0956 00 

1.0406 

00 

- 3 . 7.706 

0 1 

- 1.4856 

00 -3.779E 

01 

1.1506 00 

1.0056 

00 

-3.5346 

01 

-1.0686 

00 -3.536E 

01 

1.2C36 00 

1.76 3 6 

00 


v*> 

W 



2.25E6-C1 
3.49n-ci 

3 . 4 » E “ c l 
3.t*#ee-ci 

3.7586-01 

3. E50E-C1 
3.956E-C1 

4 . C 9 7 £ - C 1 
^-1S7S~C1 

4. Z57E-CI 
4.3576-01 
4.457E-C 1 
4.557E-CI 
4 . £ 5 7 E - C 1 
4.7576-01 
4.657E-CI 
4.557 E-C 1 
5.C5 7 E-Cl 
S. 1*5 7 E-C 1 
5.2576-Cl 
5.357E-C1 

' 5.4S7E-CI 

5 . 557 e-c i 

5. £576-Cl 
5.7? 7S-CI 
5 , 8 9 7 E - C l 

5 . 5 5 7 E - C L 

6. CS7E-C1 
6 . 1? 7 E-C 1 
£.2 5 7 E-C 1 

6.257E-CI 

6 - A 5 t S - C 1 
fc.5?te-ci 

•• t.tCtE-Cl" 

: 6.756G-C1 
; £ . E 5 6 E - c 1 
£.556 E-C l 
- 7.C966-C1 

7. I q 6 E-C 1 

7.256 E-C 1 

7.396E-CI 
7 . A 5 ft £ - C 1 
?•» 5 9 £ E - C l 

7. £5£E-C1 
7 . 7 v 6 E — C 1 

— T.e^te-ci 
7.55£E“C1 
6 .C56E-C l 
e. lCdE-C 1 
B.<9£E-Cl 

8. 24t6-C 1 

' e.45£c~Cl 

6. S96E-C 1 
8.656E-CI 

B.795E-C1 
8. 6956*01 
8.5S5E-CI 
S.C456-C1 
q. 1955-C1 
“ 9.2S5E-CI 

9 


2.5C9E 
3.6P3E 
3 . 9 A 3 E 
3. 2406 
1 .559E 

B. C82E 
1 . G £ ? C 
2.E65E 
3 . 6 C '3 £ 
2.35 6E 

C. 7C7E 
1 . 3CCE 
3 . 3 3 R E 
3 .3 m 
1.6ACE 

‘ 9. 157E 
2.556E 

3 . 5 £ 9 E 
2.341E 
9 . A £ C E 

1 .5 Ate 
- 3 . 2 £ 7 E 

2 . 19 IE 
i . e a c e 

9 . 5 2 A £ 

1 . 7 A 1 E 
3.1566 

2 , 3 C ( ; £ 

2.53 7 6 
1.552E 
1 .05 8 E 
1.422E 
2-3456 
3.C83E 
3. 2 1 7E 
2. 1576 
2 - 9 7 A E 
2 . £ £ AS 
A1AE 
2 2 56 
119 6 
06 4 6 
C £ 3 E 
1 126 
2 . 2 C A E 
2.3326 
2 . A 7 9 E 
2 . £ 1 8 E 
2.7116 
2.716E 
2 . £ C 36 
2 . 3 5 9 E 
2 . C 1 A E 
l . £ 6 7 E 
1.9526 
1*5676 
2.06 3 E 
2.6316 
2.66CE 
2.6AIE 
2.1816 


Cl 
Cl 
01 
C 1 
Cl 

cc 

C 1 
Cl 
C 1 
C 1 

cc 

Cl 
C 1 
C l 
C 1 
CC ' 

01 
C l 
C 1 
CO 
Cl 
Cl' 

Cl 
C l 
cc 
01 
C l 
C 1 
C 1 
Cl 
C 1 
C 1 
C 1 

Cl“ 
C 1 
01 
Cl 
Cl 
C 1 
Cl ' 
Cl 
Cl 
Cl 
C 1 
Cl 
Cl 
C 1 
C 1 
01 
01 
C l 
Cl 
01 
Cl 
Cl 
01 
01 
Cl 
Cl 
Cl 
Cl 


1.39 PC 
l . L A 7 E 
1 .02 1 E 
1.C926 
l . 2 7 e E 
1 . A 7 0 C 
1 . A 2 E C 
C 0 3 E 
Ct 3 E 
0966 
39 2 E 
, 3 7 £ E 
1 . 0 A C £ 

1 . 0 2 7 6 
1 . 32 7 C: 

1 .APfcE 
1. L97E 
9.R2CE 
1 . U: A E 
1 . A 1 ft E 
1 . 32 7 E 
"1.027E 
1 .Cl 1 E 
1 . 2 A 8 E 
1 . A 3 (' E 
1 . 3 ? 9 C 
1 . 0 7 3 E 
' 1.CC96 
1. 1296 
1. 322E 
1 . A 2 A E 
1 . 39 2 £ 
l. 22 IE 
“ 1 .06CE 
1 . 0 1 2 E 
1.0166 
C 5 4 E 
1 C A 6 
15 IE 
1 9 C E 
1.217E 
1 .23 IE 
23 1 E 
2 1 7 E 
1 9 C E 
1526 
l C 5 E 


1 . Q 5 6 E 
1 .01 AE 
9 , 0 9 9 E 

9 . 9 1 5 E 
1 .026-6 

1 .09 IE 
l . 17 A 6 
1.2A2E 
1.2616 
1.21CE 
1 .1356 
1.1C66 
1 . 1 A 9 £ 
1 .24 1 e 


C 2 
C 2 
C 2 
C 2 
C 2 
C? 

C 2 
C 2 
C 1 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
01 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C2 
C 2 * 

C 2 
C 2 
C 2 
C 2 
C 2 
C2 ” 
C 2 
C 2 ' 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 2 
C 1 
C l 
C 2 
02 
C 2 
CH 
C 2 
C2 
C 2 
C 2 
C 2 
C 2 


1 . C 2 0 E 0 2 
1 . 0 2 C i: 0 2 
1 . C 7 0 E 0 2 
1 . C 2 C 6 C 2 
1.C2CG 0? 
1.0296 02 
1 . C 20 C 02 
1 . C 2 C 6 0 2 

i . c 2 a r. 0 2 

1 . C 2 C E 0 2 
1.C2CG 02 
1.C20E 02 

1 . c ? c E 0 2 
1 .0206 02 
1.C2CE 02 
' I.C2CE C 2 

UC20E 0? 
1.0206 02 
1 . C 2 0 I: C 2 
1.C2CE C 2 
1.C2CE C 2 
- 1 . G 2 0 6 0 2 
1 . C 2 0 E 0 2 
■ 1.C2C5 C 2 
1 . C 2 C 0 0 2 
1.C2CG 02 
1.020 6 C 2 
• 1.C2CC 0 2 
1.C20C 02 
1 .0206 C 2 
I.C2C6 C 2 
1.C206 C7 
1.C2CE 02 
1.C2CG 02 
1.C2CC- C 2 
1.C2CE 02 
1.020c 02 
1 . 0 ? C E c 2 
1.C2CE C 2 
‘ 1 . C 2 C E C2 
1..C2CE 02 
1.0206 02 
1.C20E C 2 
1.C2CE 02 
1.C20E C 2 
1.C2CE C 2 
1.C2C6 C 2 
1.02CE 02 
1.020E 02 
1.0206 02 
1.C2CE C 2 
1.C206 02 
1.02CE 02 
1.C20E C 2 
1.C20E 02 
1.02CE 02 
1 « C 2 0 E 02 
" 1.C2CE C 2 
1*0206 C2 
1*0206 02 
l . 020E 02 


1.0706 C 2 
1.(37 0 6 C 2 
1.9706 C? 


t.n/or c 2 
1 . P 7 0 1: C 2 
1 . P 7 0 £ C 2 
1 . 11 7 C E 0 2 - 
1.R7CE C 2 - 
1.9 7 0 E C 2 . - 
1.8706 02 - 
1.8706 C 2 - 

t.nvot c 2 - 

1 .8706 C 2 - 
l . 8 7 C f: C 2 - 
1.8 TOE C 2 * 
1.8 7 OF: C 2 • 
1.R70E C 2 • 
1.8 70C C 2 ■ 
1.R706 C 2 
1.87CE C 2 
1.8706 C 2 
’ l . 8 7 0 E C 2 
1 . 8 7 0 E C 2 
1.9706 C2 
1 .8 fOE C 2 
1.8 7CE C 2 
1. ft 706 C 2 
1.8 706 C 2 
1.C706 C 2 
1.8706 C 2 
1.8,7 0 6 C 2 
1 .9706 C 2 
1 .8706 C 2 
1.8706 C 2 
1.9706 C 2 
1.8706 C2 
1.8706 02 
1.870E C 2 
* 1 .8706 
1.8 70 6 C 2 • 
1 .8706 C2 
1.8 70E C2 
1.8706 C 2 
1.C7CE C 2 
1.C70E C 2 
1.870E C 2 
1.8 706 C 2 
1.C706 C 2 

1.8706 C 2 
1 . 8 70 E C 2 
1.8 706 C 2 
1.870 6 C 2 

1.8706 02 
1.870E C2 
1.8 706 C 2 

1.8706 C2 
1.0706 02 
1.8706 C 2 
1.0706 C 2 
1.870E ,C2 
1.8706 G2 


00 

00 

CC 

CC 

c c 

CO 

cc 

cc 


1.7:91 
1 . 9 A (1 1 
? . 1 A 6 1 

2 . :s 3 1 ( 

2 ,5C 3C 
? . £ £ 2 1: 

? . 0 C 6 f- 
■ 2 . 9 2 5 T 
• 3 . c e . o r 
•3.1501: 
-3.2256 

- 2 . 2 C A 6 
-3.3596 
-3.3586 

- 3 . A 2 2 £ 

-3.4 2 IE 
-3.4266 

- 3 . A C 7 £ 
-2.3736 
-3.3266 
-2.2(66 
-3. 16 A £ 

- 3 . 1 C 9 6 
-3-0136 

- 2 . 5 C 6 6 

-2.7696 
-2.6£3£ 
-2.5286 
-2.36AE 
-2.2 2 A 6 

- 2 . C 7 7 6 
-1.51A6 
- 1 . 7 A 6 6 

- 1 . 5 7 A £ 
-1.359c 

1.2216 
1 .04 It 
8 . £ C <5 6 - 0 1 
6. EC IE “Cl 
A . 5 5 R 6 - 0 1 
3.2C86-01 
1 . A 3 8 £ - C l 
3.C67E-02 
2 .C 1 06-0 1 
£851 


3.9806 01 
3.2096 01 
3.02 31 
2.8 2 31 


0 1 
Cl 


-2.0096 01 


01 - 


oo - 
cc - 
cc * 
cc - 
CO ■ 

co ■ 
oc 
oc 
cc 
cc 
cc 
cc 
oc 
cc 
00 
cc 
cc 
cc 
oc 
cc 

GC 

cc 

oc 

cc 

oc 

oc 

cc 

CO 


2 * 


-0 1 
3C9C-C 1 

et- os-0 1 
,3546-01 
,8446-0 1 


1.1236 

1 . 2 5 A £ 

1.3776 

1 . A 5 2 6 

1.5696 

1 . 6 6 B 6 

l . 7 E 7 £ 

1.0696 

1.6396 

2 .CC 16 

2.C536 

2.0566 


?. 3 ft A i: 

2 . i A 1 1 6 Cl “ 
1.9036 01 ■ 

1.6506 01 ■ 
1.390E 01 ■ 
1.1256 Cl 
■8.5596 CO 
•5.8 8 81: 00 
-3.1016 00 
-3. 6 226-01 
2-3686 00 
5-0766 00 
7.7516 oo 

1. C3PE 01 
1.296E 01 
1 . 6 4 7 E 0 1 
1.7906 01 

2. C25E 01 
7.2506: 01 
2 . A 6 8 6 0 1 
2.6696 01 
2.8626 Cl 
3.0426 01 
3.2096 01 
3.7636 01 
3.5C76 01 - 
3.6296 01 - 
3.7396 01 - 
3. 93 5 6 01 - 

3.9166 Cl - 
3. 9 (126 01 - 
A. 0326 01 - 
4. COPE 01 ' 
4.0886 01 • 
4.0946 01 ‘ 
A . 0 R 5 E 01 ' 
4.0616 01 ' 
4.C24E 01 
3.9726 01 
3.908E 01 
3.83CE 01 
3 . 74 06 0 1 
3.6796 01 
3.5266 01 
3.403E Cl 
3.2706 01 
3.1286 01 
2.9776 01 
2.8196 01 
2.6536 01 
2.4026 01 
2.3056 01 
2.12AE Cl 
1.9396 01 
1.7516 01 
1.5616 01 


CO 

00 

oc 

CO 

oa 

cc 

co 

oo 

oc 

cc 

oc 

cc 


? . C 7 3 E 00 
7.0776 00 

1.89 06 00 
1 .8 !I7 6 CO 

2.0 j 4 n co 

2.5076 00 
3.0406 CO 
3.1836 00 

• 2.904 6 0 0 
■2.8076 00 
•3.1456 00 

• 3 . (; 4 8 6 00 
-3.5716 00 
-3.1406 00 
-3.1556 00 
-3.563E CO 
-3.7096 00 
-3.2926 00 
-3 . C*i 3 E CO 
-3.2466 CO 
-3.5566 CO 
-3.3346 00 
-2.6B7E 00 
-2.755E CO 
-2.9476 00 

-3. 0966 00 
-2.8196 00 
-2.3726 CO 
-2.1036 00 
-2.0 A 86 00 
-2.1146 00 
-2.1426 CO 
- 1.9ft 5E 00 
- 1 .6606 00 
- L . 3 2 8 6 CO 
1.0476 CO 
8.2866-01 
6.5426-01 
5.0346-01 
■3.6106-01 
•2.2076-01 
•7.6CCE-02 
7.3146-02 
2.75GE-0 1 
3.7596-01 
5. 2136-Cl 
6.5716-01 
7 .e04E-01 
8.9216-01 
9.9736-01 

1.1C5E 00 
1.2286 CO 
1.3706 00 
1,5646 00 
1.77 EE 00 
1.9786 CO 
2.C93E 00 
2.0E7E 00 
2.0136 00 
1.9636 00 
1.9906 00 


- 3 . 2 7 7 E Cl 
-3.0686 01 

- 2 . 9 1 9 C 01 
-2.7536 Cl 
-2.6416 01 
-2.4156 01 
-2.0596 01 
-1.0C16 Cl 
-1.6146 01 
-1.4566 01 
-1.2C9E 01 
-8.4926 00 
-5.3656 CC 
-3.3826 00 
-1.1.4 16 CC 

2 . 2 5 2 E 0 0 
5.0696 CC 
P.129C CO 
9.924 E 00 
1.2 A 96 Cl 
1.5866 01 
1.8506 01 
2 . 0 C 9 6 Cl 
2.1766 01 
2. A 176 0 1 
2.7056 01 
2.9256 01 
3.0566 01 
3.1416 01 
3.2896 01 
3.4566 01 
3.6456 Cl 

3.0096 01 
3.915E Cl 
3.9726 Cl 
3.9596 01 

4.0096 01 
4 .0096 01 
4.0C26 01 
3.9876 01' 
3.9636 01 
3.9306 01 
3.8856 01 
3.C3C6 01 
3.7646 01 
3.6806 Cl 
3.6066 Cl 
3 • 5 1 9 E 01 
3.4306 OL 
3.3406 01 
3.2476 01 
3.1476 01 
3.0326 01 
2.395E 01 
2.727E 01 
2.5306 01 
2.321E 01 
2.1246 Cl 
1.949E 01 
1.7846 01 
1.6116 01 


1.2556 00 
1.3066 00 

1.354 I- 00 

1 . 4 C 1 6 00 
1 . 4 4 A 6 00 
1.4896 00 
1.5296 CO 
1.5666 00 
1.6C1H 00 
1.63 A 6 00 
1.663E CO 
1.609E 00 
1.7136 00 
1.733c 00 
1.7516 CO 
1 . 7656 00 
1.7776 CO 
1.785c 00 
1.79C6 00 
1.7626 CO 
1.7916 CO 
1.707c CO 
1.70C6 00 
1.7716 CO 
1.7596 00 
1.7446 00 
1.7276 OC 
1.7C7E 00 
1.6656 00 
1.6616 00 
1.634E CO 
1.606E 00 
1.577c 00 
1.5 A 56 00 
1.5126 00 
1.4706 00 
1.4436 00 
1.4C7E 00 
1.3706 00 
1.3336 00 
1.2956 00 
1.2576 CO 
1.2196 00 
1.1896 00 
1.1426 00 
1.1C4E 00 
1.0676 00 
1.030E 00 
9.9446-01 
9.55 AE-C l 
9.2546-01 
8.9266-01 
8 .6096-01 
E.3C76-CI 
8.01 BE -0 1 
7. 7456-01 
7.488E-C1 
“ 7.2476-01 
7.0256-01 
6.3206-01 
6.6336-01. 


1„71?*E 00 
1.6506 00 
1.5976 00 
1.5 3011 00 
1.4576 00 
1 . } 0 0 E 00 
1.2 9 RE 00 
1.2UE 00 ’ 
1.1226 00 
1.0296 00 
9.3276-01 
8 . 346E-0 1 
7.3466-01 
6.333E-01 
5. 3C9E-01 
4.2eiE~01 ' 
3.2526-01 
2.2266-01 
1.2C9E-01 
2.C37E-C2 
-7.8566-02 
-U 7 5 5 6-01 
-2.7C1E-C1 
-3.6196-01 
-4.5036-0 l 
-5. 3626-01 
-6 . 1 3 16-01 
-6.9606-01 •• 
-7.6976-01 
-0.3906-01 
-9.C37E-0L 
-9.636E-01 
-l.OlftE 00 
- 1 .C68E CO ' 
-1.1136 CO 
-1.1526 00 
-1.1366 00 
- 1 . 215c 00 ' 
-1.2386 00 
-1.2556 00 ' 
-1.2686 00 
-1.275E 00 
-1.2766 00 
-1.2736 00 
-1.2646 CO 
-1.251E 00 
-1.2326 00 
-1 . 209 E 00 
-1.1826 00 
-1.15C6 CO 
-1.1156 00 
-1.0756 00 
-1.032E 00 
-9.857E-C1 
-9.362E-01 
-8.8396-0,1 
-8.2916-01 
-7.720E-01 
-7.1296-01 
-6.5206-01 
-5. 8986-01 


9. A55E-CI 

1 .760E 

Cl 

1.33CE 

C 2 

q .55 56- Cl 

1 .7 ICG 

Cl 

1 . 357E 

C 2 

q.tbse-ci 

2.C5PE 

Cl 

1.7b a 6 

02 

5. 7956-0 1 

2.6C9E 

01 

i . ibce 

C? 

9.eV;f-Cl 

2.CC 7E 

C 1 

l . 1 2 C F 

C 2 

O.S95E-C1 

2. 58 IE 

C 1 

l - 1 A 1 t 

C2 

1.CC5E CO 

2 . 1 2 5 E 

C l 

1 . 7 1( E 

C 2 

1.C15E CO 

1.7 2 2E 

Cl 

1 . 1 8. 5 f: 

C 2 

I.C29E CC 

1 . 1 1 A E 

Cl 

1..3CAC 

C 2 

I.C2AS CO 

1 . b?6E 

C l 

1 .25 IE 

C 2 

l.C'-SE CO 

2 . A A 3 E 

0 1 

1 . 1 ft A E 

C 2 

1 . C 5 A 6 CO 

2 . 6 b A h 

C 1 

1 . ICE 

C 7 

l.CtAE CC 

2 . 6 ? 5 E 

Cl 

1 . iCbE 

C 2 

1 . C 7 5 E CC 

2 . 2 3 6E 

C 1 

1.15 76 

C 2 

i . c 2 b e co 

i . b b * E 

Cl 

1 . 2 2 3 6 

C 2 

l . C b A £ CO 

1 . 7 5 7 E 

C 1 

1 . 2 7 b E 

C? 

1 . lCbE CC 

1 .71 AE 

01 

1 . 3C 16 

C 2 

1.1 l 3 l CC 

l .E 766 

C l 

1-2016 

C 2 

1.128c CO 

2 . 1 5 A E 

Cl 

1 . 2 3 7 6 

C 2 

1.1386 CO 

2 . A 2 U 

Cl 

i. iet£ 

C 2 

l.lAtt CO 

2 . 5 9 A 6 

Cl 

1 . 1 A B F: 

C 2 


1.15SE CC 2. £585 Cl 


1 . 1 * S E 

CC' 

2 . t 39E 

01 

1. 1 7 e e 

CO 

2.572E 

Cl 

i . \ ?et 

CC 

2 . A o 7 6 

C 1 

1 . 1 988 

CC 

2 . A C 2 E 

C 1 

1 . ?C2£ 

CO 

2 - 3 2 BE 

C l 

1 . 2 1 8 E 

CC 

2.27 IE 

C 1 

1.2286 

CC 

2.23 1C 

Cl 

l . 2 ; £ c 

CO 

2 . 2 C B F. 

Cl 

1.2-76 

CO 

2. 2C 1 1 

C 1 

1 . 2 5 7 f 

CC 

2.2C0E 

Cl 

l . 2 1 t £ 

CC 

2.21 (S 

C 1 

1-2776 

CC 

2.2551 

C 1 

t ,?P7E 

CC 

2 . 2 b C E 

Cl 

1 .2b7E 

CC 

2.3251 

Cl 

1 . 2 C / E 

CC 

2 - 3 ft 5 C 

C 1 

1.3 1 7C 

CO 

2 . 3 b A 6 

01 

1.3278 

CO 

2 . A C 6 E 

C 1 

1 . 33 1 E 

CC 

' 2.2btE 

Cl 

1 . 3 A 7 E 

CO 

2.357C 

0 l 

1 . 3 5 7 E 

CO 

2.250c 

01 

1.3C6E 

CC 

2 . lbbE 

Cl 

1.3 7 1 6 

CC 

2 . C b 7 1 

C 1 

1 . 3F6E 

CO 

2 .CC 3 E 

Cl 

l . 3 s ft E 

CO 

l . b A 3 E 

0 1 

1 . ACfcE 

CC 

i - s i n e 

Cl 

1 . A 1 6 E 

CC 

2.CCCE 

Cl 

1 . A 2 t E 

CC 

2.1176 

Cl 

l . A 3 d E 

CC 

2.253c 

C L 

1 . A A t E 

CC 

2 . 3 6 5 E 

Cl 

l.AStE 

CO 

2. A2 IE 

01 

1 . A c 6 E 

CC 

2 . A C 7 E 

Cl 


1 2 *6 

1 2 2 £ 

I A 6 E 

ltzr. 

II If: 

1 ^ c e 

195 6 

2 C '3 r 
zcve 

1 . 2 c ‘3 C 

l./CCE 
S.l'/lC 
lf'CE 
It 7E 
lb 1C 
IMC 
me 
1 2 7 E 
i . i ? b r 
1.120C 
1 . 1 5 5 F 
1.177E 
1.2C 1£ 
1-222C 


.A7dE 
. a 5 f 

. . A <; 5 E 
.. 5CbE 
L.515E 
L .52 5 6’ 
1.5356 
1.5A5E 


CO 

CC 
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- A . A f: 8 E - 0 1 
-A.628C-C1 
-A .76 7 6-0 1 
-A .8366-01 
-A. 9356-01 
-5-C63E-C1 

1 196-01 
1 5 A 6 - 0 1 
-5.1ftt6-01 


-5 

— 9 


A . 5 A I! 6 
A . 7 1 7 E 
A . P ft A 6 
A . 9 3 7 6 
5.0886 

5 .u(- r 

5.7726 
5 . 7 5 A C 
5 . ? ft 5 E 
5 . 2 5 A F 
5.7726 
5 . 1 <956 
5.09 5 6 
S.C07E 
A . 11 9 1 6 
A . 7 ft 1 F 


CO 

00 

QO 

00 

CO 

00 

CO 

no 

CO 

00 

oc 

00 

00 

00 

CO 

00. 

00 

CO 

00 

CO 

oc 

CO 

cc 

oc 

00 


-5.155E-C1 
-5. 1216-01 
-5-0656-01 
-A. 985 E- 01 
-A.S8AG-01 
-A.7ME-01 
-A. 6156-01 
-A.A5CE-01 
-A .28C6-U1 
-A * 0 6 ft C - 0 l 
-3.8796-01 
-3.67,06-01 
-3. A 1CL-01 
-3. I'UL-OV 

- 2 . 9 A 5 6 - 0 1 
-2 . 6976-0 l 
-2. A3 A 6-01 
-2.1726-01 
-1.9076-01 
-1 .6 ACE-01 
-l . 372F-C1 

, -1 . 1056-01 
-0.30 AC-02 
-5. 73 96-02 
-3. 1206-02 

- 5 . A 9 7 6-03 
1 .9806-02 
A . A 5 5 6 - 0 2 
6.86 1 E-C? 
9.19 AC-02 
1 . 1AAG-01 


7.8386 00 
7.9576 00 
■ 8 . 0 A l 6 00 
.8.0916 CO 
•O.IOUC 00 
-8.0976 CO 
■0.0A7C 00 

- 7.9ft 16 00 
-7.0A96 00 
-7.7C6E OC 
-7.53AE CO 
“7.335E OC 
-7.1096 CO 
-6.358c CC 

- 6 . 5 8 A 6 00 
-6.2C8E 00 
-5.972E 00 
-5.63VE 00 
-5.788E CO 
-A. 9226 00 
-A.5A3E* 00 
-A. 1526 CO 

— 3 * 7 A 9 6 00 
-3.338E 00 
-2*9196 00 
-2.A536 CO 
-2-063E 00 

— 1 . 7> 2 9 E 00 
-1.19AE 00 
-7. 556E-0 l 
-3. 2 7 66-01 

9.9896-02 
5.2066-01 
9.326E-0 l 
1.33AG 00 
1 . 7 2 2 C 00 
2.0956 CO 
2.A52G CO 
2.7916 CO 
3. 1 HE 00 
3 . A 11 6 0 0 
3 . ft 9 0 E 0 0 
3 . 9 A 7 6 OC 
A . L 0 2 F 00 
A.35AF 00 
A. 5 GAG 00 
A.750G 00 
A . 0 9 A 6 00 
5-01A6 CO 
5 . I 1 2E CC 
5.1876 CO 
5.2396 00 
5 . 2 ft9 6 00 
5.2776 00 
5 . 2 0 A 6 00 
5.2296 00 
5 . 1 7 5 C CO 
5. 1 OOE 00 
5.0C7E CC 
A . 85 A 6 CO 


9.0716-01 

2. 703E-01 

9. 153G-01 

2.777E-01 

9. 237G-01 

2.830E-01 

9. 3736-01 

2.0066-01 

9. A 10 6-01 

2.92 1 E -01 

9.A58E-C1 

?. 9A3E-0 1 

9 . 5 8 ft E - 0 1 

2.9526-01 

9. 675 E -01 

2.9A9E-01 

9.7630-01 

2-93 A £-01 

9 . 0 5 1 E - C 1 

2.9066-01 

9. 93 7 £-01 

2. 668E-01 

1.0026 00 

2. 818E-01 

1.01 IE 0 0 

2. 797E-01 

1.0196 00 

2.686E-01 

1.027E 00 

2.6056-01 

1.0 3 AS 00 

2.515E-01 

1.0A2E 00 

2. A 166-01 

1.0A9E 00 

2. 3C56-01 

1.055G 00 

2. 195E-01 

1.0626 CO 

2.0736-01 

l.OfcCi: 00 

1 . 9 A ft 6 - 0 1 

1.0736 CO 

1 ,e 136-01 

1.0706 00 

1 .67AE-01 

1.083b 00 

1.532E-01 

l.OCHE CO 

U 3866-01 

1.0916 CO 

1 . 2376-01 

1.0956 00 

1 .00ft£-01 


1.0986 00 
1 . 1C0E 00 
1 . 1C 7 6 CO 
1.1CAE 00 
1 . 1056 00 
l . 1C 6 E CO 
l.lCftE 00 
1 


9. 33AE-02 
7.7586-02 
6.2596-02 
A. 725E-02 
3*2016-02 
l .6936-02 
2.C7AE-G3 


106E 00 -1. 2506-02 
1. 1056 00 -2. 67 A £-02 
1. KA6 CC -A.C59E-02 
1.1036 00 -5. A COG*- 02 
1.1C1F CO -A.A92E-02 
1.0996 00 -7.932E-02 
1.0506 00 -9.11AE-02 
1.0936 00 -1.C2AC-01 
1.C9C6 00 -l* 129G-01 
1.0806 00 - 1 . 2286-01 
1.082E 00 -1,3206-01 
1.0786 00 -1.A05E-01 
1 . 0 7 A E 00 -1* A 826 -01 
1.0696 00 -1.5516-01 
1.0656 CO -1 .6136-01 
1 . 0 ft 0 6 GO -l. ft 6 76-01 
l , C 6 A b 00 -1.7 L 36-01 
1 . 0 A 9 6 00 -1.7506-01 


l • C A A 6 00 -L.7PC6-01 
1.0386 00 -1.8O2E-01 
1.033E 00 -1.8166-91 
1.0276 00 -1. 8226-01 
1.0226 00 -I.S2CC-01 
l.OlftE CO -L.811E-01 
1.C11E 00 -1.79AG-01 
l.OCftG 00 -1-77CE-01 
A.765E 00 1.0006 00 -1.739E-01 




1.2C3E C2 
I.2C2E C2 
|.2C'iE C2 
1.2C3E C2 
1.2C2E C2 
1.2C96 C 2 
1-2C9E C2 
1 . 2096 C2 
1.2C96 C2 
'1.2C5E C2 " 
1.2056 C 2 
1.2056 G 2 
1.2C6E C2 
1.2C66 C 2 
1.2C76 C 2 
"1.2C7E C 2 
1.2C8G C 2 
1.2086 C 2 
1.2C56 C 2 
1.2C5E 02 
1.2C56 C 2 
— 1 .2 ICG C 2 
1 . 2 1 C 6 C 2 
1 . 2 1 C 6 C 2 
1.21CE C 2 
1.2116 C 2 
1.2UE C 2 
-1.211C C 2 ~ 

1 . 2 1 l E C 2 
1.211E C 2 
I.2UE C2 
1.211E C2 
1.2U6 C2 
1.2UE C 2 
1.2U6 C 2 
1.21CE C 2 
1.2 ICE C 2 
1.21C6 C 2 
U21CE C 2 
"1.21CE C 2 
1.21CC C 2 
' 1.21CE C 2 
1.2106 02 
■ 1.21CE C 2 
1.21C6 C 2 
- 1.21CE C2 
1.21CE C 2 
1.21CE C 2 
1 . 2 1 C E C 2 
1.2 ICE C 2 
1.21CE C? 

1 . 2 1 c r ; C 2 
1.2 ICE C 2 
1 . ? I C 6 0 2 
1 .2 1C6 C 2 
1 . 2 1 C E C 2 
1.21CE C 2 
' 1.2106 C2 
1.2 ICC 02 

1 .2 ice C 2 

t.ZICt 02 


1 .0706 02 
1.0206 C 2 
I.C206 02 
1.C20E C2 
1 .0206 C 2 
1 . C 2 C E C 2 
UC20C 02 
1.0206 02 
1.C2CE C 2 
1 .0 20 6 02 
1.0206 02 
1.02CE 02 
1.C2CE 02 
1.C2CE 02 
1.0206 02 
1.0206 02 "' 
1.C2C6 02 
1.0206 02 
1 .0 206 02 
1.0206 02 
1.C2CE 02 

- 1.C2CE 0 2 
1 .0206 02 
1.0206 02 
1.C2CE C 2 
1.C2C6 02 
1.C206 02 

- 1.0206 02 
1.C2C6 C 2 
1.C2C6 02 
1 .0206 02 
1 . C 2 0 6 0 2 
1 .0206 02 

" 1 . G ? 0 E 0 2 " 

1 . C 2 0 6 0 2 
1 .0206 02 
1.C2C6 02 
1 . C ? C 6 0>2 
1.0206 02* 
-1.C20E C 2 " 
1.C206 02 
1.0206 02 
1.0206 02 
1.C20E 0 2 
1.C20E 02 
— 1 .0206 02 
1.C2CP 02 
1.C2CE 02 
1.0206 02 
1.CZ0E 02 
1.C206 02 
1.C2CE 02 
1.02CE 02 
1.0206. 02 
1.0206 02 
1 .0 206 02 
1 . 0 2 C 6 0 2 
1 . 0 2 0 C 0 2 
1.0206 02 
1 . C 2 0 6 0 2 
1.0206 02 


1.0706 C2 
1.0706 C? 

1.0 70 6 C 2 
1.0706 02 
1.070E 02 
1.0706 C 2 
1.870E C 2 
1.670c 02 
1.0 70 6 C2 _ 
1.0706 02” 
1.0706 C2 
1.0706 C 2 
1.8706 C2 
1.870E 02 
1.8706 02 
’ 1.8706 C2 
1.0706 C 2 
1.87CE 02 
1.870E C2 
1.8706 C2 
1.6706 C 2 
‘ 1.0706 C 2 
1.0706 02 
1.6706 C 2 
1.0706 C2 
1.P70E C 2 

1.0706: C 2 

- 1.B70E C 2 
1.0706 C 2 
1.0706 C2 
1.C70C 02 
1.0706 C 2 
1.0706 C2 

1.8 70 6 02 " 

1 .0706 C 2 
1.8 7C6 C 2 
1.0706 C 2 
1.8706 C 2 
1.8706 02 ' 
"1.0706 C2 
1.0706 C 2 
1.0 70 6 C 2 
1.P70E C2 
1.8706 02 
1.8706 C2 
•'1.6706 C 2 
1.0706 02 
1.8706 C 2 
1.0 706 C 2 
1.0706 C 2 
1.8706 C2 
1.0706 02 
1.0706 C 2 
1.0706 02 
1.0706 02 
1.0 7 06 C 2 
1.11/06 C 2 
l . 0 / 0 C 0 2 

1.0701. C 2 

1 . 0 7 0 C C 2 

1.0 /06 0 ? 


3 . 

93 16 

00 

-i 

1936 

00 

2. 

<59 7 6 

CO 

2. 

(996 

00 

2. 

9356 

00 

2 . 

17 16 

00 

1. 

9096 

00 

1. 

.6396 

00 

1 . 

.3(26 

00 

1 - 

.09 16 

oo 

8. 

.2016- 

-01 

5 . 

.5166- 

-01 

2 

, 86 3E 1 

-01 

2 

. 5076 

-02 

- 2 . 

.nee 

-01 

-9 

."096 

-01 

-7 

.•>376 

-01 

-9 

. r 87E 

-01 

- 1 

. i 0.56 

00 

-1 

.'•0 2 6 

: oo 

- 1 

.'ORE 

: oo 


1.3676-01 9.6196 00 

1 . 5 t 9 6 - 0 1 A . 9 5 1 6 GO 

I.uu-oi A. zMt: oo 

1.9 A 06-01 A.O016 0 0 

2.1206-01 3-0766 00 

2.2616-01 3-6976 00 

2 - a c r > C - 0 1 
2.5(96-01 
2. 65 7 6 -Cl 
2.7966-01 

2-6516-01 
2.9726-01 
3.C276-C1 
3.C926-01 
3.1316-01 
3.1576-01 
3 .1696-0 1 
3.1(76-01' 

3.1536-01 
3.1266-01 
3.C66C-01 
3.0296-0 1 
2.9U6-C1 
2.6966-01 
2.8 11E-C1 

2-7166-01 
2.(126-01 

2. A 9 06-01 - 1 . ° 0 9 6 CO 
2 3776-C1 -l. r B8E CO 
2. 2 AOC-Ol - 2 . 1 6 0 6 00 
2.1126-01 -2.^206 00 
1.9706-01 - 2 . 6 7 6 CO 
1.6286-01 ~ 2 - * CC E 00 
1.(716-01 -2.-206 00 
1.5156-01 -2.0?76 00 
1 T=5r-01 -2.9196 00 
1*1536-01 -2.9976 00 
l . C 2 06 -0 V -3.0 5 16 00 
8. (A Of. -0 2 -3-1116 00 
6.9956-02 -3.1976 OC 
5. 3 A 56-02 -3. 170E 0 0 
3.7C6C-02 -3.1706 00 
2.0626-02 -3.1736 00 
A. 01U:-C3 -3.1556 00 
- 1. 05 16-02 -3.1236 00 
-2.(296-02 -3.0806 00 
-A. 127.6-0 2 -3.02 96 0 0 
-5.570C-C2 -2.9556 CO 
-6-97R6-C2 -2.8766 00 
-0.3106-02 -2.7096 00 
-5.5556-02 -2.69C6 00 
-1.CE0E.-Cl -2.50 IE CO 
-1.15 A 6-01 - 2 . 5 6 A 6 00 
-1 ,3C 16-0 1 - 2.3 396 00 
-1.9C06-01 -2.20(6 CO 

- I . A 5 1 6-0 1 -2.0 (.76 00 
-1.5196-01 -1-9226 00 
-l.( 996-01 -1.7/26 00 
-1.7151-01 " 1 .6176 CO 

- I 7 7 36-0 1 - V . 9 5 9 f 00 
- 1.2906 oo 


1.3606-01 
1.56*6-01 
l . 76 16-0 1 
1.9956-01 
2.1176-01 
2 . 2 7 8 6r 0 1 
2-92(6-01 
2.56 1 6-01 
2.6836-01 
'2.7926-01 
2 .08/16-01 
2.9706-01 
3.0386-01 
3.0936-01 
3.1396-01 
3.1626-01 
3.1766-01 
3.1766-01 
3 . 1696-01 
3 . -13 0 6- CL 
3. 10C6-01 
3.O5CE-01 
2.9076-01 
2.91AE-01 
2.0296-0 1 
■ 2.7396-01 
■ • 2.6306-0 1 
2.51(6-01 
2.3996-01 
2.2696-01 
2 .127 6-01 
l .9096-01 
1.8356-01 
1.601 6-01 
1.5296-01 
1.3636-01 
1.2006-01 
1 .0356-01 
8.6876-02 
7 .022E-O2 
5.3610-02 
3.7116-02 
2.0786-02 
A .67(6-03 
-1.1136-02 
-2.6506-02 
-9.1626-02 
-5.6196-02 
-7.02 9 E-02 
-0.3706-02 
-9.6516-02 
-1.0066-01 
-1.2016-01 
-1 .3076-01 
-1.9076-01 
-1 .99nE-0l 
-1.5016-01 
- 1 .65(6-0 1 
-1.7226-01 
-1 .7006-01 
-1.0301-01 


9.6196 00 
9 . 95(6 00 
9.2796 00 
9.0096 00 
3.0856 00 
3.6706 00 
3.9996 C C 
3.2096 00 
2.9(56 00 
2.7196 00 

2.9876 00 

2.1956 00 
1.9296 CO 
1.6606 00 
1.3896 OO 
1.1186 00 
8.9796-C1 
5.7516-01 
3.1316-01 
5 . 0886-02 
-2 . 065E-C 1 
-9.5816-01 
-7.0296-01 
-9.3596-01 
-1.1(06 GO 
-1.387E CO 
-1.5966 00 
-1.7996 00 
-1.9816 00 
-2.1556 CO 
-2. 3176 CO 
- 2 . 9 6(6 CQ 
-2.6016 00 
-2.7236 00 
-2.0306 CO 
-2.9296 CC 
-3.0036 CC 
-3.0606 00 
-3.119c 00 
-3.1556 CO 
-3.1706 00 
-3.1866 00 
- 3.1816 00 
-3.1636 00 
-3.1316 CO 
-3.0876 00 
-3.0316 00 
-2.9(96 CC 
-2.0856 00 
-2.7956 00 
-2.695E 00 
-2.50(6 00 
-2-9606 00 
-2.3926 00 
-2.2C9E GO 
-2.0696 00 
-1.9236 00 
-1.7726 00 
- 1.(166 00 
- 1.9576 00 
-1-2956 00 


9 .951*6-01 
9.9C16-01 
9.0526-01 
9.0096-01 
9.7596-01 
9.71 56-01- 
9.(736-01 
9.6396-01 - 
9.55(6-01 - 
9.5616-01 - 

9.5296-01 - 
9.9996-01 - 
9.9726-C1 - 
9.9966-01 - 
9.92(6-01 - 
9. 9086-01 * 
9.3526-01 ■ 
9.379E-01 - 
9.3(96-01 ■ 
9.361 E-Ql 
9.3576-01 
9.3556-01 
9.3566-01 
9.3606-01 
9 .3(66-01 
9.3756-01 
9.3866-01 
9.3596-01 
9.915E-01 
9.9336-01 
9.953fc~0t 
9.9756-01 

9. 9596-01 
9-5296-01 
9.5516-01 
9.5756-01 
9.(086-01 
9.6396-01 
9.6706-01 
9.7036-01 ' 
9.7356-01 
9.7(96-01 
9.CC2G-01 
9.8366-01 
9.8706-01 
9.9096-01 
9.9376-01 
9.9706-01 
l.OCOE 00 
1 . OC 36 00 
1.0C76 00 
1.010c 00 
1.0126 00 
1.0156 00 
1.0186 00 
1.0206 00 
1.02 36 Cl 0 
1,0256 00 
1.027E 00 
1.0296 00 
1.0316 00 


1.7.02E-01 
1.6506-01 
1 . (006-0 1 
1.5536-01 
1.9926-01 
1.92(6-01 

1.3556-01 
1.2806-01 
•1 .2C1E-01 
-1.1196-01 ' 

-1. 0336-01 
-9.9596-02 
-8'. 5526-02 
>7.6326-02 
>6. 6906-02 
-5.7556-02 *" 
-9.30(6-02 
-3.8556-02 
-2.9O6E-02 
-1.9696-02 
-1.0326-02 
-1.1306-03 “ 
7.8726-03 
1.6686-02 
2.5296-02 
3.3596-02 
9.1536-02 
9.92C6-02 ' 
5.6526-02 
6.3966-02 " 
7.0006-02 
7.6126-02 
8.1816-02 
8.705E-02 ' • 
9.1836-02 
9.6196-02 
9.9966-02 
1.0336-01 
1.0616-01 
1.C856-01 ~ 
1.1036-01 
1.1176-01 
1.12(6-01 
1.1296-01 
1. 1296-01 
1.1236-01 “ 

1.1136-01 
I.C906-OI 
1.0796-01 
1.0566-01 
t .C29G-01 
9 . 909E-OZ 
9.6976-02 
9.2736-02 
8.8(76-02 
8.9336-02 
7. 973E-02 

7.9096- 02 

6.9096- 02 
(.9616-02 
5.9216-02 



v>j W ^ W iJ UJ U> l hi 


3.37CE 

CO 

2 - 1 S 2 6 

0 1 

1.2106 

02 

1.0206 02 

3 - 3 ec E 

CC 

2. IS 36 

Cl 

1 .21CE 

C 2 

l . 0 2 C 6 0 2 

3 . 3 5 C 6 

CO 

2 . 1S3E 

Cl 

1 .2 ICE 

C 2 

1.C2G6 02 

3 • 4 y C E 
3.4 ICE 

CO 

2. IS 36 

01 

1.210E 

C 2 

1.0206 02 

CO 

2. IS 3E 

C 1 

1 .2 1 C6 

C 2 

1.0206 02 

3 . 4 ? C 5 

CC 

2 . 1 S 4 £ 

01 

1 . 2 1 CE 

C 2 

1.0206 02 

3 . 4 2 S t 

CO 

2. 1 5 4 £ 

Cl 

1 .2C56 

C 2 

1.0206 C 2 

3 . 4 ) S E 

CO 

2. IS 46 

C 1 

1 .2C5F 

C 2 

1 . C 2 0 6 0 2 

3 . 4 4 S 6 

CO 

2 . 1 5 5 E 

01 

1 . 2 C 9 E 

0 2 

1.C206 02 

3. 4 5 5 £ 

CO 

2.195c 

Cl ' 

1 .2C5E 

C 2 

l . C 2 0 6 C 2 

3 . 4 t 5 £ 

CO 

2. 1S5E 

Cl 

l .2C5E 

C 2 

1.C296 €2 

3.4 75 E 

CC 

2. 1S6E 

C 1 

1.2C5 6 

C 2 

l . C 2 C 6 C 2 

3 . 4 i> SC- 

CO 

2 . 1 9 6 E 

C 1 

1 , 2 C 5 6 

C 2 

1.C20I: 02 

3.4 4 5 E 

CG 

2 . ISfcE 

C l 

1 .2C5E 

C 2 

1 . C 2 0 G C 2 

3 . 5 C 5 E 

CO 

2. 1S6E 

C 1 

1 .2C5E 

C 2 

1.C2CC 02 

3 .5 1SE 

CC 

2. 1566 

C 1 

l . 2 C 5 E 

C2 

1 .0206 C 2 

3.5256 

CO 

2 . IS 76 

C 1 

1 . 2C P E 

C 2 

1.0206 02 

2 . 5 2 S E 

CO 

2. 1S7E 

Cl 

1 . 2 C e 6 

C 2 

1 . C 2 0 E C 2 

3 . c . 4 e E 

CC 

2. IS 7E 

C l 

l .2CEE 

C 2 

I.C2CE C 2 

3 . 5 5 i 6 

CC 

2 . 1 S 7 E 

C 1 

1.2CPE 

C 2 

1 . 0 2 C E 02 

3 . 5 L E 6 

CO 

2 . 1S7E 

Cl 

1 . 2 C 8 E 

C 2 

1.C20I: 0 2 

3.5766 

CO 

2. 1S7E 

Cl 

■ i - 2 c n E 

C 2 

1 . C 2 C E 02 

2 . 5 e 6 E 

CC 

2. 1S7E 

C 1 

1.2CP6 

C 2 

1 . C 2 C S 0 2 

3 .5566 

CO 

2. 1S?E 

C 1 

1.2C86 

C 2 

1 . 0 2 C 6 0 2 

3.(06 

CO 

2 . 1 5 7 E 

Cl 

1 . 2 C C t 

C 2 

1.C206 C 2 

3-6 166 

CC 

2 . 1 S 7 L 

C 1 

1.2CCE 

C 2 

1.C2C6 C2 

?.« 2EE 

CC 

2. 1S7E 

C 1 

1 . 2 C p 6 

C 2 

l . C 2 0 6 0 2 

3 . ft 3 6 E 

CO 

2.1576 

C 1 

1.2C8C 

C 2 

1.C206 02 

3.64 CE 

CO 

2 . 1 S 7 t 

C l 

l . 2 C fl C 

C 2 

1 . C 2 3 6 0 2 


3.(5Ec CO 
3. { t ?E CC 
3.617c CO 
3.6676 CO 

3 .c<;?eco 

3.7C7E CC 
3 . 7 l 7 E CO 
• 3.7/ 7 E CO 
3.737c CO 
3. 7*1 7 E CC 
3.757E CO 
3.7(76 CO 
3.7776 CO 
3.787E CO 
3.75(E CC 
3 . P C t E CO 
3.816E CO 
3.E26E CO 
3.E36E CO 

3.84tE CC 

.E5(E CC 
. 5 6 6 E CO 

• C74E CO 
. 5 e 6 E CC 
.SS£E CC 
.SC(6 CC 

• ShE CO 
.5256 CO 

3.9 256 CO 
3 . S ^ 5 6 CC 
3.9 556 ,CO 
3.9656 CO 


IS 7E 
19 76 
15 7 6 
l S 7 E 

IS 7 E 
1 S 7 f 
15 7 E 
1 5 7 £ 
1S7E 
IS 7E 
l S 7 E 
1S7E 
1S7E 
IS 7E 
IS 7 £ 
1 S 7 £ 
IS7E 
19 76 
2. 19 7 E 
2* 1S7E 
1 S 7 E 
IS 7E 
15 7 6 
1S7C 
19 76 
157E 
1 S 7 F. 
IS7E 
15 / C 
15 7 6 
1S7E 
1 S 7 E 


Cl 
Cl 
Cl 
Cl 
Cl 
C 1 
01 
C 1 
C 1 
C 1 
Cl 
C 1 
01 
C 1 
01 
Cl 
Cl 
Cl 
Cl 
C l 
C 1 
Cl 
01 
Cl 
Cl 
01 
Cl 
Cl 
C 1 
Cl 
Cl' 
Cl 


1 . 2 t P E C 2 
1.2CKE C2 
1.20 8 6 C 2 
1.2C86 02 
1.2C PC C 2 
1.20 PE C 2 
1.2CPE C 2 
1.2CPE C 2 
1 . 2 C e 6 c 2 

1 . 2 0 p C C 2 
l . 2 C fi E C 2 
1.2C6E C 2 
U2C86 C 2 
l . 2 C f? E C 2 
1 . 2 C ? 6 C 2 
i - 2 c e e c ? 

1.2CGE C 2 


l . C 2 C 6 C 2 


1 . C 2 0 l. 
I . C 2 0 £ 


0? 

02 


I.C206 02 
1.C20F 02 
1.C206 02 
1 .020 F 02 
1 . C 2 C C 0 2 
1 . C 2 C l: 0 2 
1 . C 2 C 6 0 2 
I.C206 0 2 
1 . C 2 0 f; C2 
1.C2CC C 2 
1 . C 2 C 6 0 2 
1 . 0 2 0 l C 2 
1 . C 2 C F. C 2 
I.C2CE C 2 


1 . 

8 70 6 

C 2 

- 1 -. 

8(4 6 - 

01 - 1 . 

1 7 4 r ; 

00 

1 . 

8 70 6 

C 2 

- 1 . 

6506- 

Cl - 

■ S . 

69CE- 

01 

1 . 

8 7 06 

C 2 

- 1 . 

9 2 06- 

■0 1 - 

-« . 

0326- 

01 

U 

8 706 

0 2 

- 1 . 

5 2 r - 

Cl - 

-7: . 

"’7 2 0 

0 1 

1 . 

8 7 0 6 

C 2 

- 1 . 

5 4 2 E - 

•Cl - 

■4 . 

"’196- 

0 1 

1. 

6 7 0 6 

C 2 

- 1 . 

s 4 1 r - 

•0 1 ■ 

.. ■> 

O 8 6 - 

C 1 

1 . 

8 7 0 6 

C 2 

-1 . 

5 2 2 6- 

■0 i - 

- 1 . 

'■ ' . 7 6 - 

0 1 

1 - 

6 / 0 \- 

0 2 

-1 . 

9 1 r >6~ 

■0 l 

1 . 


02 

l . 

P 70 E 

C 2 

- 1 . 

,89 06- 

-01 

1 . 

'’0 2 6- 

01 

1 . 

8706 

C 2 

- 1 . 


-0 1 

3 . 

1 ? nz - 

01 

1 . 

8 7 0 6 

C 2 

- 1 , 

12 19 6- 

-0 1 

/, . 

^ 111: - 

0 1 

1 . 

. 0 7 0 6 

02 

- 1 . 

.7 7 3 6- 

-0 1 

6. 

’ 4 56- 

•01 

1 . 

,P 7CE 

C 2 

-1 . 

.7210- 

-0 1 

7 . 

’*266- 

■0 1 

1 . 

,8706 

C 2 

- 1 . 

,(736- 

-0 1 

6 . 

, 4 f i 6 - 

-01 

1. 

.870 E 

C 2 

- 1 , 

.5506'- 

-0 1 

1 . 

.'‘in: 

00 

1 . 

. n v 0 C 

C 2 

-1 , 

. 5 2 9 1; ■ 

-0 1 

1 . 

.7 10 6 

00 

1 . 

.8706 

C 2 

-1 . 

.4556- 

-C 1 

1 . 

. *» 4 2 E 

00 

1 , 

. 870E 

C 2 

- 1 • 

. 3766- 

-C 1 

1. 

•* 4 f C 

00 

1 , 

. 0706 

C 2 

- 1 

.2926- 

-01 

1. 

'•456 

00 

l 

. P70E 

02 

-1 

. 2 C 5 E ■ 

-0 1 

1 . 

r ' 3 4 C 

00 

1 

. 8706 

C 2 

- 1 

. 1 156' 

-C l 

1 , 

'156 

00 

l 

. 8 7 0 E 

C 2 

- 1 

. C 2 2 6 

-0 1 

l. 

4 o nz 

CO 

1 

.8706 

C 2 

-9 

.2(46' 

-02 

1 

75 36 

00" 

1 

. 8706 

C 2 

-8 

.2676 

-0 2 

1 

8 09 6 

00 

1 

.8701: 

C 2 

-7 

.25 56 

-02 

1 

f S 6 f: 

00 

l 

. 0706 

C 2 

-6 

. 25 36 

-02 

1 

89 5 0 

00 

l 

. a v o r ; 

C 2 

- 5 

.2830 

- C 2 

l 

9? 56 

00 

l 

.8706 

C 2 

-4 

.2716 

-0 2 

l 

'i 4 7 6 

00 

1 

.8706 

C 2 

-3 

.2(06 

-C2 

1 

97,0 6 

OC 

1 

. 8 7 OG 

C 2 

- 2 

.25 ME 

- 0 ,7 

1 

9 (, 5 2 

00 

1 

. 8 70 L' 

C 2 

- 1 

- 2 f. 4 f 

- C 2 

l 

9 4 1 6 

00 

l 

. R 7 Cl (': 

C 2 

- 2 

. 0 C 5 0 

-03 

l 

9491. 

CO 

1 

.8 7 C 6 

C? 

4 

. H £ 4 F 

-0 1 

1 

n \ c 6 

CO 

1 

. 8 7 C L 

C 2 

l 

.(• 3 30 

-02 

1 

•io ; s,: 

0 0 

1 

.8706 

c 2 

2 

.567 !': 

-02 

1 

8 /> 8 L- 

CO 

1 

.8706 

0 2 

•» 

.4616 

-02 

l 

82 6 6 

CO 

1 

. 870 E 

C 2 

4 

. 3 t 0 

- 02 

1 

7 ! 8 6 

00 

1 

. 8 706 

C 2 


.12 3 0 

- C 2 

l 

7 2 3 6 

00 


-l 


1 .H 70 6 
1 . 3 70 6 

i . n v o c 

1 .8706 
1 .3706 
i . p 7 o i: 
1 .P.70 E 
l. 87CC 


C 2 

C 2 
C 2 
C 2 
02 
C 2 
C 2 
C 2 


.9 1 SI -02 
.6346-02 
.3486-02 
,55 7C-02 
. 5 9 9 r - 0 ? 


i .6.61 n 

1 .59 4 6 
1.5226 


00 

00 

00 


1.4456 00 


1.2CCE 

C 2 

1 .0206 

02 

1 . 8706 

02 

1 . 2C E E 

C 2 

1.0206 

02 

1.8706 

C 2 

1 . 2C8E 

C 2 

1-C20E 

0 2 

1.8706 

C 2 

U2C86 

C 2 

1 . C 2 C E 

02 

1.8/06 

02 

1.2CPE 

C 2 

1 »C 20E 

02 

1.3706 

C 2 

1 . 2 C 6 6 

C2 

l .C20C 

02 ' 

1.3706 

C 2 

1 . 2 C P 6 

C 2 

1 . 0 2 c L- 

02 

1.8706 

C 2 

1 . 2 C P E 

C 2 

1*0236 

0 2 

1.0/0 6 

C 2 

l . 2 C p E 

C 2 

1 .020E 

0 2 

1 .-8 70 6 

C 2 

1 . 2 C R E 

C 2 

1 . C 2 C E 

02 

1 .8706 

C 2 

l .268E 

C 2 

1.C20E 

02 

l.mce 

C 2 

1 . 2 C 8 6 

C 2 

L . C 2 C 6 

C 2 

1 . 0 7 ,)C 

02 

l . 2CHC 

C 2 

l . C 7 C E 

02 

1 .8706 

C 2 

1.2CEE 

C 2 

1 .C2CE 

C 2 

1.3706 

C 2 

1 *2006 

C 2 

1.C20C 

02 

1.8706 

02 


S . ! l L - C 2 

9.es'»e-c> 
1.C11F3-01 
1 .C * 16-0 1 
1.C666-01 
1 . n 5F.-C 1 

i . i s oe~c i 
i . iso r -oi 
1.1736-01 
1 . 1F2C-01 
l . U: 6 F - C 1 
1 . 1 6 9 C - 0 l 
1 .1 JOF-C 1 
1.U7E-CI 

i . r. ;; - c i 
1. 12 26 -Cl 
1.1C 75-01 


00 

00 


i . c 7 9 e - c i 


1 . 3036 
1.2 7/6 
1 . I 8 f! E C C 
1.C95E 00 

l.oone oo 

9.0326-01 
C.CA2E-01 
7.93C6-01 
6.076:6-01 
5.0106-01 
3 . T93F-01 
2-5016-01 
1.9776-01 
<i . 0556-0 2 
9 .9R2E-04 
- 9 . 6 6 / C - 0 2 
-l.i'7fit-01 
70 3 6- 01 
-3, (.596.-01 


1.8706-01 - 
1.902E-01 - 
1.926C-01 - 
1.941001 - 
1 .9486-01 - 
1 . 9 4 ( 6 - 0 1 - 
1.9371-01 - 
l .9 196-01 
1 . f! 9 6 £ - 0 1 
1 . B6 1 OO 1 
1 .0776-01 
1 .7 756-0 1 
■1.7736-01 
■1.6696-C-l 
1 .595001 
l .9296-01 
■1.454001 
■1.3750-01 
■1*2976-01 
■1.2066-01 
■1.11 4001 
C206-01 
-S. 2976-02 
-8. 2706-02 
-7.2776-02 
-6.2736-02 
-5.24 3002 
-9.2910-02 
-3.2390-02 
-2.2366-02 
-1 .29 1 OC2 
-2 . 5 POO -0 3 
7.C93E-03 
1.66 7 0 0 2 
2 . 66 1 F.-0 2 
3 . 6 7 96-02 
6.33 6 - 0 2 
5 . 1570-02 
5 .93911-0? 
6.6796-02 
7. 37 3 £-02 
C .C22F-02 

0. 6 24!" -02 

9 . 1 '766-02 
9 . 6 13 C 6-02 

1. C13O01 
1. 0536-01 
1 .008 6-0 l 
1.1186-01 
1.1426-01 
1 . 161E-01 

1 . 1 75E-01 
1 . IRSL-Ol 
1 .1006-0 1 
1.107 6-01 
1 .1016-01 
1 . 1 7CE-CI 
1.154001 
1*1346-01 
1.1096-01 
1 .0016-01 


1.1316 00 
9.643001 
7.9026-01 

-6.3166-01 
•4.69 76-01 
■3.0 1 16-C1 
■ 1 .-3 656-01 
2.1406-02 
1,7 Rll: - 01 
3. 3C9E-01 
4.7946-01 
6 .2296-0 l 
7 . 609 E - 0 1 
0 . 9 3 1 6 - 0 l 
1.019E 00 
1.1391: 00 
1.2506 CO 
1.355E 00 
1.4 52 E CO 
1.5 4 IF 00 
1.622C CC 
1.6956 00 
1.7596 00 
1.014E 00 
1.0626 CC 
1.900E 00 
1.9306 00 
1.9516 00 
1.9646 OC 
1.9 ( BE CO 
1.9646 00 
1.9536 00 
1 . 9 3 3 E CO 
1.9C5E CC 
1,0716 00 
1.8296 00 
1.7(306 CC 
1.7256 00 
1.66 36- 00 
1.5966 00 
1.5246 OC 
1.4466 QO 
1.3646 00 
1.2796 00 
1-109E 00 
1.09/6 CO 
1.0026 00 
9.042E-01 
0.0526-01 
7.0406-01 
6.0356-01 
5.0106-01 
4.CC2O01 
2.9696-01 
1.9056-01 
9.9256-02 
1*6726-03 

-9.3526-02 

-1.0726-01 

-2.7776-01 

-3.6536-01 


1.0336 00 
1.0346 00 
1.0356 00 
1.0376 OC 
1 . 0 3 n 6 00 
1.03 PC CO 
1.0396 00 
1.0396 00 
1.0406 00 
1.0406 00 
1.04 06 00 
1.0396 00 
1.0396 00 
1.0396 00 
1.0 J 6 6 CO 
1.0376 00 
1.0366 00 
1.0356 00 
1.0346 00 
1.0326 00 
1.0316 CO 
1.0296 00 
1.0206 00 
1.0266 00 
1.0246 00 
1.0226 00 
1.0206 00 

1.0 l Re 00 
1.0146 00 
1.0146 00 
L . 0 1 2 G 00 
1.0106 00 
1.0C06 CO 
1.0066 CO 
I.0C4E 00 
1.0C26 00 
9.9996-01 
9.9806-01 
9.9616-01 
9.9426-01 
9.9256-01 
9.9086-01 
9.8926-01 

9 . 0 746-01 
9.86 16-01 
9.0486-01 
9.8356-01 
9.823E-01 
9.8126-01 
9.0C2001 

. 9.7936-01 
9.7856-01 
9.7/06-Cl 
9. 7726-01 
9.7686-01 
9.7646-01 
9. 7616-01 
9. 7606-01 
9. 759E-01 
9.7606-01 
9. 7616-01 


5. 368002 
4 .804E-02 
4.2316-02 
3.4536-02 
3.0716-02 
2.4886-02 
1 . 9 0 1 E - 0 2 
1 . 3 30E-02 
7 . 5 BSE- 0 3 
1 . 96 0 E - 0 3 
-3.5596-03 
-8.9516-03 

- 1 . 4 1 9 E - 0 2 

- W92 76-02 
-2.4176-02 
-2.8066-02 
-3.3346-02 
-3.7596-02 
-4. 15SE-02 
-4.5346-02 
-4. 0826-02 
-5.203E-02 
-5.495E-02 
-5.7596-02 
-5.9936-02 ■ 
-6.1976-02 
-6. 3 70 E -02 
-6.5146-02 
-6.6266-02 
-6. 7096-02 
-6.7626-02 
-6.7856-02 
-6. 7/86-02 
-6. 7 4 3 E - 0 2 
-6.f.80£-02 
-6.5906-02 
-6.4736-02 
-6. 332E-02 
-6 . 1 7.66-02 
-S. 9776-02 
-5.7676-02 
-5.5376-02 
-5.2886-02 
-5-C226-02 
-4.7396-02 
-4.4436-02 
-4. I 34E-02 
-3.8 136-02 
-3.403E-O2 
- 3 . 1 4 4 E - 0 2 
-2.800E-02 
-2.4506-02 
-2.0966-02 
-l . 7416-02 
-1 . 3856-02 
-1.C3 16-02 
-6.7896-03 

3. 1096-03 
1 . 1766-04 
3.4786-03 
6. 7596-03 



3. 97 56 00 
2 - r i n t. E CC 
3.SS5E CO 
A.CCSE CO 
A.C156 CO 
A.C25G CO 
A.CJA6 cc 
A . C AA t CO 
A . C 5 A 6 CO 
A . C 6 * 6 CO 
A * C /At C C 

A. CRAG cc 

a . c 9 a e co 
A . I C A G CO 
A. MAC CO 
A . 1 2 A 6 CO 
a . 1 3 A 6 CC 
A . 1 A A C C C 


l 


CC 


A . 1 1 1 6 CO 
A.172C CC 
A . t £ 2 £ C C 
A . 1 9 3 li c c 
A . 2 C 2 £ CO 
A . 2 L 2 C- CO 
A.222E CC 
A. 2 3 36 CC 
A.2A2E CC 
A . 2 5 2 £ C C 
A . 2 t 3 E CC 
A . ? 7 .1 6 CO 
A.2S2E CO 
A . 2 *5 2 E CC 
A. 2 CPC- CC 
A . 2 l 2 E CO 
A. 3226 CC 
A . 2 2 2 E CC 
A.3A2G CC 
A . 3 3 2 E CO 
A . 2 62 E CC 
A. 2 7 2 6 CC 
A.3E2E CC 
A*292E CO 
A . A C l E C C 
A. A 1 IE CO 
A . A 2 1 6 CO 
A . A 2 l E CO 
A . A A 1 E CC 
A . A 5 1 E CC 
A . A 6 1 E CO 
A .A 7 1 E CO 
A . A 8 1 £ CC 

a. as i e cc 

A.5C16 CO 
A.511E CO 
A . 6 2 C E CC 
A . 6 3 C E C C 
A. 5 ACE CO 
A.55CE CO 
A . 5 tC 6 00 
A . 5 7 C G CO 


2 . 
?. . 
2 * 
2 . 
2 . 
? . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 


1S7G 

is n. 

1 S 7 E 
1 S 7 C 
15 7 0 
5 S 7 G 
1 S 7 G 
1 S 7 E 
1 A 7 E 
i s 7 e 
, isyr 
. 1 A 7 E 
. 1 A 7 F_ 

, l S 7 E 
. 1S7G 
. 1 5 7 E 
.19 71; 

. i s 7 r 
. i s 7 c 
. IS 7E 

. twe 

2 . 1S7E 
2. 1S7E 
1 0 7 E 
IS 7 6 
IS 7 E 
15 7 6 
2. IS 76 
2 . I S 7 E 
2 . 1 S 7 E 
2 . 1 5 7 6 
2 . I S 7 E 
2 . 1 5 7 E 
2 . 1 S 7 e 
15 7 G 
1 S 7 E 
1 S 7 E 
1 S 7 E 
1 S 7 E 
1 S 7 E 
1S76 
1S7E 
1S7G 
1 S 7 E 
1 S 7 E 
, 1 S 7 E 
. 1S7G 
. 1 S 7 E 
, 1 S7 E 
. 1 S 7 6 
. IS 7E 
. 1S7S 
1S7E 
1 S 7E 
IS 7E 
IS 7G 
IS 7E 
187 6 
1 S 7 E 
1S7G 
1 S 7 C 


2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2. 
2 . 
2 
2 
2 . 
2 
2 
2 


Cl 
C 1 
C l 
C 1 
Cl 
C 1 
01 
C 1 
C l 
Cl 
C 1 
C 1 
C 1 
Cl 
C 1 

c 1 "" 

Cl 
C 1 
Cl 
C 1 

C 1 

Cl - 
Cl 
Cl 
Q 1 
Cl 

c l 

Cl "■ 
Cl 

c l 
0 1 
Cl 
C l 
Cl 
01 
Cl 
Cl 
C 1 
C l 
01 • 
C 1 
01 
01 
Cl 
01 
Cl 
Cl 

c 1 

Cl 
Cl 
C l 
Cl 
C 1 
0 1 
Cl 
Cl 
C l 
C 1 
01 
Cl 
C l 


1.2C8E C 2 
1.2CPE C 2 
1 . 20 EE C2 
1 . 2 C P E C 2 
1.2 CEE 02 
1.2CPC C 2 
l.2C(’r. C 2 
l . 2 . c r G C 2 
1 . 2 C f! E C 2 
1.2CCC C 2 
1.2C6P C 2 
1.2CP6 02 
1.2CFE C 2 
1.2 CHE C 2 
1.2CSE 0 2 
"1.2CI2E C2 “ 
1.2C8E C 2 
1.20 EE C? 
1.2CCC C 2 
1 . 2 C <’ E C 2 
1.2C2E C 2 
“ 1 .2C6E C 2 
1.2CCE C 2 
1.2CEE C 2 
1.2C8E C2 
1.2CPE C 2 
1.2CG6 C 2 
- 1 . 2 C E 0 02 
1 . ?C EE C 2 
1.2CRE C 2 
1 . 2 C P G 0 2 
1.2C0E C 2 
l . 2C2 6 C 2 
l . ?Ci?E C 2 
1.2008 C 2 
1.2C86 C 2 
1.2CPE C 2 
1.2CPE C 2 
1.2C8E C 2 
-1.2CCE C 2 
1.2C8E 02 
1.2CCE C 2 
1 . 2C C E C2 
1.2C06 C 2 
1.2CEE C 2 
1 .2C8E C 2 
1.2C8E C 2 

l.aceG C2 

1 .2CRE 02 
1.2CRE C 2 
1.2C«E C 2 
1.2CCE C 2 
1.2CRE C 2 
1.2086 02 
1.2CPE C2 
1 . 20 EE 0 2 
1.20 RE C 2 
1.2C0C C 2 
1.2C86 C 2 
1 . 2 C 0 E C 2 
1.2C86 C 2 


I.C20E 02 
I.C20E C? 

I .0 2 CC 0 2 
1 . C 2 0 6 0 2 
UC20C 02 
1 . C 2 0 1: C 2 
I.C20E C ? 

1 . C ? 0 C C 2 
1.C20E C 2 
1.C200 02 
1.C20E 02 
1 . 0 2 C. t 0 2 
1 . C 2 0 1- C 2 
1.C20E 02 
1.C20E 02 
' 1 . C 2 0 E 0 2’' 
1.C2CG 0 2 
1.C2CE 02 
1 . 0 2 0 C C 2 
1 . C ? 0 E 0 2 
1 . C 2CE C 2 
' 1.C29E 0 2 
1.C20G 02 
1.C20E C 2 
I.C2CG C 2 
1.0 2 CE 02 
1 .020!: 0 ? 

1 . C 20 E C 2 
I.C20G 02 
1.C20E 0? 
1.020E 02 
I.C2CP C 2 
1.02CC 02 
• 1.C20C 02 
l . C 2 C G 0 2 
] .C2CE C 2 
1.0200 02 
1 . C 2 0 1: 02 
1.C20F 0 2. 
’1.0206 C2 
1 . C 2 C E 0 2 
1.C2CE 02 
1.020E 02 
1.C20E 02 
1.C2CE 02 
1.02CE 02 
1.02CE C 2 
1.C2CG C 2 
1.C2CG C 2 

1 . G 2 o e 0 2 
1.020E 02 
l . a 2 C E 0 2 
1.02CE 02 
1.020E 02 
1.C20E 02 
I.C2Q0 02 
1.02QG 02 
1 . C ? C E 0 3 

1 . o ? 0 c 0 2 

1 . 0 2 0 E 0 2 
1.02 Oil 0? 


1.R70G C 2 
1 . ft 7 0 P C 2 
1.R70F 02 
1.07 0 I: C ? 

l.rnoi: c 2 

1 . 12 7 0 P C 2 
l . f» 7 0 C C 2 
1 . P 7 0 E C 2 
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1.0086 00 
l « 00 5 1: 00 
1.0056 CO 
1.0C56 CO 
1.0C5G CO 
1 .0056 00 
1.0056 GO 
1 .0056 00 
1.0C5E 00 
1.0C56 00 
1 . 0 C 5 6 00 
1 . CC 5 £ CC 
1.CC9E CO 
1.0C3E 00 
1.0C3E 00 
1.0036 00 
1.0C3E CO 
1.0036 00 
1.0C7E 00 
1.0C2E 00 
1.0C2E 00 
1.002E 00 

i.ccie oo 
l.ocie oo 
1.0016 oo 

1.0C1E 00 
1.CCC6 CC 

l.ocoe oo 

9 . 99 BE - 0 1 
9.9866-01 
9.9856-01 
9 .952E-01 
8.990E-01 
9.9096-01 
9. 9076-01 
9.9056-01 
9.9856-01 
9 .9836-01 
9.9816-01 
9 .9006-0 1 
9, 9796-01 
9 . 9 7 9 C - 0 1 


1.^.086-02 
1 . 372E-02 
1 . 3 376-02 
1 .2076-02 
1.23HE-02 
1 . 1056-02 
1.129E-0 2 
l .0676-02 
l .CO 7E-02 
9. 5 29 E-0 3 
fl. 7636-03 
n.CfiOE-03 
7.3056-03 
6.6776-03 
5.862E-03 
5.253E-03 
5.5736-03 
3.R06E-03 
3.C95E-03 
2 . 309E-O3 
1.69 5 E-03 
1.0156-03 
3. 5086-05 
-2.952E-05 

- 9 . 7 l 5 E - 0 5 

-1 .575E-03 
-2. IC5E-C3 
-2.6556-03 
-3.1796-03 
-3.675E-03 
-5. 138E-03 
-5.5706-03 
-5.9676-03 
-5.J25E-C3 
-5. 6556-03 
-5.8556-03 
-6. 199E-03 
-6.51SE-03 
-6.602E-03 

— 6 . 7 5 2 E-0 3 
-6.8666-03 
-6. 952E-03 
-6.9016-03 
-6.9016-03 
-6.9556-03 
-6. 67CE-C3 
-6.759E-03 
-6.612E-03 
-6.5326-03 
-6.2226-C3 
-*.9056-03 
-5.7226-03 
-5.537E-03 
-5.1336-03 
-5. BIO E-03 
-5.571 E-03 
-5 . 1196-03 
-3.7566-03 
-3.3836-03 
-3.C02E-03 
-2. 6166-03 


102 ! 



* . 7 5 C E CO 2 . 1 8 7 £ Cl 1 . 2 C e £ C2 1 . 0 2 0 C 02 l . 8 7 C E 02 1.3C4U-02 

5.eCCE CO 2 . 1 5 7 E 01 1.2CPE C2 I.C2CE 02 I.H70G C2 1.3C0E-02 

5 - e i C E CO 2 - 1 ^ 7 E Cl 1 . 2 C 8 E C2 1 . C 2 C C C2 1 .8701: C2 1 . 3 C 5 E - 0 2 

5 . £ ?tC E CO 2.1876 Cl 1.2CEE C2 1.C20E 02 l.OTQf? C2 1.2576-02 

5.^3CE CC 2 * 1 C 7 E Cl 1.ZC86 C 2 1.020E 02 1.8706 C2 1.264C-02 

5 „ £ 4C 6 CO 2 - 1 ^ 7 G Cl 1.2CFC C2 1.C20E 02 I-U70E C2 I - 2 fi & f*— 02 

5. CSCE CO 2.1876 Cl 1 . 2 C ? E C2 1 . C 2 C 6: C2 1.0 70C C? 1 . 2 '• 2 1- - 0 2 - 

' 5 . E 5 8 6 CC 2.1876 Cl 1 . 2 C 8 5 C2 1 . C 2 0 E 02 1.P7 0C C2 l . 2 1 .4 r > 0 2 - 

5 . 0 t ^ E CO 2 . 1 5 7 i 01 1 . ? C 0 E 0 2 l . 0 2 0 E 0 2 1 . " 7 0 6 0 2 1 . 1 E l 6 ~ O 2 - 

5 . E 7 S E CC 2 . I S 7 E Cl 1.2CRC C 2 1.C206 0 2 1 - 7 C 6 C? 1.14 4f“0; ! - 

5. asst CC 2.1576 CL 1.2CCC C2 1.C2C6 02 I.H70E C 2 1. I C 30-02 - 

5. WE CG 2.1576 CL 1.20 PC C2 1.C20E 0 2 1.P70E C? l.C6ni'-G? - 

5 . r; C 9 6 CO 2.1576 Cl l . 2 C C E C 2 1 . C 2 C f: 02 1.0701: C2 l.CC'H-C? - 

5 . s i q e c c 2 . r; 7 f ci i . 2 c f n c 2 1 . 0 2 0 e c 2 1 - n v 0 (■ 02 «; . 0 7 4 1 - c j 

CC 2.1876 Cl 1.2CP6 C2 1.C2CC 02 1.0 700 C2 9.C27r-03 

5.4! 8 [• CO ' 2.187E Cl 1.2CSC C2 1.C2CE02 1-8796 C 2 0.451C-G3 

5.9496 CO 2-157E Cl 1.2CPE C2 1.C2CE 02 1 . f! / 0 Cf C? 7. PE 6 C-OZ 

' 5.WE CC 2.1576 Cl 1.2 CP 5 C2 1.C2CH C2 1.0 7 0C C? 7.2 2110-0 3 

5.9 £56 CO 2.1S7G Cl . 1.2C8E C2 1.C20E C2 1-8706 C2 A.tC?H-03 

5 . 7 e E CO 2.187E Cl 1.2C6E C2 1.C2CE C2 1-8 7CE C2 5.951E-03 

5 . ^ E 0 E CO 2 - 1 5 7 t 0 1 1.2CPE C2 1.C2CE 02 1.8736 C2 5 . 2 E 0 E - C 3 

— -S.CCEE CO 2. 15 76 C1~"1-2C8£ C2"’1.C2CE 02" 1.8 70 E C2 4.616E-02 

E.CCEE CO 2-187E Cl 1.2CBE C 2 1.C2CE 02 1.87CE C2 3-8396-03 

6. C1EE CC 2.187E Cl 1.2CSE C2 ' 1.C2CE 02 1.P70E C2 3.269G-C3 


6 . C 2 £ c 

CO 

2.187c 

Cl 

1.2CPE 

C 2 

1 . C 2 0 E 

02 

1.P7CC- 

C? 2.573E-C3 

6.C 32E 

CO 

2 . 18 7 £ 

Cl 

1 . 2CPE 

C 2 

I.C2C-E 

C 2 

UP70E 

C 2 1.8CCE-C3 

6 . C 4 £ £ 

CO 

2.1876 

Cl 

l . ?CPE 

C 2 

1.C2CE 

G 2 

1 .G70E 

C 2 1.22 9E-C3 

~ 6.C5EE 

CO " 

2. 187E 

Cl '■ 

' 1.2C8E 

C 2 

' 1 . C 20 5 

02 

1 . fi 7 0 E 

C 2 5 .0276-04 

. C t S E 

CO 

2.157c 

Cl 

1.2CF6 

C 2 

1 -C2C6 

0 2 

1 .0706 

C 2 -9.1036*05 

6.C 7c c 

CC 

2.1876 

Cl 

1 . 2 c e E 

C 2 

l .C20G 

C 2 

1.07CT 

C 2 -7.34 36-04 

t . C E £ £ 

CC 

2 . 1 8 7 £ 

Cl 

1 .2CPE 

C? 

I.C20E 

C 2 

1 - 8 7 0 E 

C 2 -1.3 £26-0 3 

6 . C 8 7 £ 

CC 

2.18 76 

C l 

l . 2 C ? 16 

C 2 

1 . G 2 Q E 

0? 

1.8 7011 

C 2 -1 .87 3E-C3 

6 . 1 C 7 E 

CO 

2 . 1 5 7 E 

C 1 

t . 2 C P 6 

C 2 

1 . C 2 G t: 

C? 

1 . n v 0 f 

C 2 -2.5546- C 3 

‘ 6.1176 

CC 

2 . 1 9 7 E 

C 1 

' 1.2C PE 

C 2 

1 .C2C6 

C 2 

1.9 70 E 

C 2 -3.0636-03 

.6.1276 

CC 

2 . 1 8 7 i: 

C 1 

1.2CPE 

C 2 

1 . C 2 9 E 

0? 

1 . 5 706 

0 2 -3.5£4r-c:: 

6.1376 

CC 

2.1576 

c 1 

1 .2CP6 

C? 

1 .C20C 

02 

l .8706 

C 2 -4 .CC 26 -03 

6 . 1 4 7 E 

CO 

2.157E 

Cl 

1 . 2C P E 

C 2 

1 . 0 2 0 1; 

0? 

1 . 8 7 0 E 

C 2 -4.4046-02 

6.157b 

CO 

2 . 1 8 7 E 

Cl 

1 . 2 C P E 

C 2 

1.C2CE 

C? 

1 .07CE 

C 2 - 4 . 7 7 3 £ - C 3 

6.1676 

CC 

2 » 1 8 7 E 

Cl 

1.2CPE 

C? 

1 . C ? 0 E 

02 

1 . C 7 0 P 

C 2 -5. 1136 -03 

•’ 6.1776 

CO • 

2.1876 

Cl 

1 . 2 C P C 

C 2 

1 . C 20 6 

02 

1 .P.7CE 

C 2 -5.4 296-03 

6 . i 8 7 £ 

CC 

2.187c 

Cl 

l .2CPE 

C 2 

1 . C 2 0 E 

02 

1 . P 7CE 

C 2 -5,7236-0? 

,6.1876 

CC 

2.1876 

Cl 

1.2CPE 

C 2 

l.C/CC 

0? 

l .8706 

C 2 -5.980C-C3 

6.2C 76 

CC 

2.1876 

Cl 

1.2CE6 

C 2 

1 . 0 ? C C 

C 2 

l . 0 7 0 E 

C 2 -6.2426-0:1 

6.216c 

CO 

2 . 15 7E 

C 1 

1 . ZG ? E 

C 2 

1.C20E 

0 2 

1 .870E 

C 2 -6.25 46-C! 

6 . 2 2 1 E 

CC 

2.1876 

Cl 

1 .2C6E 

C 2 

l .C206 

02 

1.0706 

C 2 -A. 4 £36-0! 

'* 6.236E 

CO ■ 

2 . 1 8 7 G 

c l 

" 1 . 2 C P 6 

C 2 

" 1 . C 2 C E 

02 ' 

’ 1.0 7CE 

C 2 -6.4656-03 

6.246 c 

CC 

2. 15 7C 

c 1 

1 . 2 C P 6 

C 2 

1 .C 2 C 6 

C? 

1 8 70 6 

C 2 -6.4 1 16-03 

6 . 2 5 6 E 

CO 

2 . 1 5 7 E 

c 1 

l . 2 C fi E 

C 2 

1 .C20E 

C 2 

L . 0 70 E 

02 -6.3156-02 

6 . < fc 6 6 

CO 

2.1576 

c 1 

l .2C E6 

C 2 

1 .C206 

02 

I .8706 

C 2 -6.1PAC-0 7 

6.2 7 6 6 

CC 

2.157H 

Cl 

1 . 2 ce e 

C 2 

1.C2CC 

02 

1 . 0 7 C C 

C 2 -6. 0216-03 

6 . 2 £ 6 E 

CC 

2 . 18 7E 

Cl 

1.2CEE 

C 2 

1 . C 2 0 C 

02 

1 . P 7 0 6 

C 2 -5.8576-03 

' 6 . 256 6 

CO 

2 ♦ 1 8 7 £ 

0 1 

'1 .2C?E 

C 2 

1 . C 2 0 5 

0? ' 

" 1 .R70E 

C 2 -9.6706-0 1 

6. 2C6 6 

CC 

2.1576 

Cl 

1 . 2 C P E 

C 2 

1 . 0 2 0 6 

02 

1 . 8 l C C 

C 2 -9. 4 7 56- C 2 

t.21tc 

CC 

2 . 1 8 7 E 

c 1 

1 . 2 C » c 

C? 

1 . C 2CE 

C? 

1.8706 

C 2 - 8 . .7 7 4 C - C i 

6.2266 

CC 

2. 18 7E 

Cl 

1 .2Cpr : 

C? 

1 . C 2 0 C- 

0? 

1 .8706 

C 2 - 5 . 0 7 3 £ - C 3 

6 . 2 i 5 E 

CO 

2.1576 

Cl 

1 .2CPE 

C 2 

1 .C2CE 

02 

1.8701: 

C 2 -4.8726-02 

6.3456 

CC 

2 . 1 8 7 C 

c l 

l . 2 c e 6 

C 2 

1 . C 2CG 

02 

1.87 C C 

C 2 -4.6 746-03 

' 6 . 2 556 

CC 

2. 1876 

Cl 

l . 2 C C E 

C 2 

l .C?OE 

02 

i . P. 7 0 C 

0 2 -4.4UCr.-C3 

6.2656 

CC 

2 . 1 8 7 E 

c 1 

1 . 2CPE 

C 2 

1 . C 2 C C 

02 

1 .870E 

C? -4.28? -C i 

6.3756 

CO 

2.1576 

Cl 

1.2C 66 

C 2 

1 . C 2 C Ei 

02 

1,0706 

C2 -4.1 116-03 

6 . 3 £ 5 E 

CO 

2 . 1 8 7 E 

Cl 

l . 2 C P E 

C 2 

I.C20E 

02 

1,8706 

C 2 - 3 . 9 2 7 E - C 2 


(> . A l' 4 E - 0 2 1.306E-0?. 6.6CSE-02 9.9786-01 -2.226E-03 

5.525E-02 1.310E-02 5.526E-0 2 9.9776-01 -1.034E-O3 

4 . 3 7 517- 02 1.30 7 0- 02 4.378C-C2 9. 9 7 70-01 -1.4426-03 

9.2416-02 1.30CE-02 3.2436-02 9. 976E-0 1 -1.952E-03 

7 . 1 2 ( !* - 0 2 1 . 2 H () f: - 0 2 2. 120 £-02 9.97 61-01 - 6 • 6 4 3 E - 0 6 

1.0161-02 1-2686-0? 1.0376-02 9.97ftf>01 -2.ei76-04 

3 . C 4 2 6 - 0 4 1. 24413-02 7016-04 9.97AI-CI 9.4A5C-0S 

1.06. IE -0 2 1 .21 6 £-02 -1-060E-02 9.9746-01 4 . 4 9 2 G - 0 4 ( 

2. Of, U-02 l.lP.3r:-02 -2. 087L-02 9 . 977 , 0.-01 8.2256-04 

■> . r , I f - 0 ? 1 . 1 6 (: E - 0 2 - 3 - 0 1 7 r - C 2 9 . 9 ( 6 1: - C 1 1 . 1 7 1 £ - 0 3 

W3GC-02 1. 1046-02 - 3 . 9 3 E - 0 2 9. 97 76-Cl 1.598(3-03 

4 . 19 6 - 0 2 1 .0 55 13-0? -4. BOH 6-02 9. 9770-01 l . R 3 3 C -0 3 

5 . ' 4 C E - 0 2 1.011E-02 -5.6356-02 9.97 0U-C1 2. I 4 36-03 

(• . '■ 1 7 E - 0 2 9 . 9 ft 9 (3 - 0 3 - 6 , A 1 3 G - C 2 9.97HL-C1 2 . 6 3 8 E - 0 3 

W 6 6 12-0 2 9.06 2E-C3 -7.139E-C2 9.979G-01 2.71/6-03 

■ 7 . ° 1 7 f; - 0 2 0*7,676-03 -7.P12C-02 9 . 9 POL -01 2.979E-03 

P . 6 3 6 - 0 2 7 .1169 E- 0 3 -0.6 3LE-02 9. 9 PIC -01 3.226E-03 

0 . 9 ,J 7 E - 0 2 7.25CE-03 -P. 993E-C2 9W2C-01 3* / *51E-03 

9 .**> 0 3 £ - 0 2 6.613E-03 -9.690C-02 9.9l!3E:-UI 3.6 r jCE-03 

9 . 9 *> 0 H - 0 2 8. 96 IE -0 3 -9 . 966E-02 9 . 9 P 6 1. -01 3.367E-03 

1. C36E-01 5.2976-03 - 1.033E-C1 9.9C1.E-01 4.0156-03 

1.C67E-01 4.6266-03 -1.066E-C1 9.936E-01 6 . 1666-03 

-l -09 46-01 3.9456-03 -1.094E-C1 9.9C7E-01 4.2926-03 

-1.115E-01 3.264L-03 -L. 1156-01 9.9896-01 4.400E-03 

-1.13 1E-01 2.502E-C3 -1. 130 6-C 1 9.980E-01 4.4 876-03 

■ 1 . 14CE-01 1.903E-0 3 -1.14 06-01 9.99 LE-0L 4.554E-03 

1 - 144E-0 1 1.2 3 16-03 -1.144 6-01 9.99 36-01 4.6016-03 

1.1436-01 5 .64 66-04 -1.1426-01 9.9946-01 4.62PE-03' 

1.1366-01 -9.1P7E-05 - 1 . 13 5 6-01 9.995E-C 1 4.6 356-03 

1 . 12 36-0 1 -7. 3566-04 -1 . 1 236-0 1 9.99 76-01 4.62 3E-C3 

1.107:6-01 -1.365E-03 -1.106E-C1 9. 9986-01 A.591E-03 

1.0846-01 -1 .9 766-03 -1.0046-0 1 l- 00 06 00 4.54 IE -0 3 

1.9466-01 -2.8506-03 -1.056G-C1 l.CCOE CO 4.4736-03 
1.0246-01 -3.0876-03 -1.C246-01 l.OCOE 00 4.3986-03 

5 .4786-02 - 3 . 5 6 ?’ E - C 3 - 9 . 0 / 7 E - C 2 l.CCOE 00 4 . 2 8 P E - 0 3 

9.4796-02 -4.007E-03 -9.4746-02 1.0006 00 4.1746-03 

8 .0406-02 -4.4096-03 -9.0396-02 l.OCOE CO 4.0486-03 

P.5756-02 -4.7796-03 -0.575E-C2 1.0C1C 00 3.910E-03 

P . r HOC - 0 2 -5 . 12CE-03 -0.0«^E-02 1.001L 00 3. 7625: -0 3 

7. '.756-02 -5.436E-03 -7.5756-02 1.0C1C 03 3.6046-03 

7.C446-02 -5.731E-C3 -7.0466-02 1.0C1E OO 3.436E-03 

6.4986-02 -6. CO 6 6 -03 -6.499E-C2 l.OClE 00 3.2616-03 

5.9346-02 -6. 2506-03 -5.9366-02 1.0016 00 3-0776-03 

5.3496-0? -6.4026-03 -5.351C-C2 1.CC16 00 2.8G7E-03 

4.7486-02 -6.4 706-03 -4.7526-02 l.OCLE CO 2.694 E -03 

■4. 1466-02 -6.4716-03 -4.1496-02 l.OClE 00 2.500E-03 

3.5466-02 -7j.4K.6-03 -3. 8496-02 1.001C 00 2.J07E-03 

2.9566-02 -6. 322E-03 - 2. 9 59 6 -C 2 l.OClE CO 2.116E-03 

■2. 3 7 RE-02 -6.1526-03 -2- 3826-02 l.OClE 00 1.920E-03 

■i . P I 9 6 - 0 2 -7j.037E-C3 -1.823E-02 1.001E 00 1. 745E-03 

•1.2796-02 -5.0636-03 -1.2036-02 l.OOlfc 00 1.5666-03 

•7.7 046-03. -5.6766-0 3 -7.6576-03 l.OClE CO 1.393E-03 

-2.657E-03 -5 - 4 1! C 6 - C 3 -2.71IE-C3 1-OOU. 00 1.226E-03 

2-C5IE-03 -5.20CC-03 1.9556-03 1.001C CO L. 065E-03 

6 . r L 7 K - U 3 -8.0706-03 6.4506-03 1 . 00 1 E 00 9. 0966-04 

l. "746-02 -4.0776-03 1.060E-02 l.OClE 00 7.6C4E-04 *‘‘K« 

1.4 726-02 - 4 . 6 7 9 6 - C 3 1.4656-02 l.OClE 00 6.1726-04 

1. °4 66 -C 2 -4.4056-03 1.8406-02 1.0C1C CO 4.7996-04 

2 . ’976-0? -4.2976-03 2.1916-02 l.OOLE CO 3.403G-04 

2.^246-0? -4.116C-03 2.5196-02 1.001E 00 2.2226-04 

2. ’'3 36-02 -3.9426-03 2.827E-02 l.OClE 00 1.C18E-04 



6.2556 CO 
6.AC5C CO 
t.Uit CO 
6.A25E CO 
t.V’bf CO 
t . A A 9 CO 
. 4 5 5 6 CO 
t.AoAE CO 
0.4749 CO 
fc . 4 B 4 £ CC 
6.4549 CO 
6.5CAE CO 
fc . 5 1 4 E CC 
6.5246 CO 
6.5346 CO 
6.54A6 CO 
6.55AE CO 
6.5649 CC 
6.5746 CO 
6.5026 CO 
fc . 50 3 E CO 
6 - 6 C 2 E CO 
6. £139 CO 
6.6 2 2 E CO 
6.632c CC 
6 . 6 A 3 E CC 
6.653E CO 
6 . 6 6 3 E CC 

6 . 6 7 2 E CO 
fc.fcWE CC 
6 .6-5 2 E CO 
6 . 7C2 £ CO 

6. 7 126 CC 
6.7 2 2 6 CO ' 
6.7226 CO 
6.7426 CO. 
6.752E CC 

6 . 7 6 2 E CO 
6.772E CO 
6 . 7 6 2 E CO " 
6.7526 CC 
ft . 9 C 2 E CO 
'6.0126 CO 
6.E21E CO 
6.8316 CC 

- 6-gAie co 

6.8516 CC 
• 6.E61E CC 
6.6716 CC 
6.6816 CO 
ft . PS 1 6 CC 
~ 6.5C1E CC 
6.0116 CO 
fc. <5216 .CO 
6.5316 CC 
ft , 5 A l 6 CO 
6 . 0 5 C £ CO 

— 6.56CE CC 
6.57CE CO 
6.59CE CO 
6 - 0 <1 C 6 CC 


1579 
I S 7 E 
1C I E 

1 WE 
15 7 E 
15 7E 
15 7 6 

1 we 

2 . 1 W E 

2.1579 
15 7 6 
15 7 6 

1 WE 
15 76 
15 7E 
1WC 
L 5 7 6 

1 5 7 E 
1 5 7 E 
15 7 9 

. 15 7E 
.15 7 6 
. 15 7E 
.1WE 
. 1576 
. 1 5 7 £ 
.1576 
.15 7 6 
.15 7 8 
. 1 5 7 £ 
15 76 
15 76 
2.1576 
2.157E 
15 76 
1 5 7 E 
1 5 7 E 
1576 
1576 
1 5 7 E 
1576 

2 . 1 5 7 E 
2 . 1 5 7 E 
2. 1 5 7 E 
2 . 1 5 7 E 
2.1579 
2 . 1 5 7 E 
2.1576 
2.1576 

2 . 1 5 7 E 
2 . 1576 
2. 1579 
2.1576 
2.1579 
2.1576 
2.1576 
2. 15 7E 
2.157E 
2.157E 
2 . 1 5 7 E 

2 . 1 5 7 G 


2 . 

2 
2 . 

2 

2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 . 
2 
2 
2 
2 


01 
C l 
C 1 
Cl 
C 1 
C 1 
0 l 
C 1 
C l 
Cl 
0 l 
C 1 
C 1 
Cl 
C l 
Cl ' 
C 1 
0 1 
C 1 
Cl 
C 1 
C 1 
Cl 
Cl 
Cl 
Cl 
Cl 
Cl" 
Cl 
C 1 
Cl 
C L 
Cl 
01" 
Cl 
Cl 
0 1 
Cl 
Cl 
Cl ' 
C 1 
Cl 
C 1 
C 1 
Cl 
01 
Cl 
Cl 
01 
C 1 
Cl 
C 1 
Cl 
C 1 
Cl 
Cl 
Cl 
Cl 
C 1 
Cl 
Cl 


1.2C8E C2 
1.2CPE C ?. 

1 . 2 C « t C 2 
1.2CTE C 2 
1.2CPE C 2 
1.2CrE C 2 
1.2C8C C 2 
1 .ZCPfc C 2 
1.2CEC C 2 
1.2CPE C 2 
l.ZCPfc C 2 
l . 2 C P b C 2 
1 . 2 C !: E C 2 
U2C0C C 2 
1 . 2C 6 6 C 2 
1.2CCE C 2 
1.2C8C C 2 
1.208E C 2 
1.2C86 C 2 
1 . 2 C e E C 2 
1.2C8E C 2 
1 . 2 C 8 r C 2 
1.2C6E C 2 
1 . 2 C £ 6 C2 
1.2C6H C 2 
1.2C6E C 2 
1.2C8E C2 
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